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PULSED CYCLE OPERATION OF A FAST REACTOR* 
I, I. Bondarenko and Yu. Ya. Stavisskii 


Translated from Atomnaya Energiya, Vol. 7, No. 5, pp. 417-420 


November, 1959 
Original article submitted January 6, 1959 


For many physical experiments carried out with nuclear reactors, the decisive factor is not the average output 
delivered over a long period of time, but rather the peak output delivered over a short period of time as re- 
quired by the experiment. In such cases it appears to be most advantageous to have pulsed cycle operation of 
the nuclear reactor, i.e., a cycle during which the reactor output increases manyfold during a short interval of 
time. The advantage of work in a pulsed cycle is most apparent fora fast reactor since the lifetime of neutrons 
is a minimum in such a reactor in comparison with other types. This paper presents expressions for the duration 
of output pulses in a fast reactor, and expressions describing the time dependence of the output pulses, There 
are also described the possible characteristics of a pulsed fast reactor which could be used for carrying out 


investigations in nuclear physics. 


A reactor operating in a pulsed cycle can be used 
for: 

1) investigations employing time-of-flight 
neutron selectors; 

2) study of short-lived activities (delayed neutrons 
{1, 2], short-lived isomers, etc.); 

3) various investigations for which the main 
difficulty is background not associated with operation 
of the reactor, 

In 1955, D. I. Blokhintsev proposed the use of a 
reactor operating in a pulsed cycle in which the output 
pulses were created during the motion of a portion of 
the active material which is fastened to a rotating 
disk. In this case the reactor stays subcritical for a 
long time. It passes to a supercritical state for short 
periods of time when the moving portion of the active 
material passes through the main fixed portion, 

The change in the reactor output during the output 
pulse can be expressed with adequate accuracy by 
means of the following equation: 


w(t) (t) S 
di (1) = +=) dt. (1) 
Here w(t) is the instantaneous value of the reactor 
output (number of fissions per second); T is the mean 
lifetime of the prompt neutrons?; v is the number of 
secondary neutrons; S is the output of the external 
neutron source; € (t) is the excess multiplication factor 
for prompt neutrons in the reactor. The effect of 
delayed neutrons is not considered by this equation, 
since their lifetime under the conditions of interest is 
much longer than the duration of the output burst. 
In order to obtain a short output burst, it is ne- 


cessary to have €(t) positive for a short time. We 
shall take for our initial time point (t = 0)that mo- 
ment when €(t) has its maximum value (see figure). 
We shall consider the function €(t) to be symmetrical 
with respect to the point t = 0. This condition is 
obviously satisfied if the reactivity change is accomp- 
lished by uniform motion of the movable portion of the 
active material, At a point of time distant from the 
beginning of the pulse, €(t) has a constant negative 
value équal to €;;. The reactor output at this time is 
also constant; 


After this €(t) begins to increase. During the time that 
the rate of increase of €(t) is sufficiently small and 
yet its value is reasonably far from zero, the process 
has a quasistatic character. During this period the 
reactor output stays close to its equilibrium value: 


S 
w()= — (3) 


This situation continues to exist as long as the following 
condition is satisfied: 


dw (t d 
(4) 


* This paper was written using the results of work 


completed in 1956, Certain data used herein were 
also obtained by T. N. Zubarev, 


+ For determination of the mean lifetime of the 
prompt neutrons see the work of L. N. Usachev [3]. 
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Graph of output pulse generation for pulsed 
operation of a fast reactor. 


Finally, the output begins to increase, more slowly 
than equation (3) indicates. Calculations show that 
when €(t) becomes equal to zero, the reactor output 
reaches a value of 


1, 25S 
(5) 


where y = de(t)/dt|, - -tg After this comes the main 
increase of output associated with positive values of 
€(t). In describing the output change during this period 
we can neglect the second term of the right-hand side 
of equation(1). This does not produce significant 
errors during the pulses in question in which the reactor 
output increases manyfold, It is precisely this case 
which is of primary interest, 

Then 


w(— ty) 


dw(t) _e(t)dt 
w(t) 


(6) 


The maximum output value Win’ which occurs at t = to 
(when €(t) again becomes zero), is 


—lo 


to 


(7) 


Let.us examine the pulse shape in the vicinity 
of the maximum, At the times t, close to to, the 
dependence of €(t) on t can be expressed as follows: 


e(t)= —y(t—t). (8) 


Then from equation (6) we obtain 


[ 


Thus, in the neighborhood of the maximum pulse we 
have the shape of a Gaussian curve with a half-width 
6; 


2,35 V= (9) 


At a distance from the maximum the pulse shape 
differs from a Gaussian curve; the pulse rises somewhat 
more slowly and falls more rapidly. However, nume- 
rical calculations per equation (1) show that for large 
output pulses there is significant deviation from a 
Gaussian curve only at the edges of the pulse; this 
makes only a small contribution to the total pulse 
energy. 

The maximum value of €(t) in reactors, equal 
to € |» is always small with respect to unity. There- 
fore the dependence of € on the position of the moving 
portion of the active material can be approximated 
by the parabola 


where @ is a quantity dependent on parameters of the 
moving portion of the active material and of the 
reactor, It characterizes the steepness of the change 
in reactivity with motion of the movable portion of 
the active material. The time-dependence of € is 
expressed by means of the quantity a and of the velocity 
v of the moving portion of the active material, as 
follows; 


& = €,,— avl?, 


Using equations (7) and (10), we obtain an expression 
for the total number of fissions during a pulse: 


E= SK (€m); 
1,25 1,33e9/2 
va. vta‘/2 (11) 


We are also interested in the dependence of the 
pulse duration on the basic quantities which charac- 


| 
| 
| 
3 
Em 
t 
ty 
4,258 e (t) dt 
= ee exp 
Vv 
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terize the reactor (T, ©, v). From the relations given 
above one can obtain the following expression: 


6 = 1,74 A. (12) 


Thus the pulse duration lessens with decreasing T and 
with increasing « and vy. The pulse duration also 
depends on the pulse energy (more precisely, on the 
relation of the pulse energy to the output of external 
neutron sources). However, this dependence is quite 
weak over a wide range of pulse energy. In equation 


(12) this dependence is contained within the coeffi- 
cient A; 


A~ [ In | 
Up till now we have examined single output 

pulses, We shall now examine those cases in which 

the reactor operates periodically, yielding n pulses 

per second, For such an operating reactor cycle 

we shall have an accumulation of delayed neutron sources, 

If the reactor output has any significant magnitude, 

then the neutron sources S which are external with 

respect to the pulse appear to be delayed neutron 


sources which have accumulated during the previous 
pulses; 


S= dedi, 
(13) 
where c; is the concentration of previously delayed 
neutrons of the ith group; dj; is the corresponding decay 
constant. If the interval between pulses is small in 
comparison with the lifetime of the delayed neutrons, 


then the average reactor output W(t) can be expressed 
by the equation 


W()=EW)n+ Ene, (14) 


Here the second term corresponds to the output occurr- 
ing during the time between pulses, In those cases 
‘which are of more interest, practically the entire 


output occurs during the pulses, A necessary condition 
for this is 


eu 


< 1. 


It is precisely these cases which are examined below, 
Taking into account the accumulation and decay 
of delayed neutrons and using equation (11) we obtain 


the following equation describing the change of average 
output; 


W (t) K (8m) 
dc; (t 
(a5 


where 8 ; is the yield of delayed neutrons of the ith 
group; 2 B; =f . From these equations it is 


easy to obtain the conditions for quasistationary 


reactor operation (pulse amplitude remains time- 
constant); 


K (£,,) nb = 1. (16) 

We shall designate by €) that value of € ,, which 
satisfies equation (16), Thus a self-sustaining chain 
reaction ensues when the reactor operates with € pm = 
= € 9; this corresponds to operation of an ordinary 
reaction in a steady state with Kef¢= 1 (Keff is the 
multiplication constant), For € ,, > €9,the reactivity 
increases (pulse amplitude increases); for € , < €9, 
the chain reaction dies out. It is interesting to note 
that for small perturbations (when € rp is close to € 9) 
the time-dependence of the average output of a pulsed 
reactor is analogous to the time-dependence of the 
output in ordinary fast reactors, 

To show this, we shall write the kinetic equations 
for ordinary reactors as follows: 


dei (2) 1 17 
de; 

Fe = — dics (t)-+ W (t) vB;. (17 b) 


Here it is assumed that the lifetime of prompt neutrons 
equals zero, For reactors in which the prompt neutrons 
have short lifetimes (fast reactors) this assumption is 
always valid in terms of delayed criticality. It can 
easily be shown that for small perturbations of the 
system, equation (15) describing operation of a pulsed 
reactor coincides with the system (17) provided we 
designate 
K (&) 
aK (&m) 


d&m 


by B% 


Thus, the time-dependence of the average putput 

of a pulsed reactor turns out to be the same as that 

of an ordinary reactor having a certain effective 
portion of delayed neutrons. Therefore, a pulsed 
reactor can be regulated by the same system as or- 
dinary reactors, Furthermore, changing the regulating 
device would lead to change in € pp). 


For large perturbations the kinetics of the pulsed 
reactor differ from those of ordinary reactors, This 


difference manifests itself if only by the fact that for 
ordinary reactors, as opposed to pulsed reactors, the 
regulating influence of delayed neutrons ceases for 

€ >86; this materially changes the reactor kinetics, 

As an example,we shall use some calculate:| 

results for one model of fast reactor operating ir a 
pulsed cycle. In this model the stationary part of the 
active material is plutonium, while the movable part 
is U™, The values of the usual reactor parameters 
which determine its pulsed operation are as follows; 
a= 7,5 x 1079 r=10° sec; 0.2, For 


velocities of the moving portion of the active mate- 
rial equal to 300 m/sec, and with a pulse frequency 
10 pulses/ sec, the pulse duration is ~ 13 - 10~* sec 

for a constant value of prompt subcriticality € 9 = 2.6 ° 
107°, The effective fraction of delayed neutrons 


corresponding to this cycle is 2.4 - 10, for an average 
reactor power output of 10 kwt, the power during 
pulses reaches ~ 10° kwt. 

In conclusion the authors wish to express their 
appreciation to D, I. Blokhintsev, A. I. Leipunskoi, 
and O. D, Kazachkovskoi for constant attention and 
advice, and also to F, I, Ukraintsevy and N. V. Kras- 
noyarov for valuable discussions, 
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HEAT-GENERATION STABILITY IN A PULSED 


Translated from Atomnaya Fnergiya, Vol. 7, No. 5, pp. 421-423 


November, 1959 
Original article submitted December 19, 1958 


neutron pulse. 


In [1] the operation of a pulsed reactor is in- 
vestigated, where in order to obtain a stable heat 
output, the duration of the neutron pulses is automa- 
tically decreased or increased for increasing or dec- 
reasing reactivity in the reactor during a neutron pulse. 
We will consider the dependence of the heat output 
from a pulsed reactor on various factors characterizing 
the intensity of the neutron pulses for the conditions, 


TK () 
Ak 
ait (2) 
| Ak (t)| < 4. (3) 


Here T is the neutron lifetime in the reactor, 
Tg is the mean half-life of the nuclear fragment 
emitting the delayed neutrons, w = 2m/Tp (Tp is the 
period of one neutron pulse), 4 k(t) = ko¢¢(t) -1 and 
Ako = Kote max” 1s where kege(t) and Kore max ate 
the effective neutron multiplications in the reactor 
at time t and when the maximum reactivity is attained, 
respectively,with the delayed neutrons neglected. 

If conditions (1) and (3) are satisfied and Ak(t) < 
< 0, then the function N(t) giving the number of 


neutrons in the pulsed reactor in terms of the time 
satisfies the relation 


N N 
dN = 5 Ak (t) dt + (4) 


where N ais the equilibrium number of delayed 
neutrons in the reactor determined from the mean 
reactor power. If the heat generated during the time 
between neutron pulses is significantly less than the 
heat generated during the pulses, then 


QvpT 


(5) 


3 


The conditions for heat-generation stability in a pulsed reactor are investigated, and formulas are obtained for 
the calculation of the heat output for the cases of linear and sinusoidal laws of variation reactivity during the 


where Q and qp are the quantities of heat generated 
in the reactor during one neutron pulse and from a 
single uranium nucleus respectively, 8 is the pro- 
portion of delayed neutrons in the total number v of 
neutrons, and T,, is the time between two consecutive 
neutron pulses, 

The relation (5) is valid in the case of a constant 
heat output during the neutron pulses when no nuclear 
fragments are lost from the active zone. If the heat 
output varies gradually from pulse to pulse, then in 
(5), instead of Q, we must substitute 


Tp Qi 
Tq (6) 


where Q; is the heat generated during the ith pulse, 
and the summation extends over the pulses occurring 
during the time T,. 

If Ak(t) > 0 the sources of delayed neutrons can 
be neglected, and (4) takes the form 


dN Ak (t)de. (7) 


The quantity of heat generated in the reactor 
during a single neutron pulse will be 


Q=7\ Nod, 


where q is the amount of heat generated during the 
absorption of one neutron, 

Let us assume that the reactivity varies during a 
neutron pulse according to the law 


Ak(t)=Aky= at for (9 


Ak(t)=Ak,(1—Zot) for 
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where for the initial instant of time (t = t,) 


Ak (t;)= —Ak,, > 1. (11) 
=0 
If we use (9) and integrate (4) from tain ® 4 


we obtain the value of N(t) for Ak(t) = 0; 


Na 
Ako x 


If the conditions (1), (2), and (11) are satisfied, 
then (12) yields 


We will denote the quantity of heat generated 
when Ak(t) varies from 0 to Ako, from Akg to 0, and 
from 0 to (-Ak) by Qy, and Qs sespectively. Then 
if Ak(t) varies according to (9) and (10), the relations 
(7), (8), and (13) yield 


3/ 2 


0 
If we use the relation 


2 
0 


ex 


where M(-1/., x”) is the confluent hypergeometric 
function as in [2], then for condition (2) the relations 
(14) and (16) give the total heat output during a 
neutron pulse for the variation of Ak(t) according to 
the linear law 


2 Na 


Aho 
gor, (17) 


The greater part of the heat is generated while 
Ak(t) is decreasing << while approximately 
equal amounts of heat are generated during the varia- 
tion of Ak(t) from Ak, to zero and during the damp- 
ing of the chain reaction (Q, ~ Qs). The maximum 
neutron flux in the reactor is reached while Ak(t) is 
decreasing, at the instant of time when Ak(t) = 0. 

In contradistinction to the case we have considered 
when Ak(t) varies linearly, some pulsed reactor designs 
have Ak(t) varying relatively smoothly in the region 
of its maximum value. The corresponding correction 
to the heat output necessary in this case is significant, 
and can be easily found, 

Let Ak(t) vary according to the law 


Ak (t) = Ak, cos ot. (18) 


Then if (7) is used, the relations (8) and (13) and 
the condition (2) yield 


Ako 
No bk 


where Iy (Ako/wT) is the Bessel function of zeroth 
order for an imaginary argument. 

We will use the relations obtained,and investigate 
the conditions for heat output stability in a pulsed 
reactor. We will denote by K the ratio of the quantity 
of heat generated during a given neutron pulse to the 
mean quantity of heat generated during the pulses 
occurring in the interval of time T,. If (17), (19), 
(5), and (6) are used, we obtain 


(19) 


T 
K {(Aky, oT), (20) 


where 9 is the probability that a neutron be damped 
by a uranium nucleus, and f (Ako, wT) is a certain 
function of the neutron-pulse parameters, The form 
of the function f ( Aky, wT) is determined by the law 
of variation of the reactivity during a neutron pulse. 
If the reactivity varies according to (9) and (10), then 


Ako 


3 ol 


and if it varies according to (18) 


Ako 
3/2 @ oF 


‘Akg OT * (22) 


hin (Aky, ~2( ) 


The condition K = 1 can be called the heat-output 


Sky \2 0 
Ak, (S22) 
x {14248 ae \ dz » (12) 
where 
4 le dz, (14) 
0 
2 
(x)*/? ( | 
e - 
Q, ( dz, (15) 
0 
Sky 
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stability condition, since when it is satisfied, the heat 
output during the neutron pulses remains constant 
from pulse to pulse (if there is no external source of 
neutrons), 

In all cases that are interesting in practice Akg > 
> 0,01, 9&1, v 32, 8 ~ 0,008, In such cases, (20)- 
(22) imply that for K = 1 the condition (2) is satisfied 
if 


Tp 
T. > $6, 


From (17) and (19) we can also obtain the relation 
between the relative variation of the heat output 
during a neutron pulse and the variation of other 
parameters characterizing the pulse: 


Ako ( A(Ako) A (oT) 
oT Ako oT /'(23) 


where the coefficient « takes 7/2 and 2for the appli- 
cation of (17) and (19), respectively. 


From a consideration of (23), we see that stability 
of the heat output during the pulse ( AQ/Q = 0) can 
be maintained for variations in, Aky and wT that are 
not too large, if the condition 


A(Ako) _ A(wT) 
(24) 


is satisfied. 
The author thanks Yu. N. Zankova and A, K, Sok- 
olova for discussions on the results of the work. 
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The purity of the graphite used in graphite-moderated thermal reactors is a factor of prime importance. This 
article describes a method for testing the reactor feasibility of graphite. The gist of the method involves com- 
parison of flux attenuation of thermal neutrons in the test graphite and in a standard graphite specimen. The 
method serves in determining the diffusion length L or the effective absorption cross sectiono of graphite, where 
Lo and 09 for the standard graphite are known, The error incurred in determining o does not exceed 2.5%, The 
simplicity of the method recommends its use for large-scale testing of graphite under factory conditions, 


Methods for testing the reactor feasibility of gra- 
phite boil down in the last analysis to a determination 
of the thermal neutron absorption cross section o. In 
the simplest steady-state diffusion method, the value 
of o is found, as we know, by measuring the diffusion 
length L of thermal neutrons in a graphite prism (the 
prism method). This requires a neutron source and 
neutron detector. Suitable detectors for the purpose 
are, e.g., indium foils with facilities for measuring 
B- and y -activities, or BFs proportional counters with 
associated instrumentation. The amount of graphite 
required in the prism to provide satisfactory precision 
in the measurement of L should not be less than 6-7 
tons, 

Much less graphite is required in the danger- 
coefficient, pile-oscillator, or nonsteady-state diffu- 
sion methods [1-3], but the equipment and instrumen- 
tation is much more involved and expensive than 
that used in the prism method, and more highly skilled 
personnel are required to service and operate them. 

The prism method therefore appears to be the 
most preferable for large-scale testing of graphite under 
factory conditions. However, one drawback in this 
method is the need to use a relatively large amount of 
graphite. This places limitations on the possibility of 
grading the graphite asto quality, and renders quality 
control of small batches of the graphite more difficult. 
Below, the comparative method of graphite testing, 
which is free from this shortcoming, is discussed. This 
method was devised some time ago in the USSR, and 
has been applied successfully under mass-production 
conditions after slight modifications, 


DESCRIPTION OF THE COMPARATIVE METHOD 


Theory, The essence of the method was first 
discussed theoretically by V. S, Fursov, and found as a 
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result of a comparison of the test graphite with refe- 
rence graphite of known capture cross section 09. 

Let a source be placed at a point O on the axis of 
the prism, the latter being a long rectangular prism 
consisting of the standard graphite, and let a neutron 
detector be placed at point D, with the distance sepa- 
rating neutron source and detector such that the mea- 
sured density of thermal neutrons ng at.point D may be 
represented by only the first harmonic (Fig. 1). By 
replacing a portion h of the prism by the graphite 
specimen whose absorption cross section must be mea- 
sured, and measuring the neutron density n at the 
same point D, the condition of equality between the 
neutron densities and fluxes at the boundaries of the 
insert of the thickness h allows us to derive the follow- 
ing expression; 


Rd_ 
Rodo 


Rd \2 Rd \ 
) ; 


4 


(1) 
Here, Ro and R are the relaxation lengths of ther- 
mal neutrons in the prism, for the reference and the 
test graphite respectively, determined from the 
familiar equation 
1 


=> 


a 


where L = VIZXX,, is the diffusion length; agp = a + 
+2X 0.71A,, and begs = b + 2X 0.71 Agy are the effective 
transverse sections of the prism; A is the absorption 
length; A+, is the transport length; a and b are the geo- 
metric transverse sections of the prism; dy and d are the 
densities of the reference and test graphite specimens, 


8 n 
| 
eff ‘eff 
(2) 


Under conditions of practical interest, the values 
of 6 and do/d could not differ from unity by more 
than + 10%, Therefore, when measuring 6 to an 
accuracy not exceeding 0,5-1%, (1) may be simplified 
by reduction to the form 


(3) 


Equations (2) and (3) yield the following simple 
relationship for 9 and oo: 


1 In af 4 ) 
h batt (2) 
1 1 1 \d 
Re ™ 
eff eff 


The sensitivity of the method is determined by 
the expression 


(4) 


2 (3 2 
R, dad 5 
A = 1 j ( ) 
eff eff 


For the prism dimensions indicated in Fig. 2, a devia- 
tion of the value of o by 3-4%(0.1 millibarn) from 
Oy of the standard leads to a 1% change in 6. 

Practical execution of the method. Several 
requirements must be fulfilled with respect to the 
reference graphite, the geometry of the measurements, 
the neutron source, and the neutron detectors, for the 
comparative graphite testing method to prove conve- 
nient under mass-production conditions, and for it to 
yield satisfactory results when using the modified for- 
mula (4). Let us discuss those requirements briefly. 

The density and capture cross section of the re- 
ference graphite must be roughly equal to the mean 
values of the density and capture cross section of the 
graphite. When this condition is observed, the devia- 
tions of 6 and d)/d from unity will be least, and the 
results obtained from the approximate formula (4) 
and the exact formula (1) will be negligible. 

The effective absorption cross section of the 
reference graphite must be measured exactly, e.g., 


Fig. 1. Measurement diagram in the compa- 
rative method; O) source; D) neutron detec- 
tor; a) geometric transverse cross section of 
prism. The hatched area of the prism of re- 
ference graphite, of thickness h, is replaced 
by the test graphite. 


by the prism method, Due attention must be paid 
to any possible inhomogeneity of the graphite in 
the process, 

After measuring G9, the reference graphite spe- 
cimen should be placed, for ease of handling, in the 
form of a horizontal rectangular prism on a low, firm 
table, and covered over on all sides by a cadmium 
sheet 0.5-1 mm thick, Operating experience has 
demonstrated that the geometry of the measurements 
as shown in Fig. 2 is entirely suitable, when using 
a Ra-a-Be source. About 10-12 tons of reference 
graphite are required to set up a prism of that type. 
The test graphite, in a block amounting to some ~ 1,5 
tons, is inserted into the hatched portion of the prism. 
The thickness of the test graphite insert h must not 
be less than L, Any increase in weight of the insert 
over and above 1.5 tons will render quality sorting 
of the graphite more difficult, and the method will 
lose its feasibility. 

Neutron detectors are placed in the plane passing 
through the center of the prism at right angles to 
the prism axis. As detectors, we used one to three 
BF3 proportional counters (points D, Dy and D,), con- 
nected in parallel, with DA type instruments, The 
neutron source was placed at points Cy and C, alter- 
nately. The ratio of counting rates I, and I, with 
the source positioned at points C, and C,, respectively, 
obviously yields the ratio of neutron densities 6 = 
= n/ ng = 1,/ 1, appearing in the formulas, 

The distance between the source. and the insert 
of test graphite must be such as to provide for prac- 


Fig. 2. Reference prism for testing graphite 
by the comparative method: a) outer view; 
b) cross section through ABC, 
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tically complete slowing down of the neutrons to 
thermal energies, The distance from the edge of the 
insert to the detector must be of the order of L, so 
that the test graphite might not affect the counting 
rate I, when the source is located at point C,, How- 
ever, it was demonstrated by experimental testing 
that this distance may be reduced to 30 cm. This 
enhances the sensitivity of the method to contamina- 
tions in the test graphite, and makes it possible to 
employ a source of lower power. 

The distance from the source to the face of the 
prism and the transverse cross sections of the prism 
are chosen so that the neutrons from the source 
will have time to slow down below the cadmium ab- 
sorption cutoff before leaving the prism, When this 
condition is not observed, the background of scattered 
neutrons will adversely affect the conditions under 
which measurements are performed, 

The basic requirement concerning the neutron 
detector is constant amplification factor on the part 
of the electronic equipment used for recording neutron 
counts, although the alternating measurement of 
values of I, and I, at points Cy and C2of theprism in the 
sequence I, I, > I, > Ij... may eliminate any effect 
of monotonic drift of the amplification factor on the 
6 ratio, The requirement of high resolving power 
(short dead time) is less important, since corrections 
for spurious counts cancel out at a 6 = I,/ I, ratio 
close to unity. 

Accuracy and Capacity of the Method. With a 
neutron source and detectors as described above, 
and reference prism dimensions as indicated in 
Fig. 2, 6 may be measured with ease to an accuracy 
of 0.2-0.3%, i.e., 6 may be determined to an accu- 
racy of 0.6-1,0%(without considering the systematic 
error introduced by inaccuracies in the determination 
of Oy of the reference graphite), 

However, corrections taking into account inho- 
mogeneity in the density of the graphite, the insig- 
nificant contribution of the higher harmonics of the 
neutron flux at the left-hand edge of the insert, and 
the use of the modified equation (4) limit the accu- 
racy attainable in determining o to a value of 2-2.5% 
so that 5 may be measured to an accuracy of 0.5%, 

The effectiveness of the method depends on the 
power of the source and on the number of neutron 
detectors and their efficiencies. For example, when 
a Ra- &-Be source of 300-500 mC activity and a 
proportional counter 300-400 cm’ in volume, loaded 
with BFs to a pressure of 500-600 mm Hg, are used, 
20-30 min are required for a measurement of o for 
one insert (~ 1.5 ton graphite), with an error not 
in excess of t 2.5%, 


By using less efficient detectors or increasing 
the number used, it is possible to reduce testing time 
for a single graphite insert, and to reduce the needed 
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source power. In particular, the use of two or three, 
or a larger number of detectors placed on the plane 
passing through the center of the prism and perpendi- 
cular to the axis of the prism, has the further advan- 
tage that sensitivity to contaminations of the graphite 
in the peripheral regions of the insert is enhanced, 
producing a slight effect on the value of 6 when a 
single detector is employed at point D, Several 
detectors may be used to advantage, particularly in 
the case where the test graphite is inhomogeneous 
with respect to density or purity, By placing detectors 
of different sensitivities at pointsD, Dy, and Dz it is 
possible to achieve the effect of having all contami- 
nations in any part of the insert exert pretty much 
the same effect on the value of the 6 ratio, The 
possible inhomogeneity of the graphite must also be 
kept in mind when comparing the results of measure- 
ment of o by the comparative method described 
herein, and by the exponential prism method. We 
carried out such a test for several batches of graphite, 
and the results of measurements of o proved to be the 
same within the limits of experimental error. 


ADVANTAGES AND DISADVANTAGES OF THE 
METHOD 


The comparative graphite testing procedure has 
the following advantages; 1) simplicity on the part 
of the formula used in calculations, and also on the 
part of the instrumentation, thus not requiring highly 
skilled operators; 2) high capacity, satisfying the 
demands of mass production; 3) moderate weight (1.5 
tons) of the test graphite, which makes it possible to 
exercise quality monitoring of the graphite, reducing 
scrap to a minimum, The latter factor is especially 
important for industry’s needs, in mastering the tech- 
nology of the production of low-ash graphite for 
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Fig. 3. Analysis of results of graphite tests by 
the comparative method, 
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nuclear reactors, At the same time, the graphite 
obtained differed appreciably in quality from batch 

to batch, depending not only on the raw material used, 
but also on what part of the graphitizing furnace the 
raw material was placed in during the graphitization 
process, Grading of the graphite as toquality provided 
the possibility of deploying the graphite more ration- 
ally in nuclear reactors. 

Fig. 3 shows the results of tests of several hundreds 
of tons of graphite using the comparative method 
for the use of the first nuclear reactor in the USSR [4]. 
The entire lot of tested graphite was divided into five 
groups in increasing values of 0, obtained from (4), 
and the mean value 0 was determined, along with the 
ash content q (ratio of weight of residue after ashing 
of the graphite specimen to the weight of the specimen 
in percentage figures) for each group, Points plotted 
on the graph correspond to those mean values at which 
the absorption cross section 6 is approximately directly 
proportional to the ash content q. The range over 
which the averaging with respect to o was carried out 
is indicated by the vertical lines drawn at the points, 
The absolute values of o were reduced by approxi- 
mately 15-20% to account for inaccuracies in the 
values of 69 and ),, available at the time for the 
reference graphite. As the technology was improved, 
the quality of the graphite was improved, the ash 
content and the value of the absorption cross section 
were reduced, The diffusion length was measured in 
. @ graphite cube having sides 6 meters in length along 
the edges, The cube was composed of medium-qua- 
lity graphite prepared for the first reactor. The 
value of L was 50 4 1 cm (brought into correspondence 
with the density 1.65 g/cm’). 

Note that the value of o to be calculated on the 
basis of (4) is sensitive not only to impurities in the 
graphite, but also to the microstructure of the latter, 
since the scattering cross section o, and the transport 


cross section 0,, may both be dependent on the 
dimensions and orientation of the crystalline grains. 
However, measurements of 0, for several dozen 
batches of graphite yielded results which showed 
agreement within the limits of experimental error, 
which bears evidence of the absence of any difference 
in the microstructure of those batches of graphite. 
When the need arises to test graphite featuring 
some appreciable difference in microstructure, this 
disadvantage may be obviated by replacing (4) by the 
corresponding relationships for diffusion lengths and 
by the additional measurement of transport length. 
Inconveniences in the use of the method under 
factory conditions consist in the need to observe 
approximately identical temperature conditions in 
both test and reference graphite specimens, and to 
maintain constant atmospheric conditions (par- 
ticularly humidity) in the rooms where the graphite 
is being tested. Humidification of the graphite by 
absorption of water vapor from the surrounding air 
may result in a pronounced increase in the 09 of the 
reference graphite, Periodic monitoring of the value 
of Oo of the reference graphite is therefore mandatory. 
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The design and building of reactors requires a knowledge of the interaction cross sections of neutrons with diffe- 
rent materials over a broad energy range. Interactions between neutrons and nuclei (radiative capture, fission, 
scattering) are viewed in the light of current concepts as processes involving the formation of a compound nucleus 
and displaying a clearly pronounced resonance response; the interaction increases sharply at certain neutron ener- 
gies. The position and quantitative value of the resonance levels still defy theoretical prediction, and direct 
measurements of the cross sections concerned are consequently the sole sources of the required information. 

This article comprises a survey of papers appearing between the first and second Geneva conferences (1955 
and 1958) on the peaceful uses of atomic energy, dealing with the study of the properties of resonances of fission- 
able nuclei. Various methods of describing the interactions between neutrons and fissionable nuclei are consi- 
dered in the light of the available data, The experimental data have yet to be interpreted fully and unambigu- 
ously, though nonetheless some characteristic aspects of the fission process are already apparent. A brief descrip- 
tion is given of several refinements in technique employed in the measurement of resonance level parameters, 


and the results thereby obtained are discussed. 


Introduction 

The study of interaction between resonance neut- 
rons and fissionable nuclei is being conducted inten- 
sively by many teams of scientists in different coun- 
tries. The interest in these investigations is due,to an 
appreciable extent,to the creation of intermediate 
neutron and fast neutron reactors, where a knowledge 
of the interaction cross sections of neutrons with ener- 
gies reaching tenths of a kiloelectron volt proved 
crucial in their design. A considerable portion of the 
neutron spectrum falls within the resonance region in 
such reactors, and the properties of the resonances 
play a large part in determining such reactor para- 
meters as critical reactor dimensions, reactivity 
margin, temperature coefficient of reactivity, etc. 

A knowledge of cross sections in the resonance region 
is also necessary for designing thermal reactors. 

In addition, the application of the collective 
model to fissionable nuclei [1, 2] also stimulated a 
whole series of experiments related to verifying the 
existence of such effects as interference between 
levels, mass yield of fission fragments as a function 
of the spin of the compound nucleus, and still others. 

The significant refinements introduced into the 
measurements procedure in the period spanning the 
two Geneva conferences have made it possible to 
investigate in detail the behavior of fission and 
capture cross sections at low neutron energies to deter- 
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mine the parameters of a large number of levels, 


and to obtain an answer to some theoretical questions, 

An attempt is made below to systematize the 
data obtained in recent years on interactions between 
neutrons and fissionable nuclei in the resonance region, 
Several problems considered in detail in previous 
survey papers [3, 4] (effect of asphericity of nucleus, 
repulsion between levels, etc.) will not be touched 
upon in this article. 
Applicability of the Breit-Wigner For- 
mula: The Generalized Formula 

The best description presently available for 
interactions between neutrons and nuclei is the model 
of the semitransparent nucleus [5]. It is assumed in 
this approach that the interaction may be divided 
into two stages [6]: the neutron first interacts with 
the nucleus as a whole (this interaction is described 
by the complex potential taking neutron capture into 
account), and then with discrete nucleons or groups 
of nucleons within the nucleus, It is only in some 
cases that a compound nucleus, i,e., a nucleus in 
which the energy of the captured neutron is distri- 
buted over the entire nuclear volume, is formed at 
the termination of the process, In other cases, collec- 
tive and single-particle energy levels of the system 
(to use the terminology in [6]) are excited. 

Questions related to the first stage of interaction 
with a nucleus (potential scattering, transparency of 


the nuclear surface, etc.), have been studied in detail 
in the literature [3, 4]. The present article will con- 
sider data relating to reactions involving the formation 
of a compound-nucleus system. 

The theory of nuclear resonance reactions has 
been developed in a series of papers [5-10], and is based 
on a minimum number of assumptions and is charac- 
terized by great generality; however its practical 
use requires different degrees of approximation because 
of the large number of undetermined parameters, 

The cross sections for radiative capture of slow 
neutrons by nonfissioning nuclei are well described by 
the Breit-Wigner relation [11] for an isolated level, 
as shown with great accuracy for silver [12] and gold 
[13]. After the discovery was made that fission cross 
sections show a clearly defined resonant character, 
these were also described by the analogy of the 
Breit- Wigner formula, e.g., in the following form: 

Ty 
TT oy) 

4(E—E,)?+T? 


Here o, is the absorption cross section; o,¢ is the 


Onf Onyx; Onf = 4nk? 


Ony = 40k? 


fission cross section; ony 
cross section; % is the wavelength of the incident 
neutron divided by 2m; g is a statistical factor equal 


is the radiative capture 


to 2 2 7-+-1,/ for slow neutrons; I is the spin 
of the target nucleus; E is the neutron energy; Ep is 
the energy of the resonance level;I),, are the 


neutron, partial, and total level widths, respectively. 

Various methods exist for determining the level 
parameters, based on the use of this formula [4-16], 
and are widely used for processing experimental data; 
a large portion of the level parameters of fissionable 
isotopes known at the present time was obtained through 
just this approach, 

However, as the accuracy of the measurements 
was improved, it was found that the cross sections of 
fissionable isotopes could not be successfully described 
by the use of (1), because several resonances displayed 
an explicitly asymmetric shape. An indication of this 
difficulty is reported in a paper by Sailor [17], who 
studied the cross sections of U™*, This author reached 
the conclusion that interference effects between 
levels had to be looked into. Other authors [18, 19] 
came to a similar conclusion, Figure 1 shows the trans- 
mission of neutrons by a Pu” sample in the resonance 
region embracing energies of 10.95 to 11,90 ev [19]. 
Here, we see the results of calculating the transmission 
of a given specimen by the Breit-Wigner relation; the 
parameters of the levels are taken from [19]. The 
measured cross section value between resonances 
differs considerably from the theoretically predicted 
value, The deviation of the shape of the resonance 


peaks from the Breit-Wigner shape may also be ex- 
plained by the presence of forbidden levels. Such an 
assumption appears to be a bit exaggerated although 
such levels are actually found to exist in some cases, 
Interference effects may play a crucial part only 
in a small number of reaction channels; it is demons- 
trated in one paper [1] that fission of a nucleus by 
slow neutrons must proceed in just that fashion, But 
the fact that the fission reaction has a small number of 
channels available to it also flows from a consideration 
of the distributions of resonance level widths, Porter 
and Thomas [20] derived equations relating the distri- 
butions of partial level widths to the number of degrees 
of freedom of a particular reaction (Fig. 2), On the 
basis of these relations, it is possible to justify the con- 
clusion that fission of the nucleus probably proceeds 
through,at the most,one or two exit channels, The 
presence of interference necessitates the use of a 
formula for many resonance levels in describing cross 
sections of fissioning nuclei. The most general ex- 
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Fig. 1. Neutron transmission by a sample of Pu” 


in the 10-13 ev energy region, The dashed curve 
was plotted from the Breit- Wigner formula. 


Fig. 2, Distribution of resonance level widths as 


a function of number of degrees of freedom of reac- 
tion 
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pression for the reaction cross section may be obtained 
by means of the scattering matrix, The scattering 
matrix Sy u is indicated in the following notation [21]: 


= exp[— i(k, + hy) RI] Shy, 


= 7 
1—iRjp 


where the matrix is related 


to the derivative matrix Ry by the equation R} hes 
(kk) Ryu The derivative matrix Ry may be 
indicated in turn by the notation; 


Yen, Yeu 
Ryu Eg E 

8 
The following conventions are used in these formulas; 
R is the nuclear radius (in the channel); A, # are sub- 
scripts indicating the reaction channels; ky, ky are the 
wave numbers of the channels for neutron energy E; s 
indicates the resonance level; y,, and ys, are quan- 
tities associated with the partial level widths by the 
relation =2 ky 


A formula derived in a previous paper [5] takes 
into account interference: between resonances, for 


the case of widely separated levels (r/ D << 1, where 
D is the level spacing), Resolving matrix Ru into 


two parts 


Bey B Bia B 


\2/3 
where B,, = (ky) (+ ) the authors 
then expanded the matrix sy in powers of 


as the order of the magnitude of I/D) and, discarding 
all but the first three terms of the expansion, derived 
the following expression for the reaction cross section 


Wid 


= 
8 


(assuming Ty = 0, we obtain the familiar Breit- 
Wigner formula), But a formula of this form cannot 
accurately describe the energy dependence of the 
cross sections in those regions where level widths are 
of the same order as the spacing between them. Since 
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Fig. 3. Total cross section and fission cross section of U™*, Solid-line curves were plotted from the 


generalized formula, 
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the inequality T/D << 1 is fulfilled only very appro- 
ximately for fissionable isotopes and overlapping of 
levels is not infrequent, the range of applicability 

of (2) is rather restricted. : 


Reich and Moore [22] have derived a generalized 
formula free from assumptions requiring that I'/ D 
be small, with interference between many resonance 
levels accounted for. They proceeded according to the 


following reasoning: a) there is only one channel 
available for scattering, since only neutrons with 

orbital moment J = 0 interact strongly with the 
~ nucleus in the resonance energy region; b) radiative 


capture takes place through a large number of channels; 


c) the fission reaction takes place through only one 
channel; the levels interfere fully with one another. 
Below, we have the expression for the fission cross 
section presented by Reich and Moore [22]: 
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An expression was also obtained for the radiative 
capture cross section, formally differing from the sum 
of terms computed on the basis of the Breit- Wigner 
relation, but calculations performed on the basis of 
both formulas gave practically coincident results, This 
stands to reason, since radiative capture takes place 
via a large number of channels n (30< n < 500 for 
o”, as evaluated in [20]), and interference terms 
cancel out in the cross section. 

Figure 3 shows the total and fission cross sections 
of U™* measured over the energy interval from 1 to 
11 ev, and cross sections plotted from (3) by suitable 
choice of parameters [23]. It is evident from the 
diagram that empirically plotted points show satis- 
factory agreement with the theoretically predicted 
curve, The results of similar calculations for U™ 
(energy of primary neutrons ranging 0.1 to 10 ev) are 
taken from [24], and agree qualitatively with empi- 
rically obtained values (Fig. 4). The generalized 
formula provides more accurate agreement between 
the experimental and predicted curves, without any 
assumption as to the presence of additional forbidden 
levels, but its use is restricted to the region of good 
spectrometer resolution and involves cumbersome 
computational labor, As previously indicated, a 
large portion of the data on level parameters was 
obtained by methods based on the validity of the Breit- 
Wigner relation for an isolated level. Comparison of 
the parameters of strong resonance determined by 
choosing the parameters in conformity with the 
generalized formula [23, 24] with the results of earlier 
papers [25, 26] demonstrates that the discrepancies 
in the values reported do not exceed the limits of 
experimental error (Table 1). 


Interaction Cross Sections 

Averaged results of measurements of the inter- 
action cross sections of resonance levels with fission- 
ing nuclei are reported in the literature [27]. Figure 5 
shows, as an example, the neutron energy dependence 


(3) 


of the total cross section of Pu™ in the region from 1 
ev to 10 kev. 

The resolution of neutron selectors now available 
enables us to reliably separate the resonance levels 
of plutonium to energies of the order of 30-40 ev. 

At high energies, only separate groups of resonances 
may be successfully separated. The small distances 
between levels in U™* and u™ compared to Pu” 
(0.9, 0.63, and 2.9 ev, respectively,[28, 29]) make it 
still more difficult to effect a separation of the re- 
sonance levels of these nuclides, 

The many-level formula (3) may be enlisted to 
successfully describe, with adequate accuracy, the 
appearance of cross sections down to energies of the 
order of 10 ev, but an additional nonresonance term 
must be introduced into the cross sections (~ 70/7E 
for U™ and 20/7E for Pu and U™ [30}) to obtain 
agreement between the predicted and experimental 
curves at low neutron energies. The necessity for 
introducing such a term is related to the anomalously 
large thermal cross sections of all three ot these 
nuclides, The presence of a nonresonance term in 
the cross sections may be explained, for instance, by 
the existence either of a “negative” energy level, 
or the existence of weak overlapping levels, 

There is no basis for rejecting the possibility of 
the existence of “negative” levels in the neighborhood 
of zero energy, but values far exceeding the average 
values taken over the resonance region are obtained 
when the parameters of these levels are found, For 
example, for the resonance level of U™ at energy 
Eo = -1.4 ev, we get the quantity 2g 4 = 3.056 Mev 

[24] at Tt 0,18 Mev, while analysis of the cross 
section of U™* yields Ey = -5.115 ev and 2g T = 6.48 
Mev [27] at 0.10 Mev (Ij) = being the 
reduced neutron width), 


The assumption as to the existence of broad over- 
lapping levels is related to the hypothesis of a strong 
dependence of the level parameters on the spin of the 
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TABLE 1. Comparison of Parameters of Strong Resonance Levels 


290,9-106, EV 


Ny. 103,ev 


ry. 103,eV 


formula 


Breit-Wigner | Generalized 


formula 


Breit-Wigner |Generalized 
formula formula 


Breit- Wigner 
formula 


Generalized 
formula 


U283 [23,26] 


U235 [24] 
27,9 
41,9 
37 
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71,8 
92,1 
45 
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Fig. 4. Absorption (1), fission (2), and radiative capture (3) cross sections for um, 
curves were plotted from the generalized formula. 
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compound nucleus [1}. A comparison of the curves 
of fission fragment yield for U™ fissioned by thermal 
neutrons and neutrons of energies 1,8 and 2,3 ev, 
corresponding to resonances (with measurements car- 
ried out by radiochemical techniques) is performed 
in [23]. The ratio of the number of fission fragments 
at the peak of the curve to the number of fragments 
at the trough agrees, within the limits of error of the 
measurements, for resonances and shows a slight 
difference for the thermal spectrum (Table 2), but is 
not at variance with the assumption assigning the 
resonance and nonresonance terms to various spin 
states, The inadequate accuracy achieved in these 
measurements must be noted. On the other hand, de- 
termination of the spins of the Pu compound nucleus 
for resonance levels of energies Ey = 41.5 and 52.6 
ev demonstrated that these spins are different [18], 
although the partial widths of these levels do not 
differ greatly. 
The Neutron Width 

The neutron width Ij, more accurately the 
combination gI’,, was measured for fissionable nuclei 
in the same manner as for nonfissionable nuclei, viz., 
by the method of comparing the areas of transmission 
of neutrons by specimens of different thicknesses, 
Because of the nonresonance term in the cross section, 
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TABLE 2. Comparison of Yields of Fission Fragments 
of U™* in Fission by Resonance and Thermal Neutrons 


Yield ratio 
Neutron energy, ev 
Mo99/Agill Mo99/Cd115 
Thermal spectrum 4 1 
1,8 1,14+0,02 
2,0 41,05+0,04 1,09-L0, 11 


interferences between levels and difficulties encoun- 
tered in resolving different levels where instrumental 
resolution fell short of the task, divergences between 
the results reported by different authors were far in 
excess of the standard errors indicated in individual 
papers. However, despite the large spread in reported 
values, it is apparent that the neutron widths of all 
three fissionable nuclides vary over very broad inter- 
vals, 

The distribution of reduced level widths makes 
it possible to reach some conclusions on the number 
of channels in the reaction [20]. Figures6a and 6b 
represent the distributions of reduced neutron widths 
for U™ and Pu as reported in two papers [28, 31], 
and the distributions calculated for a different number 
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Fig. 5. Total cross section of Pu, 
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TABLE 3, Radiation Widths of Resonance Levels, Mev 


233 U235 Pu239 
Source of Source of data | Source of data 
ata | 
by Ep, €V 
0,20|40+4| — | 0,282] 28,8+5 |31+10| 3946 |39+6 | 27,9 [3243 [3545 | 0,297/30-2 |40+5 [39-2 
1,75 | 44+5 |48+13] 2,036] 30,8+8 |27+10| 315 |31+5 | 34,6 |31+5 |33+9 |10,93 
2°32 3848 3,17) — — | — | — — [57+23)11,90 [42+-12/40,9+5 
3.611 54+5|434+9 13,60; — — | — | 37 [3744] — |17,6 [3047] — |39,1+5 
4,70] — | — | 4,847) — — |25+9/25+9 | 25,5 |25+9 | — |22,2 |35+6 
5.75|50+5| — | 6,40} — — — | — |55413/41,4 — |46,8+9 
6,79(50+5| — | 7,10 — -- | — — | — |404+16/44,5 |40+10) — |40,4+10 
9,05/56+6! — | 8,795] — — [38419] — | — |78415)52,6 | — |34,1+5 
10,33}60+6| — |11,65) — -— | — |35+16/67+21 
12,40; — — | — — |39+16) — 
19,4 — | — — 


of degrees of freedom in the reaction, It is clear 
from inspection of Figs.6a and 6b that the experimen- 
tal curves are in good harmony with the distribution 
for v =1, 

In addition to the measurement of the para- 
meters of individual levels, the value of 44 D, 
according to the latest reported measurements com- 
prising (1.04 0.3) + 10° for [23], (1.0 + 0.2) 104 
for [32}, (1.0 + 10 for Pu™ [29] was 
determined, This value was determined on the basis 
of the slope of the curve of the transmission of a thick 
sample as a function of energy. This method is free 
from errors associated with the possibility of trans- 
mission of weak levels and insufficient resolution, 
which appear when averaging neutron widths for 
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discrete levels with accompanying calculations of the 
total number of levels, 

It should be noted here that the value of m/ D 
for U™ found from the first ten levels is almost an 
order of magnitude smaller than that obtained for the 
entire resonance region [31]. It is indicated in [33] 
that the statistical probability of such a group of 
levels appearing is less than 1%, But the anomalously 
low value of [,,/D and the presence of resonances 
having Z # 0 have not yet been successfully interpreted 
[31]. 


The Radiation Width 

The radiation width may be obtained from the 
radiative capture cross section. In the region of 
neutron energies where resolution is rather high the 
capture cress section is obtained in turn by subtract- 
ing the fission and scattering cross sections from the 
total cross section, However, even in this energy 
region the total errors proceeding from lack of know- 
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Fig. 6. Integral distributions of reduced neutron widths for U®*5 (a) and Pu®® (b). Dashed curves: Porter- Thomas 


distribution for vy = 1. 
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ledge of the amplitude. of the nonresonance term, 
presence of interference effects, indeterminacy in the 
normalizing of the fission cross section, etc., turns 
out to be quite large and casts suspicion over 

any interpretation of the experimental data. Table 3 
gives the values of radiation widths of resonance 
levels for the nuclides yr. us. and Pu”, as mea- 
sured by various authors, Although the accuracy 
achieved in determining Ty is not particularly high, 
there is still room for certain conclusions, The ra~- 
diation width varies in a relatively weak manner 
from level to level in a single nuclide, and has the 
following values (given in millielectron volts): (33 4 
4 10) Mev for U™* (31); (42 + 8)Mev for U™ and 
(40 + 8)Mev for Pu”? [27]. It is indicated in one 
paper [20] that the value of 'y is constant since 

the process may proceed through several channels. 

In computing the parameters of those levels for which 
Ty cannot be determined directly, constant values 
are usually assumed for Ty. However, it is stated 

in [32] that a grouping of levels which must be ascribed 
to the effect of a spin dependence is manifest in the 


case of Eu”! radiation level widths. It is also clear 
Filter Photomultiplier 
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TABLE 4, Radiation Widths of Resonance Levels of 
Pu” as a Function of the Spin of the Compound 


Nucleus 
ry, Mev 
Eo, 1 3 
0,30 40+3 40+3 
7,83 39+4 41+4 
10,93 28+8 32+10 
41,90 39-+5 41+5 
17,6 37+5 40+5 
22,2 34+5 
41,5 46+9 1629) 
44,5 30-+10 42+10 
52,6 38-++4 


from Table 3 that fluctuations in Ty exceed, however 
negligibly, the indicated values for U™, 

The relative constancy of Ty may be used to ad- 
vantage [18] in finding the spin of the compound nuc- 
leus, where the statistical factor g differs greatly 
for the two possible spin states of the nucelus, as. for 
example, in the case of Pu (g = ¥, for J = 0 and 
g = % for J = 1, where J is the spin of the compound 
nucleus), This may be done with success for reso- 
nances at energies of several tenths of an electron volt, 
having an appreciable neutron width. By measuring 
Oo, Xo, I’ and it is possible to find two values 
of “rom the relations +Ty and oT = 
=4mk3T,. Furthermore, by requiring that Ty be 
constant for all levels of a particular nuclide, it is 
possible to choose one of the two possible values of g 
and thus to assign the spin J, The spin assignments 
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Fig. 8. Integral distributions of partial widths of Pu™® (a) and U™®* (b), Porter-Thomas distributions for v = 1, 
v = 2, and v = 4 indicated. 
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Fig. 9. Value of o, ¥E for U™*, averaged over several resonances. 


found in this manner for levels of Pu™ having E, equal 
to 41.5 and 52.6 ev show a value J = 0 for the first 
level and J = 1 for the second level [18] (see Table 4). 


The Partial Width. 

The partial width Iy of a level may be found by 
measuring the fission cross section and absorption cross 
section with the aid of the equation ~ ~ 


~ az (with the integration taken 


over the entire level). Gas scintillation counters have 
been employed with success in recent experiments to 
measure the fission cross section [35, 36]. The design 
of one such counter used is shown in Fig. 7 [29]. 
Fission fragments are recorded by the phototube mul- 
tiplier from the emissions of a noble gas, e.g., xenon 
or argon, which they elicit. Owing to the short 
fluorescent lifetime (e.g., T= 2 ° 10° sec for a mix- 
ture of A + N,), a rather large amount of the fission- 
able material may be placed within the counter 
chamber without increasing the background due to 
particles, 

Another technique for determining I's is based 
on the assumption that ',, remains constant for all 
levels of the nuclide. It is stated above that the 
mean values of [y are 33 + 10 Mev for U™, 42 4 8 
Mev for U™®, and 40 + 8 Mev for Pu™, Since the 
accuracy achieved in measuring total widths of levels 
is not too high either in most cases, significant errors 
arise when I; is found by this method, 

There remains another method for determining 
I's, based on the use ot the constant ‘value of T, and 
the number of neutrons per fission event v. This 
method consists of measuring the effective number of 
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neutrons Ve¢s at the energy corresponding to the middle 
of the resonance peak, and using the relatiop 


“eff ry 
9— eff. 


This method was first employed by Farley [37], Approxi- 
mately the same method was employed to determine 

the partial widths of Pu™ in other experiments [18, 

29}. In addition to possible fluctuations in the value 

of [',, this method allows for some arbitrary judgment - 
in assigning the value of Vegr. 

Figure 8 gives the integral distributions -of partial 
widths based on data reported in [29, 31]. Comparison 
with the Porter- Thomas distributions reproduced here 
shows that the fission reaction proceeds via a small 
number of channels (possibly one or two at the. most), 
The lower limit of the number of channels available 
to the reaction may be evaluated with the aid of the 
quantity 2nT;/D [2], v = 1.2 is obtained for 
by this approach, It is interesting to note that a 
periodicity was observed in the value of of JE averaged 
over many resonances, as a function of energy (Fig. 9), 
for U™ {4, 31]. This points to a correlation between 
I's and D over a broad energy range; the reason for this 
correlation is as yet obscure, 


Number of Neutrons Per Single Fission 
Event 
The number of neutrons per single fission event 
v may vary from level to level and may depend 
on the spin of the compound nucleus, for example. 
The region of energies extending from thermal ener- 
gies to 0.3 ev have been intensively studied by many 
authors; it is shown in several papers [38-41] that the 
value of v is constant over that range, within limits 
of 1-2%, The proof of the constant value of v (with- 
in limits of 3-4%) for resonance peaks of Pu™ may 
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Fig. 10. Counting device 
for measuring number of 
neutrons per fission event 
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counter with 114 mg Pu” 
specimen, 


be found in [42]. Figure 10 shows an experimental ar- 
rangement used in the associated research, The ar- 
rangement is capable of simultaneously registering 
counts of fission fragments in the gas scintillation 
counter and coincidences between pulses in this 
counter and in the two liquid scintillator chambers 
due to secondary fission neutrons, Table 5 presents 
the results of the measurements, 

It has been established [43-55] that, as the total 
excitation energy of the nucleus increases 1 Mev, V 
increases by 5-6% on the average. This is in agree- 
ment with calculations based on the concept of eva- 
poration of neutrons from the fragments. At the same 
time, theoretical estimates [1, 2] indicate that the 
first excited states of a nucleus found at the top of the 
fission barrier differ from each other by an energy 
increment of ~ 1 Mev. An improvement of 1-2% 
in the accuracy of the measurements reveals a possible 
grouping of energy levels according to the spin of the 
compound nucleus, — 


Effective Number of Secondary Neutrons 

The quantity V.¢¢ is related to the fission cross 
section O¢ and to the absorption cross section 0, by 
the equation 


Veff =V =%V 
A direct measurement of the energy dependence of 
Ver based on recording secondary fission neutrons [46, 
47] enables us to obtain several level parameters (par- 
tial and radiation widths, where the total width is 
known) and to verify the presence of interference bet- 
ween levels, This method is not quite as sensitive 
to the resolution of neutron selectors as the method of 
determining parameters directly from cross section 
measurements, since the value of UV. ¢¢ is far less 
sensitive a function of energy than the interaction cross 
section. Measurements are performed as a rule with 
specimens of large thickness; this significantly enhan- 
ces the effectiveness of the technique. The absence 
of background due to & particles, which puts severe 
limitations on the use of fission chambers, is to be 
noted, 


The scintillation detectors used to record the se- 
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TABLE 5, Results of Measurements of vy for Pu in 
the Resonance Region 


Neutron V,relativeunits {Measurements 

energy, jerror, 

ev 

Thermal 1,008. 1,6 
7,8 1,009 2,6 
10,9 0,972 1,9 
14,3 0,999 3,3 
14,6 1,001 13,9 
15,5 0,975 3,0 
17,6: 1,083 4,0 
22.2 1,006 3,6 
26,6 0,983 3,4 
59 0,954 6,1 
66 1,048 3,0 
75 1,041 3,4 
86 1,023 3,2 

440—1700 0,994 


condary neutrons were both zinc sulfide scintillator 
instruments [19] and liquid scintillation counters [29]. 
The latter features high efficiency (of the order of 8- 
10%), with a drawback being their gamma sensitivity. 
The qualitative agreement between the results obtained 
in these experiments shows that the use of the liquid 
scintillators. is fully justified. At the present time, 

the possibility has arisen of devising highly efficient 
scintillation detectors for neutron fission, with remar- 
kably low sensitivity to gammas. Several organic 
scintillators (stilbene and some liquids) have two 
fluorescing components differing significantly in fluo- 
rescent lifetimes, It has been demonstrated [48, 49] 
that the intensity of the slower-acting component varies 
sharply with the nature of the incident radiation; this 
makes it possible to distinguish neutron counts from 
pulses due to the background of gammas, 

The energy dependence of Vo¢¢ in the resonance 
region has been investigated most thoroughly in the 
case of Pu”, This is because of the special interest 
in this nuclide on account of its small spin (I = 42) 
and because Pu has the largest mean level spacing 
of all the fissionable nuclides, The value of V.¢¢ 
averaged over individual resonances has been deter- 
mined [50]; the detailed behavior of v ff within the 
range of the resonances has been studied [37, 23, 19]. 

Measurement of Ve¢s is attended by some diffi- 
culty in the region of small cross sections, because 
of scattering corrections, and the accuracy of the data 
for regions between resonance peaks is therefore not 
too high. To gain an independent determination of 
Veg in that case, the total cross section 0; and the 
fission cross section Of are often measured, with a sub- 


sequent calculation of the value of Ves from the for- 
mula 
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Fig. 11. Neutron energy dependence of the ratio V.¢¢/v in the resonance region for Pu, ©) direct 
measurements [29]; X) computed from cross sections [29]; ———) from data in [37]; @) from data 


reported in [30] (results renormalized by Vore/Yote T= 0.58% [19] for resonance at 7.8 ev energy); 
4) from data reported in [19]. 
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Fig. 12, Energy dependence of Var, for U™*, in the energy region from 0,01 
to 10 ev: @) preliminary data obtained in the MTR reactor; A) data from 
Brookhaven National Laboratory. The solid curve is plotted from the genera- 
lized formula, Arrows on the diagram denote the position of resonances. 
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where 9, is the scattering cross section, Figure 11 
reproduces the results of papers [19, 29, 50, 37], with 
those points calculated from measured cross sections 
also plotted on the diagram [29]. The value of Vogf 

for U™ has been measured over an appreciably narro- 
wer energy range than that for Pu, Figure 12 repro- 
duces the results of an experiment [23] performed in 
the neutron energy region from 0.01 to 10 ev. Detailed 
measurements of the energy dependence of Veg¢ for 
Pu™ and U™ show that Vere Varies greatly within the 


range of individual resonance peaks, This may be an 
indication of the presence of pronounced interference 
between resonance levels, However, as is evident 
from Fig. 12 where the solid-line curve indicates the 
results of calculations of the energy response of Vor, 
based on the generalized formula, only a qualitative 
agreement with empirical values is attained, 

The region of thermal energies, not physically 
delineated, belongs to the resonance region; appreciable 
variations in Veg¢ associated with the effect of reso- 
nances having different parameters are observed to take 
place in the immediate neighborhood of that region. 
Measurements of Uo ¢¢ in that region are of practical 
significance, and are pursued with great care, 

Figure 13 reproduces theresults of direct measure- 
ments of Vog¢ for Pu™ performed by different authors, 
as well as the energy dependence of U,r¢ computed on 
the basis of measured cross sections [27]. In the diag- 
ram, we have the quantity a related to Ve¢r by the 
following equation: 


1 Ca — OF 


The energy dependence of v,¢¢ measured in an ex- 
periment [19] using a thin sample of Pu is found to 
coincide practically with that calculated from the 
measured cross sections, The value of v.-- at the 
resonance at energy 0.3 ev is reduced by 20% in all of 
the papers mentioned (with a spread of roughly 1% in 
the results), compared to the value at thermal energy; 
at energies near 0,07-0,15 ev, a certain discrepancy 
is observed between the findings reported in different 
papers, 

A comparison of the results of measurements of 
Voge for U™® is presented in Fig. 14 [26]. Also indica- 
ted is the energy dependence of Vers, computed from 
cross sections, 


Conclusion 

Despite the comparatively high accuracy in the 
measurements, all aspects of interaction between 
neutrons and fissionable nuclei have not been success- 
fully described with sufficient rigor of analysis. 

The generalized formula [22] enables us to ex- 
plain some effects rather convincingly, but no exact 
agreement is obtained between calculations and ex- 
perimental data [24]. This may quite possibly be 
related to the assumption undertaken in deriving the 
formulas as to the presence of a single channel in the 
fission process; the results reported in [55] are of inte- 
rest from this vewpoint. A comparison of experimental 
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Fig. 13, Neutron energy dependence of o in the energy region from 0.001 to 1.0 ev for Pu: O) BNL data 
[41]; O) Harwell lab data; X) BNL data [41]; e) Hanford lab data [51]; @) Argonne Nat'l Lab, data [52]; 
4) Harwell lab data [53]; +) data reported in [46]; solid curve indicates calculations from cross sections. 
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Fig. 14. Energy dependence of « for U™ in the energy region from 0.001 to 1.0 ev: O) BNL data [41]; 


4) Hanford data; @) MTR data [54], 
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data was carried out in [24] and the generalized formula 
was derived by the use of a different approximation from 
that used in [22]: the interference between a smalt 
number of levels was taken into account, but no res- 
trictions of any nature were placed on the number of 
channels available to the reaction, Satisfactory agree- 
ment with experiment was attained; exact quantitative 
agreement was not, however. 

The question of the origin of the nonresonance 
term in the cross sections remains an open one as 
before. The assumption that it might find an auto- 
matic explanation with the aid of the generalized 
formula has not been validated [24]. In the thermal 
region, the presence of a nonresonance term has to 
be called upon to explain the overwhelming portion 
of the U™ cross section and a significant portion of the 
uU™ and Pu™ cross sections, It does not constitute 
the major portion of the cross section lying between 
resonances, In Figs. 15 and 16, we have the values 
of ogvE for U™ and Pu, the value of the nonre- 
sonance term is indicated by the dashed-line curve 
[31]. Direct measurements of Voge for U™* demonstrate 
that the fission cross section accounts for over 90% 
of the absorption cross section between resonances [23]; 
in the resonance region the ratio of /og falls off 
considerably, A similar pattern is observed for the 
cross section of Pu”. in the region between resonances 
with energies 0.3 and 7.85 ev, the ratio o¢/o, reaches 
0.9 [29], while of all the energy levels in the 0-30 
ev interval, only one of them, an unusually broad 
resonance at 15.5 ev energy (I'~ 1 ev) has a value of 
V ore approximating to that value, We may quite possibly 
be dealing here with a direct fission process, similar 
to the stripping process, for example. 

The study of the mass dependence of fission ' 
fragment yield for different excitation states of the 
nucleus, the study of fission gammas, the determina- 
tion of level spins, etc. make it possible to reach a 
better understanding of the fission process and of 
gamma emission. But only the first few papers [56, 
57] containing preliminary findings have appeared 
along this line. 
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Fig. 15. Value of ogvE for U™ in the 0-10 ev energy 


region. Dashed line indicates contribution from the 
Nonresonance term. 
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Dozens of resonance levels have been studied; an 
increase of the resolving power of neutron velocity 
selectors in the near future would produce parameters 
of a still larger number of resonances to be obtained, 
This is particularly important in clarifying the exact 
distribution pattern of partial and radiation widths, 
data regarding which are quite meager. 
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Determination of the half-life of Ra™ has great value because of the wide use of this isotope as a radioactive 
standard. in medical therapy, in industry, and in the preparation of neutron sources. This paper presents both 


data from the literature and original experimental work concerning the accurate determination of this quantity 


by a calorimetric method. 


Measurements of the caloric effect for three oquittbebam radium preparations yielded a Ra specific activity 
value of (3.71 + 0.02) - 10” dis/sec per gram of Ra™™, which corresponds to a half-life of 1577 + 9 ee. The 
calorimetric measurements were carried out in a double static calorimeter having a sensitivity ~ 2 - 10°‘ w. 

The radium content in the preparations was determined by direct weighing of the salt RaBr2, the purity of which 


was established by spectroscopic means. 


Discovered in 1898 by Pierre and Marie Curie, 
radium occupies a special place among the radioactive 
elements, For a very long time radium seemed to be 
almost the only radioactive substance being used for 
numerous and varied measurements and investigations 
in the field of radioactivity. If we disregard the appli- 
cation of artificial radioative isotopes, radium prepa- 
rations are used widely even today: in physics, for the 
preparation of laboratory neutron sources and strong 
a-particle sources; in medicine, for internal and ex- 
ternal radium therapy; in industry, for gamma defect 
detection and radioactive core sampling. 

Since Ra™ has a small disintegration constant, it 
is convenient to use it as a standard source, Up to the 
present time the majority of ionization measurements 
use radium-compound standards, Unfortunately the 
half-life T of Ra™ has not been determined with the 
required accuracy. This makes it difficult to correct 
for the natural decay of the numerous available primary 
and secondary international sandards, Accurate know- 
ledge of the half-life of Ra” is also desirable for mea- 
surements using differential ionization chambers (in 
which it is quite often used as a source of constant 
radiation), and likewise for other measurements and 
calculations, 

Much work has been devoted to the determination 
of the half-life of Ra™ and the closely allied mag- 
nitude of its specific activity (number of a particles 
emitted by 1 gram of Ra™ in 1 second) Z = 1.847 - 
-10°"% 7"! The majority of this work was done bet- 
ween 1908 and 1928, Table 1 shows the results of 
work which either determined or permitted calculation 
of the values of Z and T for Ra”, As can be seen from 
the table, the values of Z and T obtained by various 
workers vary within wide limits, which significantly 
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increases the errors cited by the individual authors, The 
deviation reaches 15-20%, * 

Within Table 1 the results are grouped by the me- 
thods wilined and by publication dates, Since the half- 
life of Ra™ is too long to permit carrying out direct 
measurements of the decay curve of this element by any 
means, various indirect methods have been used in the 
works cited. Generally, direct or indirect counts were 
made of the & particles emitted per unit time from a 
given weight of radium, i.e., in actuality the authors 

In 2N 
r= 

“dt 


use the relationship 


Measurement of the specific activity of Ra” j 
associated with specific experimental difficulties; on 
the one hand, with the absolute « count, and on the 
other hand, with the quantitative determination (weigh- 
ing) of the measured amount of radium. Each of the 
methods enumerated in the table has, as a rule, its 
own characteristic shortcoming. Thus for example, 
methods 1 and 5 are associated withthe difficulty of 
obtaining a thin homogeneous layerof RaC and of the 
quantitative evaluation of its activity, and also with 
errors dependent on the relatively rapid decay of the 
preparation itself and on inaccurate knowledge (in the 
case of method 5) of the ionization constants. The 
results of all work using the Ra™ accumulation rate in 
a solution of ionium (method 2) apparently contain 
some sort of systematic error which lowers the values 
obtained by this method. Perhaps this may be explained 
by saying that the minerals used in this work were 


* The results for Z and T in Table 1 have been recom- 
puted using more accurate values of physical constants, 


TABLE 1. Specific Activity Z and Half-Life T of Ra” From the Data of Various Authors 


Year of |_. Values obtained 
Authors Literaturel7 10% dis / sec Method used 
rs 
Citation |y 4 g Raz y 

| J. Danysz and W. Duane 1912 {1 3.15 1860 1. Measurement of the 

E, Rutherford andH.Geiger | 1914 [2] 3.45 1695 

H, Iedzeyowski 1928 [3] 3.48 1680 tira by a given quantity 

H, Braddick and H. Cave 1928 (4] 3.65 1605 of Ra2”, 


B. Boltwood 1908 [5] 2.94 1990 2. Determination of the 
_B, Keetman 1908 [6] 3.22 1820 | Ra™ accumulation rate in a 
St. Meyer and E, Schweidler} 1913 {7] 3.36 1740 solution of ionium (Th), 
E, Gleditsch and E, F8yn 1916 [8] 3.17 1845 
» 3,54 1655 
R. Lawson and St. Meyer 1916 [9] 3.36 1740 
E, Gleditsch 1919 [10] 3.46 1690 
3.54 1655 
E, Gleditsch and E. Fdyn 1932 (11) 3.46 1690 
E. Gleditsch and E, 1935 3.46 1690 
J. Dewar 1910 [13] 3.60 1625 8. Determination of the 
3.87 1510 volume of helium liberated 
B. Boltwood and E, Ruther- | 1914 [14] 3.49 1675 from a given quantity of Ra”, 
ford 
L. Wertenstein 1928 [15] 3.49 1675 
P, Gunther 1939 (16] 3.65 1605 
E, Rutherford and H. Geiger | 1908 [17] 3.57 1640 4, Direct measurement of the’ 
V. Hess and R. Lawson 1918 [18] 3.72 1575 number of & particles emitted 
H, Geiger and A. Werner 1924 [19] 3.40 1720 by Ra™ or by an equilibrium. 
H. Geiger and A, Werner 1924 [20] 3.48 1680 Ra preparation. 
F. Ward, C. Wynn-Williams} 1929 [21] 3.66 1600 
and H. Cave 
P, Kohman, D. Ames, and 1949 [22] 3.61 1622 
J. Sedlet 
W. Sebaoun 1956 [23] 3.62 1617 
I, Curie and F. Joliot 1928 [24] 3.65 1605 5. Measurement of the ion 
current produced by © par- 
ticles from a Ra” preparation. 
V. Hess 1912 [25] 3.77 1555 6. Calorimetric method of 
S. Watson and M. Henderson} 1928 [26] 3.64 1610 measurement. 
L. Meitner and W. Ortmann} 1930 [27] 3.68 1590 
G. V. Gorshkov and N, S. 1956 [28] 3.72 1575 
Shimanskaya 
Current work 1957 3.71 1577 


insufficiently aged or tha: the equilibrium between 
radium and ionium was disturbed (for reasons associated 
with the conditions of their preservation in the earth's 

crust). 


The determination of Z and T by helium accu- 


mulation from radium preparations (method 3) is asso- 
ciated with the difficulties involved in separating and 
quantitatively measuring small gas volumes, 

The values of Z and T shown in the latest tables 
of nuclear constants as the most trustworthy were ob- 


tained by the direct counting of & particles emitted 
by radium preparations (method 4), Precisely this 


method was used in the most recent references[22]- 


[23]. The authors of [22] used radium preparations in 
which they determined the amount of radium present 

by direct weighing of the natural radium in the aliquot 
portion, The determination of the activity of individual 
preparations was carried out with a plane-parallel 27 a- 
chamber, A large part of the radon and of the short- 
lived decay products was removed by a constant flow 


of gas over the sample. The measured results underwent 
a rather complex treatment, taking into account the 
difffusion of radon and the loss of recoil nuclei. The 
final results, obtained by averaging the separate results 
of a large number of trials, gave for Z the value 3.61 

dis/sec g Ra™, This corresponded to a 
half-life of 1622 years, The probable error for this 
result was taken as 0.5%, It is precisely on the basis 

of these results [22] that the majority of recent workers 
computing the number of disintegrations in 1 g Ra” 
have started to use (perhaps, in some cases, prematurely) 
the value of 3.61 - 10” in place of the value Z = 3.7- 

- 10° which was adopted by the International Committee 
in 1930. 

A close value of Z was obtained and published in 
1956 [23]; the author used a scintillation spectrometer 
to measure an equilibrium preparation of Ra” enclosed 
by a thin mica window and obtained values Z = 3.62 ° 
- 10" dis/sec * 1 g Ra™ and T = 1617 years (+ 0.7%). 
The small-geometry o-particle detector used in this 
work permitted the measurement of relatively strong 
sources of the order of 0.15-0.20 uC and eliminated the 
necessity of difficult corrections for o-particle rebound 
from underlying material. The possibility of discrete 
energy measurements for the & particles from Ra” and 
its daughter products, and likewise measurements of the 
time constants of equilibrium Ra” preparations appear 
to be essential advantages of the method of [23] in com- 
parison with the method of [22]. 


On the other hand, the indirect determination (by 


y radiation) of the quantity of Ra present in the 
preparations in question could be a source of signifi- 
cant errors. 


The essence of the calorimetric method of deter- 
mining T for radioactive isotopes ( method 6) consists 
in determining the caloric effect Q of a given weight 
p of a given isotope. In this case we obtain from the 
basic radioactive decay equation 

No 
T= In 2N 
dt 
where No is Avogadro's constant; A is the atomic weight 
of the specific radioactive element; E is the energy 
released within the calorimeter by a single decay event. 
Since E is usually calculated with sufficient accuracy 
for many isotopes from their known decay schemes, 
a single weighing and a single calorimetric measure- 
ment is sufficient for an equivalent value of T, because 
both the operations can be carried out quite accurately. 
It is true that this requires very high chemical purity 
of the preparations with respect to any other radioac- 
tive or inactive elements, since small admixtures of 
these, and especially of heavy elements, can distort 
the results to a significant degree. It is evident that 
this method is only applicable to long-lived radio- 
active isotopes (T> 10-10° years) for which it is 


possible to obtain a weighable quantity of pure isotope. 
Such a method was used to obtain the disintegration 
constants of [30] and [31]. Several authors 
have used this me od to obtain the specific activity 
and half-life of Ra . They measured the caloric 
effect of a given quantity of Ra™ freed of its short- 
lived daughter products [25]; likewise, Ra” in equi- 
librium with its products [26, 28]; or Po” accumulating 
in the radium preparation [27]. 

It is to be regretted that in all these experiments 
the determination of the radium content in the mea- 
sured preparations was carried out by indirect methods 
(usually be means of y measurements), Furthermore, 
the majority of these measurements was of insufficient 
accuracy with respect to that of the calorimetric 
measurements (1-4%); this was likewise true for the 
calculated values of E. 

Attempts were undertaken in 1954 to evaluate the 
half-life of Ra by means of calorimetric measure - 
ments of the caloric effect of two equilibirum prepa- 
rations of radium, the radium content of which was 
determined by ionization measurements [28], Correct- 
ing for accumulation of RaD, RaE and Po”? in the pre- 
paration, for absorption of y rays in the salt, ampule 
walls and the calorimeter itself, and approximating 
a correction for the difference in absorption and self- 
absorption of y rays in the preparation and in the stan- 
dard used for ionization measurements, showed that 
the value of the caloric effect of « andB rays per 
unit weight of Ra” in equilibrium with its short-lived 
daughter products, through RaC inclusive, was 130.0 + 
+ 1.5 cal/hr + 1g Ra™ [28]. Using a calculated value 
of the total decay energy E = 25,351 Mev/dis gave 
the values Z = (3.72 + 0.06) - 10” dis/sec - 1 g Ra™ 
and T = 1575 + 25 years. 

This evaluation of the .Ra”™ half-life could not 
claim better accuracy in view of the indirect deter- 
mination of the radium content of the preparation (the 
error of this determination was taken as 1.2%), There- 
fore, in the present work T in seeking a more accurate 
determination of the Ra™ half-life, we used prepa- 
rations in which the radium content was deter mined by 
direct weighing of radium bromide (Ra Br2) and whose 
purity was checked by spectroscopic means, Such pre- 
parations were prepared by the usual technological 
method used in radium production, but with supplemen- 
tary crystallization of the bromide. The preparations 
of Ra-I and Ra-II underwent ten supplementary crys- 
tallizations; the preparation of Ra-III, fifteen crystalli- 
zations. 

Table 2 shows the results for weighed RaBr salts 
of these preparations, The absence of water of crys- 
tallization from the RaBr, salts was guaranteed by in- 


tensive heating (to 200° C)before sealing the ampules, 


TA short preliminary account of the results of this work 
was published in the JETP 34, 576 (1958); Soviet Phy- 
sics JETP 7, 397 (1958), 


TABLE 2. Characteristics of Radium Preparations 


Internal ampule 


Dimensions of glass ampules, 
mm 
internal external 
a ampule ampule 
Date 
ampule 
sealed 


Height:to which ampule 


was filled, mm 
Weight of RaBr,, mg 
ity of 

by industrial 
Relationship of 1 mg- 
equiv of preparation to 
RaBr, by weight 


by VNIIM 


Ra-I 


July 4, 1955 


70 


o 


274,8] 2,08 | 158 | 157,0+0,7 | 0,5742 


Ra-Il 


Sept. 17, 1955 6,7 


71 


1,47 175,6£0,9 | 0,5713 


Ra-II] }Oct. 12, 1955|°4,4 | 3,6 | 6,35 | 5,35 


71 


18 1,45 147,7+0,7 


0,5760 


In order to verify that the stoichiometric equation 
holds constant for the salts of RaBr, during the making 
__ of the radium preparations, one of the preparations 
(Ra-III) was converted to radium sulfate and the re- 
weighed after completion of all the basic measurements, 
There was less than 0.5% difference in the weight of 
radium in the bromide and the sulfate. The stoichiomet- 
ric constancy of the RaBr, was also shown by special 
experiments using barium bromide, Several weighed 
samples of BaBrg were roasted to 600° C, and the bro- 
mine was precipitated with silver nitrate. The weight 
ratio of Br/ BaBrg was the same for roasted and unroasted 
BaBr, samples to within an accuracy of 0.1%, 

In addition to the results of weighings, Table 2 
contains other data and parameters characteristic of the 
various preparations; data concerning the sealing and 
internal and external dimensions of the glass ampules, 
and sample density values of the radium bromide 
estimated on the basis of ampule filling height. The 
table also shows the results of ionization measurements 
of these preparations, produced by industrial methods 
and by the All-Union Institute of Metrological Scien- 
tific Research (VNIIM). The last column shows the 
weight relationship of 1 mg-equiv of each preparation 
J (as measured at VNIIM) to radium salt. Since this 
relationship must be equal to 226.05/ 385.88 = 0.5858 
for a pure radium preparation, the lower values actually 
obtained show that there were either significant im- 
purities in the radium bromide or else that there was a 
decreased radium content in the preparation when 
ionization measurements were taken. 

The calorimetric measurements were made in a 
static double calorimeter which is used in the Radiation 
Institute of the Academy of Sciences SSSR for mea- 
suring the activity of radioactive preparations [32]. 

The operating principles of such a calorimeter can be 
understood by examining Fig. 1, which shows its prin- 
cipal parts. Two identical calorimetric cylinders A 


and B are arranged inside a heavy jacket C, whose 
temperature is kept constant. Into one of the cylinders, 
for example A, there is placed the heat source to be 
measured, which heats the cylinder, The other cylinder 
B serves as a datum temperature point of temperature 
Ty. When thermal equilibrium is established between 


Fig. 1. Drawing of principal parts of the double 
static calorimeter, 


the surfaces of these calorimetric cylinders, the tem- 
perature difference T-T, can be determined; because 
of the construction of the calorimeter, this difference 
is due only to the caloric output of the source, Using 
several simplifying assumptions, one can show that 
this temperature difference is related to the source 


output W and the external surface area of the calori- 
meter S as follows; 


vj W 
r—T,=a-. 


(1) 


where the coefficient a, usually called the heat 
emission coefficient or the surface heat conductivity 
coefficient, depends on the nature of the calorimeter 


surface and onthe thermal characteristics of the surround- 
ing medium, 
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Fig. 2. Drawing of the calorimetric equipment; A and B are the calo- 
rimetric cylinders; C and D are copper shields; E is a thermostat, Elec- 
trical diagram: g, and g» are galvanometers; Ry and Rp are KMC-6 
100-kilohm resistance boxes; rg is a kor 10-ohm standard resistor; ry 
is a 20-ohm wire resistor; rp is a 500-kilohm wire resistor; PPTV-1 is 
a high-resistance potentiometer; mA is an MC-16 milliammeter; Vo is 
a standard Weston cell; V is a storage cell or filament battery of 2-20 


Equation (1) is obtained from another more general 
equation (2), which specifies the temperature diffe- 
rence in such a calorimeter at an arbitrary time t; 


9 
T-Ty=a%(1-e 


(2) 


where c is the heat capacity of the system consisting 
of the calorimetric cylinder and the preparation being 
measured, It follows from the equation, moreover, 
that the time required to establish thermal equilibrium 
and consequently the duration of the individual mea- 
surements depends on the heat capacity of the calo- 
rimeter and on its surface area, We call the quantity 
Xo = S/ac the calorimeter inertial constant, which is 
determined in every case from the experimental curve 
for heating of the calorimetric cylinder; 


T—T,=] (0). 


In our calorimeter A» was ~ 0,04 min”}, con- 


sequently, each measurement lasted 2-3 hr. The other 
principal characteristics of the calorimeter used in the 
measurements described below were as follows; the 
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volts; Ky, Ka, Kg and Ky are switches. 


internal and external diameters of the calorimeter 


cylinders were 18 and 16 mm respectively; their height 
was 85 mm. These cylinders were placed inside a 
massive copper shield C mde up of two thick-walled 
(wall thickness ~ 6 mm) copper cylinders joined to- 
gether by one of the forms. The space between the 
calorimeters and the shield was 3 cm. 

The inside surface of the shield and outside surface 
of the calorimetric cylinders A and B were silvered 
and polished to reduce the heat loss by radiation. The 
shield C was contained in a massive copper box D which, 
in turn,was placed in a thermostat in order to eliminate 
rapid thermal fluctuations, Figure 2 shows a diagram 
of the calorimeter used for our measurements, The 
temperature of the room in which the calorimetric 
equipment was used was maintained constant to within 
+ 0.2°C, This was done by a very simple device con- 
sisting of a relay operated by a contact thermometer 
and one or two conventional air conditioners, Further- 
more, a “compensating” heat capacitance was placed 
in the second “neutral” calorimetric cylinder; this 
significantly reduced the influence of the temperature 
variation of the surrounding medium. 

The temperature difference T-Ty of the calorimetric 
cylinders was measured thermoelectrically with a copper- 
constantan thermocouple, Figure 2 also shows the 
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general electrical circuit diagram of the calorimeter 
apparatus, including the calibration circuit and the 
actual measuring circuit. 

The key K, was made entirely of copper, which 
eliminated the possibility of generating parasitic 
thermoelectrodynamic forces within it, The cross- 
shaped connections in this key eliminated the influence 
of thermoelectric and contact potential differences in 
the galvanometer leads, _ 

The key K, permitted us to establish the null 
position of the galvanometer and to control its sensi- 
tivity by means of a standard Weston cell. This sen- 
sitivity and the whole measuring range could be 
changed by a series-connected resistance box Ry. The 
calorimeter was calibrated electrically by means of a. 
calibrating coil made of manganin or constantan wire. 

The current intensity in the calibrating coil was 
determined by the voltage drop across the standard 
resistance rg connected into the circuit of a high-re- 
sistance potentiometer of type PPTV-1. During the 
calibration of the calorimeter we determined the 
dependence of the galvanometer deflection on the heat 
output produced by the electric heating element. 
Figure 3 shows the calibration curve for the thermal- 
effect region which is of interest to us (18-24 cal/ hr), 


d.mm, 
130 
jon > “4 
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Fig. 3. Calibration curve. 


This region contains the values of thermal energy 
emitted by our radium preparations. As can be seen 
from the curve, the sensitivity of the calorimetric 
apparatus is equal to 2° 10 w/mm., Using the same 
galvanometer and a voltage constant of 3 - 10™' v/mm- 
-'m, it was possible to have a quite stable operation of 
the calorimetric apparatus with an increased sensitivity 
equal to 2-4 - 10-5 w/mm., In this case, several (5-10) 
lead-constantan thermocouples were placed in each 
calorimetric cylinder; for electrical isolation, the 
thermocouples were placed in thin glass capillaries 
and fastened with shellac into corresponding apertures 
in the bottom of the cylinders, 

Numerous special trials were made to determine 
the characteristics of the calorimetric apparatus, to 
determine the inertial constant A, to study the influence 
of the surface reflectivity of the calorimetric cylinders, 


and to study the heat transfer in the calorimeter. We 
also investigated the influence of the space between 
the calorimetric cylinders and the shielding on the 
calibration curve. It was established that the equilib- 
rium temperature difference T-T) does not depend on 
the heat capacity of the preparation being measured, 
nor on the heat source distribution within the calori- 
metric cylinders, nor on the heat transfer conditions 
within the cylinders, To eliminate completely the 
last possible influencing effect, the calibrating heating 
coil reproduced as closely as possible the heat emission 
and heat exchange of the measured preparations 
themselves, Special attention was given to establishing 
the accuracy with which the various thermal effects 
could be measured with this apparatus, It was shown 
that under the existing thermostatic conditions it was 
possible to measure the thermal effects, which exceeded 
5-10 cal/hr, with errors ~ 0.3-0.5%, - The error was 
caused generally by inaccuracy in reading the galva- 
nometer deflection and by insufficient stability of the 
galvanometer sensitivity. The error in determining the 
thermal output generated in the calibrating heating 
element, as taken for the calibration curve, was less 
than 0.1%, 

Three series of calorimetric measurements were 
made on radium preparations, with about two months 
time between series, Table 3 shows the results of 
various calorimetric measurements of Q, the corres- 
ponding values Qp (magnitude of the thermal effect at 
the time of sealing up the corresponding preparations), 
and the average values Qy for each preparation. When 
Qowas calculated, there was taken into account the 
thermal effect Aq; for RaD, RaE and Po”? which arose 
because of the lapse of time from the moment of seal- 
ing the ampules until the measurements were made. 

In making these corrections, the following values were 
used from the latest available data; for the o-particle 
energy of Py”, 5.306 Mev [33]; for the average energy 
of the B spectrum of RaE, 320 kev [34]; for the disin- 
tegration energy of RaD, 48.5 kev [35]; the correspond- 
ing half-lives T were taken as 138.4 days [36], 5.013 
days [37], and 19.4 years [38]. 


Next we calculated the thermal effect Ady ° 
which depends on the absorption of a certain part of 
the y radiation both in the salt preparation itself (self- 
absorption) and in the glass ampule (with its protecting 
brass container 0.5 mm thick), and in the several 
millimeters of copper comprising the calorimetric 
cylinder. It was assumed that the effective mass ab- 
sorption coefficient 4 was 9, for the y radiation from 
the equilibrium preparation of Ra in RaBry salt, 
equal to 0.65, These computations showed that 1.6- 
1.7% of the total y radiation was absorbed in the pre- 
paration themselves, 

The absorption of y radiation in the glass, brass 
and copper was evaluated from the basic experimental 
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TABLE 3, Calorimetric Measurements of Radium Preparations (cal/ hr) 


Preparation 
Ra-I Ra-II Ra-JII 
Experiment 


Series 1. 12-15 


21 

Dec, 1956 0,144 

21 


Series 2. 6-9 
Feb, 1957 


21,08 20.90 | 23.42 
Series 3, 23-30 24.12 0.183 20.94 23.42 
Mar, 1957 21,12 20. . 
° 24.12 20.94 


23.23 19.63 
23.23 19.60 0.159 | 19.44 
| 0.190 23.33 19.72 19.56 


20.96- 


Average +0.03 


19.55-+- 
£0.03 


y -absorption curves for these materials with a radium 
source, These curves showed that the correction for 
absorption in glass was 1%, while for the copper calo- 
rimeter walls and brass protecting container it was 3.5%, 
Since y radiation comprises only 7% of the emitted 
energy [39], the increase in thermal effect to correct 
for y -absorption is only 0.43%, The next to last 
column of Table 4 shows values of Q, 28 = Qo - Ady 
i.e., the thermal effect of the corresponding radium 
preparations depending only on the easily absorbed 
radiations (y rays, recoil nuclei, B particles, and con- 
version electrons) and on a major part of the x rays 
accompanying the conversion process, In order to 
obtain we also calculated the correction 
for the braking radiation not absorbed in the calori- 
meter; in our case this was negligibly small ( ~ 0.03%), 

The value of Qq.g was used to calculate the 
thermal effect q'y,g per unit weight (for 1 g Ra™), 
The corresponding values are given in the last column 
of Table 4, 

Moreover, the quantitative spectroscopic analysis 
of all the preparations was taken into account. This 
analysis, done by the copper spark method, showed 
that all three preparations contained only insignificant 
amounts of elements of the second group of the perio- 
dic system, The total impurities in no case exceeded 
0.2 weight per cent, The o-active impurities found 
for each preparation were taken into account in calcu- 
lating Qgig and q'g4g. It can be seen from Table 5 
that this correction was almost identical for all three 
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preparations, The deviation from the average q', +B 


value of 129.7 cal/hr per g Ra” did not exceed 
0.4%, The error for this value was taken to be 0.5%, 
This was composed of calorimetric measurement error 
(0.4%), weighing error (0.2%) and errors in estimating 
the various correction coefficients (0.1%). 

The spectroscopic analysis error was taken to be 
50%, which produced a 0.1% error in the final value of 

Comparison of the value p = 0.5858 P (1-c) with 
the corresponding value for 1 mg-equiv of preparation 
J (Table 4) shows a marked excess (on the average 
order of about 2%) of the first value over the second**, 
Furthermore, it was evident during the ionization mea- 
surements that radium content in the preparations had 
decreased appreciably, This was to be expected, since 
the primary standard XI containing 14.45 mg radium 
used at the VNIIM for y measurements was (1) signi- 
ficantly smaller than our preparations, and (2) was not 
enclosed in the second glass ampule. The latter cir- 
¢The value of ('y4g that we obtained for the thermal 
effect of a and 8 radiations per g of Ra™° is some- 
what less than the result q'q4, =130.5 cal/hr per g 
rd** [39], Another paper [41] reports a microcalori- 
metric method of comparing three international radium 
standards; for the thermal effect of the radiation from 
1 g Ra”®, filtered through 1.84 mm of lead, the value 
130.2 cal/hr is obtained, which corresponds to a value 
of q 128.3 cal/hr per g 


ey 20.92 | 23.20 23.06 | 19.66 19.54 
2108 | 23.26 23.12 | 19.70 19.58 
21.09 | 23.29 0.144] 23.15 | 19.72 | { 9-116] 19.60 
21.41 | 23.42 23.28 | 19.92 19.80 
21.16 | 23.53 23.37 
| 
20.34 20.65 | 23.09 22.92 | 19.57 19.43 
: 21.44 0.164 | 20-95 | 23.11 22.94 | 19.67 0.1414 | 19-58 
| 21.07 : 20.91 | 23.51 0.168 | 23.34 | 19.77 : 19.63 
21.05 20.89 | 23.47 93.30 | 19.76 19.62 
: 23.37 23.20 
| 


TABLE 4. Thermal Effect of 1g Ra” 


| Thermal effect, 


| mcal/br fectofig Ra%6 
Weight of | Quantity of | Weight of mg-equiv Ons | 
RaBro, g additives. |Ra™, mg of prepara- Co (Qe4g)/ P, 
%o ion cal/hr «g 


Thermal ef- _ 


Ra-I 
Ra-ll 


160.7 
178.8 


157.0 
175.6 
147.7 


20.96: 
23.19 
19.55 


20.88 
23.10 
19.48 


129.20 
129,90 
129,95 


Average 4g YP=4 a4 = 129-70 


TABLE.5, Average Disintegration Energy Emitted in One Disintegration by Elements of the Ra?” Series 


Element 


Average energy per disintegration, kev 


a radiation, 


B particles, 
| 


' Conversion 
electrons.and 


x-rays, 


y quanta, 


+E3 +E, 


+Ey 


Ra 
Rn 
RaA 
RaB 
RaG 
RaC’ 
Rac” 


S| 2811 | 
OS 


7825.3 
1.2 


24 376.5 


+y 


870,2 


=138,34 cal/hr per C Ra, 


88.6 


1784.14 


25 335,3 


Note: Eq+g=25335,3 Mev/dis q a+ p= 129,24 cal/hr per C Ra2®, B+y =20119.4 Mev/dis Ra ; 
q 


27 119.4 


close to the experimentally observed value, 


It was possible to calculate the half-life of Ra 


directly from the related value of q'y, = 129.7 


cal/hr per g Ra”, 


Table 5 shows values of the energy Ej emitted in 
one disintegration for the various elements involved 
in radium decay, since the distribution of this energy 
among the various radioactive radiations is important 
for the calculation of the thermal effect, (for estima- 


cumstance, taking into account the variations of the 
wall thickness of the inner ampules, lowers the radium 
content of our preparations by 0,5-0.7%, Our calcu- 
lations also showed that there was a significant diffe- 
rence in y self-absorption between 15 and 150 mg of 
radium preparation, equal to 0.8%, Thus the calculated 
value of the difference between J and p (1.3-1.6%) is 


tion of the radiation fraction absorbed in the calori- 


meter), Table 5 also shows the values of the energy 
Ei, Which depend on & radiation (taking nuclear re- 


coils into account), 8 radiation, conversion electrons, 
and x-rays accompanying the y-conversion process. 


The table also shows the total y -radiation energy for 


each element of the radium series, calculated from the 
basic spectroscopic data concerning the relative inten- 


* «A similar difference between 1 mg-equiv of prepa- 
ration and the weight of radium contained therein oc- 
curred in the RaSO, radium preparation prepared sub- 

sequently from Ra-III (see Table 2). Taking additives 
into account, the weight of radium for this preparation 
was 148.2 mg; the weight of 1 mg-equiv preparation 
J was 145.3 mg-equiv; consequently, the difference 


between these two values was about 2% as mentioned 
before. A more particular problem concerning the sys- 


tematic lowering of the radium content in radium pre- 
parations in ionization measurements, which depends 


on the different conditions for absorption of y radiation 


in the government radium standards of the USSR and 


in measured preparations, is discussed in [42]. 


274.8 0.17 
Ra-lIII 256.3 0.18 149.9 
E 

4854.4 5.3 5.8 4859.4 4865.2 | 

5985.3 = 0.4 5585.3 5585.3 
6110.1 6110.14 6110.1- 

— 70.4 179.4 290.1 469.5 

43.2 1598 664.9 2262.9 
7825.3 — — 7825.3 
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sity of the y lines of Ra, RaB, RaC, and RaC’, and 
experimental values obtained in 1935 [39] for the 
thermal effect of y radiation from an equilibrium pre- 
paration of Ra™ (9,1 cal/hr per g Ra”), Table 5 
differs only slightly from the corresponding tabular 
values of E; [29]. The above-mentioned references 
furnished data on the average energy of the RaB and 
RaC spectra, data on the number of conversion electrons 
from RaB and RaC, and some changes in other values. 
For the energy of easily absorbed radiation from 1 
Curie. of Ra” in equilibrium with its short-lived 
daughter products there were obtained the values 
= Mev/ dis and = 129.24 cal/hr 
per C Ra™ (1 C = 3,7+ 10” dis/sec). The error in 
the values of Ey,g and qq4g was obviously no greater 
than 0.3%, since the value of the o-particle energy, 
which comprised a major part of the energy, was found 
with an accuracy no worse than 0,1 or even 0.05%, 
The error in determining the average B energy and the 
energy of the conversion electrons was significantly 
greater. This error was in fact not less than 3-5%, but 
since these values comprised only 3.8% of the total 
energy, the over-allerror due to these inaccuracies 
did not exceed 0.2%, 

Using the above value dq.g = 129.24 cal/hr 
per C Ra™’, we obtain for’Ra?™ the value Z = (3.714 


+ 0,02) - 10% a-particles/ sec per g Ra?™ and T = 15774 


+ 9 years. 

“Our values of Z and T ake close to those adopted 
by the International Committee (Z = 3.70 + 10" dis/sec 
per g Ra*™, T = 1580 years), but differ significantly 
from the values shown in [22] and [23] (Z = 3.61-3.62)> 
- 10” dis/sec per g Ra”, T = 1617-1622 years). In 
order to explain such differences. it would be interest- 
ing to repeat the determination of these important 
values by the methods specified in the references and 
by other methods, 

The authors wish to express their thanks to Yu. S, 
Martynov, who participated in the calorimetric mea- 
surements and in certain calculations, We also thank 
E. K, Smirnov for dependable preparation of the radium 
samples, I. E, Starik and D. M. Ziv for valuable 
suggestions and discussions concerning the work, 
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SECTIONALIZED REACTOR SYSTEMS* 


Letters to the Editor 


Translated from Atomnaya Energiya Vol. 7, No. 5, pp. 456-457 


November, 1959 
Original article submitted January 6, 1959 


Sectionalized reactor systems represent a reactor 
combination consisting of a high-energy loaded cri- 
tical reactor (a source of preliminary neutron irradia- 
tion) followed by noncritical sections — neutron mul- 
tipliers. At the boundary between the critical reactor 
and the first subcritical section and also at the boun- 
daries between the critical sections there are neutron 
valves which permit the passage of neutrons in one di- 
rection only. The neutron filters consist of a "black", 
for thermal neutrons, layer of vy, a layer of cadmium 
(it is more advantageous to use a layer of thorium or 
strongly depleted uranium) and a layer of retarding 
material. The uranium layer is directed toward the 
critical reactor or to a preceding subcritical section. 
The coefficient of multiplication of the subcritical 
section is determined by the degree to which it is 
subcritical. 

The overwhelming majority of neutrons emerging 
from the subcritical section cannot pass through the 
neutron valve in the reverse direction since they are 
retarded within the retardation zone and are absorbed 
in the absorption layer. 

Independently of the data presented in Borst's 
paper t [1], the neutron valve described for separating 
the subcritical sections into convergent systems must 
exclude the possibility of an uncontrolled increase in 
the reactor power should the reactor break down, It 
is expedient, in order to increase the effectiveness of 
the neutron valve, to place a thin layer of water or of 
some other water-containing substance near the cad- 
mium layer, Inasmuch as the scattering cross sections 
of hydrogen for the fast neutrons of the separation 
spectrum are significantly less than for the retarded 
neutrons, the hydrogen layer permits the passage of 
the fast neutrons in the forward direction and retards 
the passage of the slow and thermal neutrons by means 
of the absorption layer. Thus the layer of water during 
the one-sided fast neutron irradiation acts as a valve 
and substantially decreases the number of neutrons 
absorbed by the absorption layer. Calculations carried 
out by V. V. Orlov have shown that for a water layer 
~ 2 cm the coefficient of neutron passage through the 
valve may reach 0,9. It is necessary to note that the 
neutron valve is useful in "separating" the various 


types of reactors into the fast and intermediate neutron 
types. The method of developing energy by the use 

of subcritical sections permits us to basically increase 
the “depth” of burning of the fissionable isotope. In 
order to illustrate this let us present the necessary 
parameters of a sectionalized reactor system having 
the same grid as the reactor of the first atomic elec- 
trostation [3]: 


The coefficient of multiplication of the grid. . .1.38 
Total neutron leakage, %..... 28 
Geometrical neutron leakage through the 
Percentage of critical reactor neutrons pass- 
ing through the valve to be multiplied in 
the subcritical 10 
Effective coefficient of multiplication of sub- 
critical section necessary for the section to 
operate with power equal to the power of 


We may obtain an effective coefficient of mul- 
tiplication of a subcritical section of 0.9 by decreasing 
the load by taking into consideration the influence of 
the form of the assembly of the subcritical reactor 
(ring section) by loading the grid with depleted heat- 
emitting elements. In this process the “depth” of com- 
bustion of the fuel may be increased two to three times, 
inasmuch as the reactive reserve to be used up in the 
course of the reactor usage in the atomic station is 0.03. 
The fact that a subcritical section can emit energy when 
Kefg < 1 permits us to use a sectionalized reactor sys- 
tem to obtain the limiting high flow of thermal neutrons. 
Actually, for the same values of the parameters of 
neutron multiplication and neutron loss which exist in 
a critical reactor, the coefficient of utilization of ther- 
mal neutrons in subcritical sections can be decreased 
according to the value of k,-- due to the amount of 
construction materials; this gives us the possibility of 


* Presented at the Physics Section of the Scientific 
Soviet of the Physics Institute of the Principal Directo- 
rate on the Use of Atomic Energy during the Soviet 
Conference of Ministers of the USSR in 1957. 

T Neutron valves are also described in paper [2]. 


increasing the flow of thermal neutrons by about 3 times, 


for example, for the same loading. 
In order to obtain a limiting high neutron flow, it is 


advantageous to apply a “descending” sectionalized reac- 


tor system, in which the loaded starting critical reactor 
of the ring irradiates with neutrons from the adjacent 
inner subcritical section, which in turn irradiates the 
adjacent central section. In such a “descending” sec- 


tionalized reactor system the neutron valves permit neut- 


rons to flow from the periphery to the center. The app- 


lication of sectionalized reactors is very much in the fore- 


front of the prospective picture for realizing pulsed ope- 


ration in which there is a large neutron flow in each pulse, 


In conclusion, I wish to thank A. K. Krasin, I. I. 
Bondarenko and L. A, Gerasev for their interest in my 
work and for their valuable advice. 
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EXPERIMENTAL STUDY OF THE DENSITY OF HEAVY 


WATER* 
S. L. Ribkin 
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Original article submitted December 19, 1958 


The density of heavy water was studied with a 
constant-volume piezometer up to temperatures of 
300° C and pressures up to 100 kg/cm*, The piezo- 
meter of known volume was filled with the liquid being 
studied at a known temperature and at atmospheric 
pressure, i.e., under conditions for which the density of 
the liquid is well known (to an accuracy of ~ 0.003% 
[1}). The piezometer volume and liquid density were 
used to determine the initial quantity of liquid present 
in the piezometer. Then the piezometer with its 
liquid contents was placed in a high-pressure thermos- 
tatic container, The volume of the piezometer was 
established with mercury, sealing the liquid in a thin 
capillary. On being heated the liquid expanded. Upon 
increasing the pressure to the required value and per- 
mitting part of the liquid to escape from the piezo- 
meter, it was possible to establish the initial level of 
the mercury which sealed in the liquid. By weighing 
the quantity of liquid escaping from the piezometer, 
one could find the density of the liquid under study. 
Appropriate corrections were introduced for the change 
of piezometer volume with temperature, and also for 
the change in the quantity of liquid in the piezometer 
capillaries located outside the thermostat, i.e., in the 
so-called parasitic volumes. 

Figure 1 shows the schematic experimental appa- 
ratus. The autoclave, consisting of a lower block and 
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an upper thick-walled flanged tube, contained the 
stainless-steel piezometer which has the shape of a 
cylinder with a cover at the top and a cone-shaped 
bottom. Two stainless-steel capillaries of 0.3 mm 
internal diameter are sealed into the bottom of the 
piezometer, One of the capillaries is connected 
through valve 5 to the water inlet system, and through 
valve 7 to the water outlet tube, The other capillary 
is sealed to a glass capillary (for recording the mercury 
level) through an adapter with platinum seals, The 
glass capillary dips into a glass of mercury which 
serves to seal the liquid under study in the piezometer. 
The platinum seals are located outside the thermostat, 
so that all the mercury is kept close to room tempera- 
ture. 

The heavy water is compressed in the piezometer 
by the nitrogen released into the autoclave from the 
heat intensifier system (gas cylinders M and B) through 
valve 8 and observation section 16, which serves to 
regulate the nitrogen input rate. The nitrogen entering 
the autoclave pushes mercury out of the glass 9 into 
the piezometer, compressing the heavy water to the 
required pressure. Thus,the piezometer is relieved of 
unilateral pressure, 


*This work was reported in 1957 at the research sym- 
posium of the All-Union Power Engineering Institute. 


The upper part of the autoclave contains the 
thermostat, which consists of a main copper block and 
two copper end blocks, provided with electric heaters 
which are connected to a voltage stabilized source of 
current, 

The volume of the piezometer, found by averag- 
ing twelve independent measurements, was determined 
to be 118.160 + 0.005 cm? at a temperature of 20° C, 
A second piezometer used in several of the experiments 
had a volume of 104.448 + 0.005 cm’, 

The volume of the piezometer which is not ther- 
mostatically controlled (parasitic volume) and con- 
tained D,O comprised ~ 0.12% of the above measured 
volume. The major portion of the parasitic volume 
was in the region which was at room temperature; the 
change in quantity of D,O within this volume under all 
experimental conditions was less than 1%, Thus the 
fact that this change was not taken into account did 
not introduce any errors> 0.001%, The parasitic volume 
in the zone of transition temperatures (from thermostat 
temperature to room temperature) comprised 0.02% 


of the piezometer volume and could safely be discoun- 
ted. 
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The temperature of the liquid under study, being 
in an equilibrium state, was taken to be that of the 
thermostat, and was measured with a platinum resistance 
thermometer to an accuracy of + 0.01°C, The pressure 
was measured with a standard spring manometer of 
class 0.2 and a piston manometer, connected to the 
autoclave via the relief reservoir and the oil pressure 
chamber. 

The main experiments were preceded by tests of 
the apparatus using H,O. The measurements of the 
density of H,O at 100° C gave results agreeing with 
known data to within 0.01%, The density of D,O was 
measured on.isotherms at about 50, 100, 150, 200, 

250 and 300° C with pressures up to 100 kg/cm*, All 
D,O density isotherms in the investigated pressure in- 
terval appear to be straight lines, as in the case of 
H,O. The experimental results showed that within the 
stated experimental accuracy the effect of pressure was 
the same on D,O and H,O densities, i.e., the ratio of 
these densities is independent of pressure within the 
investigated parameter interval. 

Figure 2 and the Table show the temperature 
dependence of the ratio of D,O and H,O densities. The 
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Fig. 1. Schematic drawing of experimental apparatus; 1) lower autoclave biock; 2) autoclave tube; 
3) flanged connection; 4) piezometer; 5) valve; 6) water inlet system to piezometer; 7) outlet valve; 
8) nitrogen inlet valve) 9) glass with mercury; 10) main thermostat block; 11) lower thermostat block; 
12) upper thermostat block; 13) spring manometer; 14) piston manometer; 15) relief reservoir; 16) sec- 
tion for observation of nitrogen input rate into autoclave; 17) oil pressure chamber; M) and B) gas 


cylinders of the heat intensifier. 
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value for H,O density is taken from [2] and [3]. An 
analysis of probable errors shows that the maximum 
error in D,O density measurements lies in the range 


0.01-0.03% (the latter value applies to the 300° C 
isotherm). The deviation of experimental points on 
the isotherms with respect to the smooth curves is of 
the order of 0.005%, this shows the good reproducibi- 
lity of the results and the smallness of the accidental 
errors. 
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0 60 100 10 20 20 t,t 
Fig. 2. Dependence of the D,O/H,O density’ ratio on 
temperature; @) data from present work; @) data from 
[1] (for P p,¢ at 25° C equal to 1.10448 g/cm’); O) data 
from [4]; 4) data from [5] ; —-— —) extrapolation of 
data from [6]. 


Dependence of the D,O/H,O Density Ratio on Temp- 
erature. 


t, °C t, °C 

50 1, 10904 180 1, 10862 

60 4 ,10930 190 1, 10834 
70 4,10944 200 1, 10802 

80 1,10950 210 1,10772 

90 4, 10953 220 1,10740 
100 4, 10955 230 1,10704 | 
110 1, 10955 240 1, 10660 
120 1,10954 250 1, 10610 
130 1, 10943 260 1, 10552 
140 1,10933 270. 110484 
150 1, 10920 280 1, 10405 
160 1, 10905 290 1, 10315 
170 1, 10855 300 4, 10240 


Figure 2 also compares our results with other work 
at 95-160° C [4] and 35-250° C[5]. The divergence 
of our data from the data of [4] does not exceed 0.01%, 
The results of [5] differ from ours significantly; the 
divergence increases with rising temperature and 
reaches 0.35% at 250°C, Such divergence can be 
explained by the low accuracy attained in the density 
measurements of [5]. Asforevaluation by the authors 
themselves, the measurement errors are on the order of 
from + 0.2% at room temperature to + 0.6% at high 
temperatures, Reference [6] similarly notes the unre- 
liability of the data of [5], and considers more correct 
the D,O/H,0 density ratios obtained by extrapolating 
the data of [4]. The data of [5] were taken into account 
during the extrapolations done in [6], and the density- 
ratio curve was drawn parallel-to a curve smoothing 
out these experimental data. It can be seen from Fig. 2 
that the extrapolation of the data in [6] fits closer to 
our data than does that of [5], but the divergence never- 
theless reaches 0.25% at 250°C, Thus our data does 
not confirm the extrapolated data from [6]. Figure 2 
does not show the experimental data of [7] and [8] 
(taken at temperatures ~ 100° C), which agree very 
well with our data, 

Figure 2 and the Table can be recommended for 
use in determining the density of D,O for pressures 
up to 100 kg/cm” and temperatures up to 300° C with 
the use of an HO density table. 
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Very few investigations have been made of the vis- 
cosity of heavy water. Of these, the work [1], in which 
the greatest range of temperatures is considered and 
which uses the method of Lavachek, does not appear to 
be particularly reliable, In [2], which uses the capil- 
lary method, the most reliable results are obtained, 
but the temperature range considered is relatively 
small — up to 125° C, 

The present work was carried out by using a me- 
thod developed by D. L. Timrot which had previously 
been applied in an investigation of the viscosity of 
H,O [3, 4]. The viscosity was measured by the capillary 
method, A difference in the pressure at the two ends 
of a capillary was created and measured by using a ring 
balance, which also served to measure the flow of the 
liquid through the capillary. 

A diagram of the experimental setup is shown in 
Fig. 1, The ring balance is formed from a tube of 30 
mm diameter, bent in the arc of a circle of 400 mm 
radius. The ring balance is joined by two tubes to a 
central cylindrical block, to which a central tube is 
attached by flanges; the capillary is located inside this 
tube. All the principal parts of the apparatus are 
made of stainless steel. The part of the central tube 
where the capillary is located is enclosed within a 
copper cylinder which is used in maintaining a constant 
temperature. On this cylinder are mounted three 
nichrome heaters; the main heater (along the whole 
length of the cylinder), and two guard heaters (to com- 
pensate for the loss of heat at the ends of the tube). . 
The heaters are supplied with a stabilized variable 
current, The uniformity of the temperature along the 
length of the experimental section is controlled by five 
thermocouples. The temperature is measured to an 
accuracy of +0.1° C, The constancy of the tempera- 
ture along the length of the capillary and the unifor- 
mity of the temperature during the experiment were 
maintained with a variation of less than +0.2° C for 
temperatures up to 100° C, and with a variation of less 
than £0.5° C for higher temperatures, 

Before the experiments the balance was compen- 
sated, i.e., the center of gravity of the balance was 
made to coincide with the axis of the rotating system, 
special regulating weights being used for this purpose. 
In order to determine the size of the correcting weight, 
an experiment was carried out with two different weights 


at the same temperature and pressure. This compen- 
sating weight did not exceed 5% of the main weight. 

Before starting the experiments, the apparatus was 
calibrated for a known pressure drop produced by the 
balance and a known rate of liquid flow. The correc- 
tion for the resistance of the junction tubes was deter- 
mined in special experiments carried out without the 
capillary, and did not exceed 0.5%, The size of the 
weight creating the drop in pressure was chosen so 
that the flow of liquid in the capillary was laminar, 
and with this condition the weight could vary from 400 
to 60 gm. 

The stainless steel capillary was 329.5 mm long 
with an internal and external diameter of 0.546 and 
2mm respectively. In order ‘to prevent any possibility 
of bending, it was located in a 6 mm diameter thick- 
walled tube of the same material. 

The coefficient of viscosity of the liquid was 
calculated from the experimental results by using 
Poiseuille’s equation with the Hagenbach-Cuette correc- 
tion for kinetic energy. This correction did not exceed 
1.5% of the coefficient of viscosity. In the calculation 
of the coefficient of viscosity the following corrections 
were taken into account; the correction x, for the 
natural travel of the balance, the correction x, for the 
resistance of the junction tubes, and the correction xg 
for the change in capillary dimensions resulting from 
variations in the temperature. 

The last correction was calculated by using the 
formula L; = Lo (1 + 0,151 t + 0.95 - 10°78 t? -0.5- 

10°" 

When the results of calibration and the corrections 
are used, the formula takes the form 


v= 0,03769 [G+ —467,8 4-25, 
PT; 


where v is the kinematic viscosity (m?/ sec), G, is the 
size of the main weight in grams, T is the time taken 
for the liquid to flow a distance of 1 mm on the scale 
(sec/mm),p the density of the liquid being investiga- 
ted (kg - sec?/m*), a =Yp/Yo the coefficient that takes 
into account the compressibility of D,O, y, the specific 
weight of D,O at a temperature of 20° C and at the 
pressure prevailing during the experiment, and y the 
specific weight of D,O at 20° C and under a pressure 

of 1 atm (the conditions under which the calibration of 
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Fig. 1. Diagram of the experimental apparatus; 1) Capillary; 2) ring tube of the balance; 3) balance 


junction tubes; 4) central tube; 5) cooler; 6) copper disc; 7) screen; 8) central cylindrical block; 9) frame; 
10) prism; 11) weight pan; 12) copper washer; 13) copper cylinder; Ty, Tz, Ts, T4, and Ts) thermocouple 


locations, 


TABLE 1, Experimental Values of the Viscosity of D,O 


P, tT 
kg/em? kgesec*/ | Gu 8 sec nm 
m4 
15,2 46,0 113,0 400 8,84 200 17,94 130 147 117 1,256 
16,0 47,5 113,0 100 36 ,88 50 80,00 128 144 114 1,263 
16,5 47,9 113,0 400 8,50 200 17,52 123 139 113 1,230 
18,1 97,2 113,2 200 15,68 100 30,27 122 138 109 1,266 
18,3 97,5 113,2 400 7,78 200 15,46 117 132 107 1,234 
18,3 97,5 113,2 200 15,48 100 30,32 118 134 107 1,252 
18,8 99,0 113,2 200 15,32 100 29,86 118 134 106 1,264 
19,1 59,0 113,0 400 7,88 200 16,18 115 130 105 1,238 
19,2 m1,5 113,41 400 7,62 200 15,17 114 - 130 105 1,238 
20,6 57,0 113,0 400 7,38 - 200 14,70 110 125 101 1,238 
47,4 70,8 113,41 200 8,45 100 16,92 62,4 70,6 59,0 1,197 
56,0 75,8 | 112,0 200 7,69 100 15,97 | 54,9 61,5 54,3 4,199 
68,3 89,7 141,2 200 5,95 100 11,43 45,5 50,6 42,7 1,185 
96,9 100 ,0 109,2 100 8,34 50 16,26 31,4 34,3 30,0 1,143 
98,7 103,14 | 4109,4 100 8,54 50 17,11 | 31,0 33,8 29,4 1,150 
101.0 120,0 | 409,0 100 8,75 50 18,23 | 30,9 33,7 28,7 41,174 
109.2 100,0 | 108,2 45 4,92 75 9,16 | 29,0 31,4 26,9 4,167 
148 ,7 86,0 104,3 80 7,10 40 13,53 21,5 22,4 19,4 1,155 
201.2 114,2 98 3 80 5,38 40) 10,10 | 16,0 15,7 14,0 4,124 
205.3 242,,0 998 80 6,14 40 13,20 | 16,4 15,9 14,1 1,128 
206, 4 53,6 97,14 80 6,00 40 12,94 15,6 15,2 13,6 1,118 
246 ,7 321,8 93,7 80 5,44 40 11,30 14,4 13,5 12,1 1,116 
247,2 325,0 93,6 80 5,43 40 11,30 14,4 13,5 12,1 1,116 
288 ,0 111,5 84,7 60 5,85 30 10,74 13,2 11,2 10,1 1,109 
*Experiments with a smaller weight G2 were carried out, but only for the calculation of the correction x. 
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TABLE 2. Smoothed Values of the Viscosity of D,O for Saturation Pressure 


/m? 


HD,0+10°, kgesec/ 


kg-sec/ 
/m? *H20 


30 193,3 
40 80,7 
50 67,1 
60 56,9 
70 48,8 
80 42,4 
90 37,4 
100 33,4 
110 30,4 
120 27,8 
130 25,4 


1,258 
1,247 
1,228 
1,212 
1,199 
1,188 
1,179 
1,172 
1,165 
1,159 
1,153 
1,148 
1,143 
1,139 


160) 20,0 1,133 
170 18,8 4,130 
180 17,6 1,128 
190 16,5 4,125 
200 15,6 1,123 
210 14,9 1,122 
220 14,2 1,419 
230 13,6 1,118 
240 13,4 4.117 
250 12,5 1,115 
260 12,0 1,413 
270 11,5 1,110 


0 50 


were taken from [5]. 


be the same as in the case of H,O. 


100 


the balance was carried out). The densities p of D,O 


The apparatus described was used to measure the 
viscosity of D,O (100%) at temperatures varying from 
room temperature up to 283° C, and for pressures from 
45 to 325 kg/ cm’, The experiments showed that the 
influence of pressure on the viscosity of D,O, in the 
temperature and pressure range used, could be taken to 


In Table 1 are given the experimental values we 
obtained for the viscosity of D,O and in Table 2 are 
given the smoothed values of the viscosity of D,O for 
saturation pressure. An analysis of the experimental 


150 


Fig. 2. Experimental results obtained by various authors for the viscosity of D,O. QO) results of the 
present work; @) results of [2]; A) results of [1]; X) results of [7]. 


260 250 


errors showed that the maximum error in our measure- 
ments amounted to +2%, It is evident from the 
results that the ratio of the viscosity of D,O to that of 
H,O decreases with increasing temperature; thus it is 
equal to 1,24 at 20° C and equal to 1.11 at 280°C. 


[6]. 


The viscosity of H,O was obtained from tables in 


In Fig. 2 a comparison is given of our results with 


the results of other investigations. The values of the 
ratio Hp io! HH,0 obtained in [2] agree with our results 
with an accuracy of up to 1%, The values in [1] are less 
than ours by 2% on the average, while the divergence 
increases for increasing temperature. The results given 


in [7] are also shown in Fig. 2. 


15 140,5 
20 127,7 
: 280 | 11,4 1,109 
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The heat capacity of, heavy water has been in- 
vestigated by means of an adiabatic circulating calo- 
rimeter with a closed circulation loop.t A diagram of 
the experimental apparatus is shown in Fig. 1. A gear- 
type circulation pump forces D,O into the flow meter 
from which it enters the inter-tube space of the heat 
exchanger; it is heated in this system to the desired 
temperature and then enters the measurement calori- 
meter which is located in a thermostatically controlled 
bath, From the calorimeter the water enters the heat 
exchanger tube, where it is partially cooled by an 
opposing flow, and then enters a cooler from which it 
returns to the circulation pump. The pressure in the 
device is due to nitrogen which is admitted from a 
chamber in a system with a mercury seal. By means 
of the latter is is possible to obtain uniform pressure 
of the liquid in the apparatus and the gas in the high- 
pressure container in which the circulating pump and 
its associated apparatus and drive (asynchronous motor) 
are placed; this eliminates any pressure asymmetry in 
the pump. 

The use of a circulating pump with a drive which 
is driven by an asynchronous motor makes it possible 
to obtain a constant flow of D,O in the circulating 
system (variations of 0,1-0,2%) without special flow 
regulators, The measurement calorimeter is in the 
form of a bifilar coil, at the end of which thermocouple 
housings are soldered. Inside the coil, on a quartz 
frame, there is a nichrome heater which is supplied 
from a high-capacity storage battery. Two copper 
cylinders are mounted on the frame and serve to equa- 
lize the temperature; in the gap between these cylinders 
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there is a multi-junction differential thermocouple 
thermometer which surrounds the entire surface of the 
calorimeter. Thus, the thermometer readings are deter- 
mined by the integrated heat exchange of the calo- 
rimeter with the surrounding medium. The calorimeter 
is made adiabatic by regulation of the thermostat so 
that the differential thermometer readings are equal 

to (or close to) zero. The thermometer also serves 

for measuring the very small uncompensated heat flow 
through the shell of the calorimeter. The increase in 
the temperature of the D,O in the calorimeter is mea- 
sured by a differential chrome-alumel thermocouple 
which is located in the measurement housings. This 
thermocouple is also used to measure the temperature 
of the water which enters the calorimeter. 

The flow of liquid is measured by a calorimetric 
technique which is essentially the following. The 
adiabatic calorimeter is inserted into the circulation 
system and the D,O is admitted at a temperature close 
to room temperature. Inasmuch as the heat capacity 
of D,O is fairly well-known at this temperature, while 
the effect of pressure on the heat capacity is small 


* This work was reported in 1957 at the Scientific- Tech- 
nical Conference of the All-Union Heat Engineering 
Institute, 


t This method is used at the All-Union Heat Engineering 
Institute for studying the heat capacity of H,O and 
steam at high pressures and temperatures [1, 2]. In 
carrying out the present work a new experimental appa- 
ratus was built; this system is considerably different 
from any carrier apparatus. 


and can be taken into account, the flow of D,O can be alumel thermocouple. In order to determine the flow 


determined by a measurement of the amount of heat of water,use is made of reliable experimental data on 
which enters the water in the calorimeter and the the heat capacity of D,O at atmospheric pressure in 
corresponding increase in the temperature ot the water the temperature range 10-50° C [3]. The effect of 

in the calorimeter. The flow-meter calorimeter is pressure (0.4% at 20° C and 100 kg/ cm?) is taken into 
set up in the same way as the measurement calorimeter. account using the data for HO. The pressure is mea- 
In order to reduce to a minimum the exchange of heat sured by a class 0.2 standard spring manometer. 
between the calorimeter and the surrounding medium Before the experiments are carried out the diffe- 
the latter is covered by an insulating material and placed _rential thermocouples are calibrated (by means of a 
in a Dewar. The increase in the water temperature standard platinum resistance thermometer) as are the 


in the flow meter is measured by a five-junctionchrome- thermometer, flow meter, and measurement calori- 


Fig. 1. Diagram of the experimental apparatus; 1) adiabatic measurement calorimeter; 2) thermostat; 
3) heat exchanger; 4) flow-meter calorimeter, 5) refrigerator; 6) gear-type circulation pump; 7) mercury 
trap; 8) U-tube; 9) chamber; 10) standard spring manometer; 11) nitrogen chamber; 12) control valve; 
13 -20) high pressure valves; 21) high-pressure chamber. 
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50 100 150 200 
Fig. 2, Experimental data on D,O heat capacity: x) data obtained 
(for 50 kg/ cm’); O) data obtained (for 100 kg/ cm’); A) data obtained 
[3] (for 1 kg/cm’), 
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Relationship in Actual Values of D,0 
Heat Capacities. 


CpD20 
CpH20 


Cpp,o 4t 
100 kg/cm" 


Cpp,o at 
56 kg/cm® 


1,002 


meter. The adjustment and checking of the apparatus 
are carried out in experiments with H,O, the physical 
properties of which have been studied with a high 
degree of accuracy over a wide range of pressure and 
temperature. The values of the heat capacity of H,O 
have been measured at temperatures ranging from 60 
to 250° C and pressures of 50 kg/cm? and found to be 
in agreement with the data of [4] to within an accuracy 
of 0.2-0.3%, The heat capacity of D,O has been 
studied at the isobars 50 and 100 kg/cm? and tempera- 
tures from 50° C to a temperature close to saturation 
at the corresponding pressure. 

The temperature of the D,O admitted to the 
measurement calorimeter is maintained constant to 
within an accuracy of 0.01-0.02° C, The difference 
in temperature at the input and output of the calorime- 
ter varies from 4 to 8° C,while the amount of heat 
generated by the heater of the measurement calorimeter 
varies from 40 to 70 kcal/hr; the correction for heat 
exchange between the measurement calorimeter and 
the surrounding medium varies from 0.01 to 0.1 kcal/hr. 
The flow of D,O through the calorimeter is approxi- 
mately 8 kg/hr. 

In analyzing the experimental data we have 
computed the ratio of the heat capacity of D,O at the 
mean temperature in the measurement calorimeter 


and the given pressure to the heat capacity of D,O at 
the mean temperature in the flow meter at the same 
pressure: 


(Cy)p, t —(Qu+ qu) Atp 
(Cp)p, ty (Qp+9p) Aty’ 


_ where Q, and Q, are correspondingly the amount of 


heat generated by the heaters of the calorimeter and 

the flow meter; qy and q, are the corrections for heat 
exchange between the calorimeter and the surrounding 
medium (as determined from the thermometer readings); 
At,, and At, are the differences in temperature between 
the input and output coils of the calorimeter. This 
ratio is then reduced to the heat capacity of D,O at 

20° C by introducing corrections (less than 0.2%) from 
the data of [3] for the departure of the mean tempera- 
ture of the water in the flow meter from 20° C, 

The experimental results are shownin Fig.2. It 
is apparent from this figure that the spread in the ex- 
perimental points with respect to the smoothed curves 
is less than 0.2%, 

Starting with the experimental data of [3] and 
assuming that the effect of pressure on the heat capacity 
of D,O at 20° C is the same as for H,O (0.2% at 50 kg/ 
/ cm’) we have computed the true heat capacity of 
D,O (Table). We then determine the ratio of the heat 
capacities of D,O and H,0O at pressures of 50 and 100 
kg/cm’ as a function of temperature (Fig. 3). 

The ratio Cyp,0/ CpH,0 is found to be indepen- 


dent of pressure (over the range which has been inves- 
tigated), This ratio decreases as the temperature in- 
creases (from 1.006 at 20° C to 0.956 at 250° C) and 
reaches a minimum close to 250° C (or remains essen- 
tially constant up to 300° C), Fig. 3 can be used for 
computing the heat capacity of D,O from the data for 
the heat capacity for H,O. An analysis shows that the 
maximum error of the present measurements of heat 
capacity is less than 0.8%, 

In Fig. 2 we compare the results of the present 
investigation with the data of [3](10-50° C), It is 
apparent from this figure that the temperature behavior 
of the heat capacity indicated by the data of [3] is 


ne 


50 


150 200 20 6,0 


Fig. 3. Relationship of D,O and H,O heat capacities. 


20 1, 1,000 
40 1 0,998 0,996 
60 4; 0,996 0,993 
80 0, 0,996 0,995 ee 
100 q 0,995 0,992 
120 0, 0,996 0,992 
140 0,980 0,998 0,994 
160 0,972 1,004 1,000 
180 0,966 1,014 1,008 
200 0,962 1,030 1,023 
220 0,958 1,051 1,043 
240 0,957 1,083 1,073 
; 260 0,957 1,130 1,116 
280 0,957 — 1,178 
| 
930 


in very good agreement with the temperature behavior 
found by us, 


The values of the heat capacity determined in [5] 
at temperatures of 10-125° C are found to be lower 
than the present values, The discrepancy is 1-1.2%, 
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Up to the present time the thermal conductivity 
of heavy water has been investigated only over a very 
narrow temperature range. The first experimental data 
on the thermal conductivity of D,O are reported in [1] 
and these data refer to the temperature range 10-60° C. 
In [2] results are reported for measurements of thermal 
conductivity at temperatures ranging from 2°C (super- 
cooled liquid) to 80°C. The experiments indicate that 
the thermal conductivity of D,O increases with increas- 
ing temper ature and that this increase is smaller than 
the increase in the thermal conductivity of H,O. The 
present work has been undertaken to investigate the 
thermal conductivity of D,O in a wide temperature 
range, up to temperatures close to the critical temper- 
ature. This data is necessary in order to carry out cal- 
culations involving thermal effects in nuclear reactors 
in which D,O is used. 

In making these measurements of the thermal 
conductivity coefficient X we have used a hot-wire 
method such as that which was used earlier for making 
relative measurements [3]. In the present work this 
technique has been used for absolute measurements, 
The basic element of the experimental apparatus is 
a quartz measurement tube (Fig. 1) which is placed in 
an autoclave made from stainless steel. Along the axis 
of the tube is stretched a thin platinum wire which 
serves simultaneously as a heater and as a resistance 
thermometer. The measurement tube and the volume 
inside the autoclave are filled with the liquid being 
studied. The heat flow from the heater through the 


liquid being investigated is in the radial direction, 
toward the wall of the tube. On the outer surface of the 
tube a platinum resistance thermometer is mounted, 

By means of surrounding heaters along the length of the 
measurement section it is possible to set up isothermal 
conditions in the autoclave. All conductors from the 
measurement tube located inside the autoclave are 
insulated from each other by quartz capillaries, which 
provide good electrical insulation, In order to eliminate 
convection flow of heat,the thickness of the layer of the 
liquid being investigated 6 is made rather small. The 
inner diameter of the tube D = 1.044 mm, the diameter 
of the platinum wire d = 0.096 mm, and 6 = 0.474 mm, 
The length of the measurement section (the distance 
between the voltage leads) 1 = 79 mm, The outer 
diameter of the tube is 2,35 mm. The platinum heater 
is centered inside the measurement tube by short quartz 
capillaries which are located outside the measurement 
portion. The centering is checked by means of a 
special optical bench which is provided with micros- 
copes, The average eccentricity of the platinum heater 
is 0.02 mm. Calculations show that if the quantity A 

is determined by the use of the formula for coaxial 
cylinders, 


QulAt 


and eccentricity is neglected the error in X is only 
0.07%, 


A, kceal/m-hr-C° 


U,d5 


0,50 


S 
> 


= 


035 


Fig. 1. Diagram 

of the quartz mea- 
surement tube; 1, 2) 
current leads for the 
internal resistance 
thermometer (hea- 
ter); 3, 4) voltage 
leads for the inter- 
nal resistance ther- 
mometer; 5, 6) 
rent leads for the 
external resistance 
thermometer; 7, 8) 
voltage leads for the 
external resistance 
thermometer; 9) 
quartz capillaries 
for centering the 
wire; 10) tungsten 
spring. 


Because the electrical conductivity of water increa- 
ses with increasing temperature it is necessary to 
determine whether or not this effect distorts the results 


of the measurements since the heater and the resistance 
thermometer at the wall of the tube are placed in the 
liquid being studied. In order to investigate this prob- 
lem we have carried out measurements of X in water 
with different degrees of electrical conductivity (by 
dissolving NaCl in the water), These experiments 

were carried out at 20°C, The results obtained for the 
values of X in these experiments for ordinary water of 
different electrical conductivity, up to 6 10™ ohm™*. 
cm * (this value corresponds to the electrical conduc- 
tivity of water att ~ 360° C), show essentially no diffe- 
rence (within 0.5%) from the values of A for distilled 
water, Thus, these experiments show that the design 
of the present device, based on the use of a hot wire, 

is completely suitabl e (from the point of view of elec- 
trical conductivity) for measuring X in water up to tem- 
peratures close to the critical temperature and, in any 
case, up to t ~ 360° C, 


In this work we have used pure D,O, The appara- 
tus is filled with D,O which is degassed in vacuum. All 
the electrical measurements are made with a PPTV-1 
potentiometer, The pressure is measured with a stan- 
dard manometer. The platinum heater, which serves 
simultaneously as a resistance thermometer, is kept 
under pressure by means of a spring; for this reason it is 
calibrated against the resistance thermometer which is 
located at the outside wall of the measurement tube. 
This calibration is checked many times during the ex- 
periment. In the majority of cases the discrepancies 
with the initial calibration are less than 0.02° C. Agreat 
effort is made to eliminate convection, Generalization 
from experiments on heat exchange in a gap between 
coaxial cylinders shows that at values GrPr < 1000 the 
transfer of heat takes place only by virtue of thermal 
conduction [4, 5]. These conclusions are based on ex- 
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Fig. 2, Experimental data on the thermal conductivity of D,O along 
the saturation line. 
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TABLE 1, Experimental Data on the Thermal Conductivity of D,O Along the Saturation Line 


Temperature} Tempera- 
difference at) ture differ- 
the wall Pa ence in the Reniarks 
2 fhe tube At, | layer 
Pkg/cm °C Atyig °C Q, keal/hr °C 
1 1,16 7,87 0,8282 0,506 25,3 
,86 5, ’ ’ ’ 
4 76 ‘ easurements c out, 
4 0,88 5,90 0,6281 0,512 28,8 } after the experiment t=220°C 
| oe | | | | camels 
1°19 7,92 08461. 0,514 Measurements carried out 
4 4,19 7,94 0,8470 0,513 34,3 after the experiment t=355°C 
4 4,20 7,92 0,8772 0,533 64,1 
15 1,03 6,62 0,7510 0,546 85,9 
15 1,38 8,94 41,0104 ,044 88,2 
45 1,25 8,16 0,9269 0,546 118,7 
17 1,26 8,50 0,9522 0,539 1494 
33 1,26 8,92 0,9732 0,524 4185,7 
33 1,26 9,25 0,9812 0,509 204,0 Pressure correction 0,2% 
33 0,79 6,13 0,6258 0,490 220,4 
33 1,24 9,62 0, 9835 0,491 221,8 
37 0,79 6,64 0,6364 0,463 243,6 
37 0,94 7,88 0,7537 0,460 246,9 , 
74 0,54 4,55 0, 4187 0,441 264 ,4 Pressure correction 0,4% 
67 0,52 4,83 0, 4290 0,427 279,14 
114 0,54 5,32 0,446 299, 
129 0,54 5,36 0,4343 07385 | 310/3 | Pressure correction 1.0% 
133 0,44 5,04 0,3734 0,355 325,3 Pressure correction 0, 47% 
207 0,27 3,40 0, 2393 0,319 344,3 Pressure correction 5.8% 
208 0,11 1,52 0,0950 0,282 355, 2 | Pressure correction 6,0' 


TABLE 2, Smoothed Values for the Thermal Conductivity of D,O Along the Saturation Line 


| A, A, A, 
cal/mehr- kcal/mehr> kcal/mehr- kcal/mehr- 
0,492 
20 0,504 120 0,547 220 0,493 320 0,365 
30 0,510 130 0,545 230 0,483 330 0,346 
40 0,518 140 0,543 240 0,472 340 0,325 
50 0,525 450 0,540 250 0,464 350 0,300 
60 0,532 160 0,536 260 0,449 360 0,265 
70 0,538 170 0,530 270 0,437 — pete 
80 0,543 180 0,524 280 0,424 — _ 
90 0,546 190 0,518 290 0,414 — eg 
0,548 


periments with cylinders of relatively large diameters 
and relatively large gaps. More recent investigations 
indicate that when the gaps are small (with the cylin- 
ders in the vertical position) the value of GrPr lies inthe 
limits 1700-2000 [6, 7]. The values of GrPr in the 
present experiments (from 1 to 24 of Table 1) are less 
than 1000; in the remaining experiments the values are 
1070, 1130, 1460, 1880 and 1760. This means that it 

is not necessary to introduce corrections for natural 


convection; as indicated by the results in [5] this 
correction is less than 1%, In order to investigate the 
effect of pressure on AD 0 We have carried out addi- 


tional experiments at t = 100° C and P = 150 kg/cm’ 
to t= 225°C and P= 200 kg/ cm’, The results of 


the experiment show that the dependence of An.0 


Fe) 


pressure is the same as for H,O [8]. In determining 
Xp so we have introduced corrections for the loss of 


| 
| | 
: 


The values of \ obtained in [1] are somewhat higher 
(since these authors used water with 95% content of 


% DO). The data obtained in [2] are very close to the 
i a present data, with the exception of the experiment at 
E la t = 80° C (in the present experiments this is a limiting 
3 | Ss temperature), On the basis of the experimental data 
> sas Eel we have drawn up Table 2 for the smoothed-out data 


} for the coefficient X po on the saturation line in the 


045, temperature range 10-360° C. The dependence of 
AD pressure can be obtained from the results of 


the experiments with H,O (3, 8]. 


Fig. 3, Comparison of the experimental 
data on the thermal conductivity of 


D,O along the saturation line: O) present LITERATURE CITED 

data; A) data of [2]; X) data of [1] (95% 1, F, Meyer and M, Eigen, Z, Naturforsch.8a, 800 

D,O + 5% H,0). 500 (1953). 

2. A, Challoner and R. Powell, Proc, Roy. Soc. 238, 
heat from the ends of the platinum heater which, in 1212, 90 (1956). 
the present case, amounts to 0.2%, In all these experi- 3. D.L. Timrot and N, B, Vargaftik Zhur. Tekh, 
ments the pressure P was larger than the saturation Fiz, 13, 1063 (1940). 
pressure Pe All the experimental values of \ in 4, F. Kraussold, Forsch. Geb. Ingenieurwesens 5, 
Table 1 refer to the saturation line. In certain experi- 4, 186 (1934), 
ments corrections are introduced for the excess pressure 5. M, A, Mikheev, Fundamentals of Heat Transfer 
(P > P). [in Russian] (Gosénergoizdat, Moscow, 1956). 
The calculations show that the maximum error 6. E. Schmidt and W, Leidenfrost, Forsch, Geb. 

in the experimental data is less than 0.8%, The results Ingenierwesens 21, 176 (1955). 
of experiments for the saturation line are given inFig. 2. 7, R. V. Shingarev, Dissertation, All-Union Heat 
The deviations of the experimental points from the Engineering Institute (1952), 
average curve lie within the limits of 0.4%, The 8. Tables of the Thermodynamic Properties of Water 
maximum value of \ is found att = 110°C. In Fig. 3 and Water Vapor [in Russian] (Gosénergoizdat, 
we compare our results with the data of other authors, Moscow, 1958), 


ELECTROLYTIC SEPARATION OF RADIOACTIVE 
ISOTOPES OF RARE EARTH ELEMENTS 
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Electrolytic deposition of rare-earth elements is For the study of the electrolytic separation of rare- 
usually carried out from alcohol [1], acetone-alcohol earth elements we used co. found in equilibrium with 
[2] and other nonaqueous electrolytes [3]. We have Pr™ and Ho™, The purity of each isotope was verified 
previously worked out a quantitative method for elec- by measuring its half-life (T). The oxides were dissol- 
trolytic separation of actinide group elements (thorium, vedin hydrochloric, nitric, chloric, sulfuric, and oxalic 
uranium, neptunium, plutonium, americium) from acids, The concentration of Ce was 4 + 1074 pg/ml, 
weak. acid aqueous solutions [4]. The conditions found that of Ho was 10 © to 107° ug/ml; the carrier 
appeared to be a starting point for the present work on content was as high as 0,0025-0.04 mg/ml. The con- 
electro-deposition of radioactive isotopes of rare-earth centration of hydrogen ions was 107! to 1074 M; in the 
elements, 


presence of buffers, it was 10 ° to 10°° M. 
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Fig. 1. Dependence of electrolytic deposi- 
tion of cerium and holmium on pH of solu- 
tion, Current density 100 ma/ cm?, Dura- 
tion of electrolysis 2 hours. 
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0 2 50 7 100 125 150 175 200 225 250 
Current density, ma/cm 
Fig. 2, Dependence of electrolytic deposition of 
cerium and holmium on cathode current density. 
Hydrogen ion concentration 10-* M. Duration of 
electrolysis 2 hours. 


ce™ 


0 2 40 62 80 100 120 10 160 


Time, minutes 
Fig. 3. Dependence of electrolytic deposition 
of cerium and holmium on duration of the ex- 
perimental run, Hydrogen ion concentration 
1073 M, Current density 100 ma/ cm?, 


The electrolytic separation of rare-earth elements 
was carried out in the usual electrolytic cell consist- 
ing of a small platinum cup (anode) and a platinum 
petal 0 25 mm thick of elongated shape with 1 nf 
effective surface (cathode), In such a cell, in contrast 
to glass or plexiglas electrolyzers, one can operate 
with very small concentrations of radioisotopes and 
with any value of electrolyte pH without significant 


1’ 


> 
3 80 > 4, 
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0 0 20 30 40 50 60 70 60 90 100 
Time, minutes 
Fig. 4. Dependence or electrolytic deposition rate of 
holmium on pH of solution. Current density 100 ma/ 
/cm?*; pH: 1) 3.05; 2) 2.0; 3) 5.3; 4) 6.5; 5) 7.3; 


6) 8.1. 
100 
80 

40 

0 2 $60 60 80 100 120 10 160 


Time, minutes 
Fig. 5. Dependence of electrolytic deposition 
rate of cerium on pH of solution. Current density 
100 ma/cm?; pH: 1) 2.8; 2) 2.1; 3) 1.8; 4) 4.7; 
5) 6.5; 6) 8.0. 


Electrolytic Deposition (Ce™ + Pr) and Holmium 
From Various Acids. 


Yield, % 
Acid 

(Ce™ + pr™) 
HCl 98.7 100.2 
HC1O, 97.8 99.6 
HNO 100.6 97.4 
H2SO, Experiments not car- 

ried out 96.9 
Experiments not car- 

ried out 90.4 


loss of the radioactive element through adsorption by 
the walls, The solution was stirred by rotating the 
cathode at 80-100 rpm. The radioactive element 
deposited out on both sides of the petal-shaped cathode 
simultaneously during the electrolysis, The8 activity 
of the electrolytic deposit was measured with the usual 


B -counting equipment of type B-2, being careful to 
maintain uniform geometry. 
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In order to discover the optimum conditions for 
electrolytic separation of rare-earth elements, we 
investigated the dependence of the yield of each 
radioisotope on the cathode current density and on the 
electrolyte pH. We also studied the electrolysis kine- 
tics of these elements using various acids and various 
electrolyte pH values. 

It is known that rare-earth elements in aqueous 
solutions hydrolyze in the region of the cathode. There- 
fore it was most important to establish the influence 
of pH on the electrolytic process; pH of 1-4 was ob- 
tained by diluting the original solution with distilled 
water, while pH 5-8 was obtained by neutralizing the 
original solution with ammonia, Figure 1 shows the 
results obtained, As can be seen, electrolytic deposi- 
tion of holmium from the solution is quantitative for 
pH 2-5, while for cerium this occurs only for pH 2.5- 
-3,2, At higher pH values, cerium evidently hydro- 
lyzes, 

Figure 2 shows data on the dependence of rare- 
earth element yield on cathode current density. As 
can be seen, holmium deposits out quantitatively over 
a very wide range of current density (practically from 
25 to 200 ma/cm’), while cerium and praseodymium 
do so only in the range 100-150 ma/cm*, A current 
density of 100 ma/ cm? was selected as the optimum 
condition for electrolysis of the given isotopes. 

Figure 3 shows results of study on the kinetics of 
the electrolytic deposition of cerium and holmium, It 
is interesting to note that the electrolytic deposition of 
holmium on the cathode is practically completed in 
1.5-2 minutes, while it requires 1.5-2 hours for quan- 
titative deposition of cerium and praseodymium, We 
can hypothesize that on electrolysis of cerium an oxi- 
dation-reduction reaction occurs in the region of the 
cathode with the transition 


Ce (TV) Ce (Ill), 


which slows the process of its electro-deposition onto 
the cathode, The rate of holmium electro-deposition, 
on the other hand, is determined simply by ion mobi- 
lity. 


In connection with the large difference in the 
kinetics of electrolytic deposition of holmium and 


cerium, it is of interest to pursue the variation of their 
deposition rates with pH. Figurés 4 and 5 show the 
corresponding results for holmium and cerium, The 
curves obtained show that even for pH = 5,3 the greater 
part of the holmium (~ 80%) is found in the ionic 

state and deposits onto the cathode, while cerium hy- 
drolyzes almost completely (> 80%) for pH = 4.7. 

Thus the method described permits us to observe the 
behavior of a radioactive element in very dilute solu- 
tions, 

It is possible to use oxalic acid and any of the 
mineral acids (hydrochloric, chloric, nitric, sulfuric) 
as an electrolyte, since holmium and cerium electro- 
deposit out of them quantitatively under the optimum 
conditions; hydrogen ion concentration of 107° M; 
current density 100 ma/cm/’, duration of electrolysis 
1.5-2 hours, Reproducibility of the experiments under 
these conditions is + 3%, 

The Table shows the average results of three to 
five runs for various acids. This electrolytic deposi- 
tion method can be used not only for analytical chem- 
ical analysis of rare-earth elements, but also for 
preparation of their plate-deposits for study of isotopic 
composition. 
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PHASE DIAGRAM FOR THE NIOBIUM-RHENIUM 


SYSTEM 
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Niobium and rhenium are transitional metals 
with an unoccupied d-shell. Some comparative data 
on these two metals are given in the Table below. 
The following references may be consulted in 
the literature on the niobium-rhenium system: [1] 
reports on x-ray and microscopic investigations of 
cast alloys in this system. The formation of a x -phase 
(of the &-Mn type) with a = 9.67 A was established 
in alloys containing 60-82 at.wt. % rhenium, annealed 
for 72 hrs at 1200° C. In alloys close to 50 at.wt. % rhe- 
nium, annealed for 120 hrs at 1000° C, evidence was 
found for still another intermediate phase: the o-phase 
with a = 9.72, c = 5.07 A; c/a = 0.52. It was found 
that a solid solution of niobium and rhenium contains 
up to40 atwt.% rhenium, X-ray analysis of cast speci- 
mens of niobium-rhenium alloys revealed still only 
the x-phase, as reported in[2]. The lattice constant 
for the alloy with 75 at, wt.% rhenium was 9.67, A. The 
sites of niobium and rhenium atoms in the unit cell 
were verified on the basis of assumptions that the 
structural formula Reg,Nb7 and the coordinates of the 
atoms in the compound are equal to the coordinates 


Comparative Data for Niobium and 


enium. 
Characteristics Niobium Rhenium 
Atomic radius, A 1.45 1.37 


Lattice constants, 
As 


Body-centered 
cubic (a= 3,294) 


Hexagonal close- 
packed (a=2,755- 
-2.76; c=4.4493- 
-4,458:; c/a= 
= 1,614-1.615) 

Melting point, 

Specific weight, 
g/cm® 8.4 

Difference in 
atomic radii, 
compared to 
rhenium, in % 


2415 3170 


21.02 


+ 5,84 


of &-Mn for the ordered and disordered arrangements; 
the conclusion then followed that the ordered position 
might be considered validated by evidence, X-ray 
analysis of cast niobium-rhenium alloys in the alloy 
with 75 at. wt.% rhenium recently [3] revealed a x -phase 
with a = 9.638 A, and in the cast alloy having 50 atwt.% 
thenium, annealed at 1000° C, revealed a o-phase 
having a = 9.78, c = 5.11 A and a niobium-base solid 
solution. One author [4] takes the view that the o - 
-phase in the niobium-rhenium system is stable from 
the melting point to a temperature above 1200° C. 

Rhenium powder of purity 99.8 wt. % obtained by 
hydrogen reduction of KReQ,, and a forged niobium 
cermet whose composition (wt. %) was; 0.8 tantalum, 
0.04 titanium, 0.04 silicon, 0.02 iron, 2 + 10™ lead, 
2-104 bismuth, and the remainder niobium were 
used in preparing the alloys. 

A charge was prepared from lumps of niobium and 
rhenium powder, previously compacted and sintered 
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Fig. 1. Phase diagram of the niobium-rhenium 
system: O) single-phase structures; 4) two-phase 
structures; @) melting points, 
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Fig. 2. Microstructures of niobium-rhenium alloys; a) 53 wt.%Re, annealed 600 hr at 
1000° C, a-phase; b) 60.3 wt.%Re, annealed 50 hr at 1600° C, ( % +X)-phase; c) 71 
wt.%Re, annealed 50 hr at 1600°C, x + eut. (@ +X)-phase; d) 84wt.%Re, annealed 60 
hr at 1000° C, x-phase; e) 95wt.% Re, annealed 600 hr at 1000° C,(x +8)-phase; f) 
100 wt.%Re, as cast. 
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at 1700° C to ensure compactness and degassing nece- 
ssary in further smelting. The procedure followed in 
preparing the alloys is described in another article [4]. 

The melting point determination was performed 
by the “droplet method” after filling a hole drilled in 
the specimen with liquid metal. The temperature was 
determined by means of an optical pyrometer gradua- 
ted by reference to pure metals under the same ambient 
conditions. The procedure and arrangement resorted 
to are described elsewhere [5, 6]. 

Phase analysis of the alloys was carried out by the 
x-ray structural analysis and microscopic methods, also 
by measuring the hardness of the alloys and the micro- 
hardness of the structural components on cast alloys and 
alloys annealed at 2200° C for one hour, at 2000° C for 
5 hr, at 1500° C for 50 hr (in a high- vacuum furnace) 
and at 1000° C for 600 hr (with subsequent quenching in 
water), 

The x-ray structural analysis was performed using 
CuKq- and VKq- radiations in a RKD type diffraction 
camera. Polished thin sections were prepared by con- 
ventional grinding and polishing techniques, followed 
by etching [etchand; HNOs (conc) + HF (conc) 1: 1 
mixture], 

Hardness measurements were made on a Vickers 
instrument with 10 kg loading; a load of 5 kg was used 
for brittle alloys in the region of the chemical com- 
pound, Structural components were identified by 
measuring their microhardness on a PMT-3 instrument, 
under loads of 100 and 200 g. It is fitting to note that 
the hardness values for a load of 5 kg showed good 
agreement with the microchardness values for a load 
of 200 g. 


The phase diagram of the niobium-rhenium system 
obtained as a result of the research described is shown 
in Fig. 1. The microstructures of the alloys are rep- 
roduced in Fig. 2, 


The system features a large region of solid solu- 
tions of rhenium and niobium extending to roughly 
64 wt. %(or46,5 at.wt.%)rhenjum at the melting point 
and narrowing down as the temperature drops, to about 
56wt. %(or 38.5 atwt.Zorhenium at 1000°C, According 
to data secured by microscopic analysis, alloys con- 
taining up to 53 wt.Jorhenium possess a single-phase 
polyhedral structure for the solid solution, both as cast 
and as-annealed (Fig. 2a). The alloy with 60,3wt. % 
rhenium(42.6 at.wt. %) is a single-phase solid solution 
at temperatures close to the melting point, and breaks 
down into two phases at lower temperatures (Fig. 2b), 
whereupon the amount of the second phase present 
increases appreciably as the annealing temperature 
diminishes. With increased rhenium content, the 
melting points of alloys in that region drop from 2415° 
C for the pure niobium 2340° C for an alloy with 65 
wt.%orhenium, Primary crystals of the a-solid solu- 
tion form directly from the melt. 


There is a X-phase (of the a-Mn type) which 
forms out of the melt and the solid solution, based on 
thenium, at 2520° C, via the peritectic reaction L+8= 


=X. The homogeneity range of the x—phase hardly 
shows any variation with temperature, and compri- 


ses roughly 82-93 wt.% rhenium (or 69.5-86,.5atwt. %) 
The lattice parameter of the x - phase varies with the 
rhenium content, being 9.67 for the alloy with 83.4 

wt, %(71.5 atwt. %) rhenium and 96,1 A for the alloy with 
90.9 wt.% (82 at. wt.%) rhenium, which shows some appro- 
ximate agreement with reported data [1, 2]. The mic- 
rostructure of the alloy from the homogeneity region 

of the x -phase is reproduced in Fig. 2d. The structu- 
ral formula assigned to the x -phase corresponds to a 
compound RegNb7. Note that a dendritic intrinsic to 
X-phases *2 cast is only partially destroyed as a result 
of even very protracted annealing (600 hr at 1000° C) 
and this is related to the slow pace at which diffusion 
processes progress. The microhardness of the x - phase 
is 1080 kg/mm? under a loading of 200 g, Vickers 
hardness is 1050-1150 kg/mm? under a loading of 5 kg. 


A two-phase region & + x is to be found between 
the region of the x-phase and that of the niobium- 
base &-solid solutions, A eutectic forms in this region 
at 2340° C withrhenium content ~ 67-68 wt.%(50-51 
at. wt. %) (Fig. 2C). 


A highly interesting fact is the presence, in that 
region, of a o-phase, detected by x-ray analysis only 
in cast alloys having 64,5 wt.% rhenium together with 
the o-solid solution, and in alloys having 71.5 and 
74.3 wt,%orhenium accompanied by the x -phase and 
a small quantity of the o-solid solution. As already 
noted, there are contradictory views regarding the 
existence of the o-phase in the niobium-rhenium 
system. These disagreements possibly revolve about 
the varying degree of approximation to the equilib- 
rium state of the alloys studied by certain authors 
[1, 3]. The assertion stated in [3] to the effect that 
the o-phase is a high-temperature phase and exists 
in the range from the melting point to 1200° C cannot 
be taken as well grounded, In alloys annealed at 
2200, 2000, 1500, and 1000° C, we failed to detect 
the o-phase by any of the research techniques employ- 
ed, The alloys possess the two-phase & + xX structure 
in that region. A more thorough investigation of the 
stability region of the o-phase is being pursued. 

The single-phase region of the solid solutions of 
niobium and rhenium has been tentative ly defined. 

It is apparent that the solubility of niobium in rhenium 
at temperatures close to the melting point is not in 
excess of 2-5 wt.% niobium. 

A two-phase B + xX region lies between the 
region of the x -phase and the rhenium-base solid 
solution. The microstructure of the alloy with 95 wt.% 
(90.5 at. wt. %) rhenium, for that region, is seen in 
Fig. 2e. The melting point of the alloys in this region 
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rises from 2520° C (peritectic transformation L+8 = 
*=X) to 3170° C for pure rhenium (Fig. 2f). The pri- 
mary crystals of the rhenium- base solid solution form 
directly from the melts, The hardness of the alloys 
drops from 1100 kg/mm’ for the alloy having 90,9 wt.% 
rhenium to 300 kg/mm’ for pure rhenium. 
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It is well known that the yield of products of nuc- 
lear reactions from a thin target at a particular charged- 
particle energy Ep is consistently less than that obtained 
from a thick target. For example, in the Cu® (d, 2n) 
Zn® reaction at Ey = 10 Mev and energy loss AE, = 0.1 
Mev in the thin target, the ratio of the yields is AB/B™ 
= 0.05. In both cases, the decrease in the number of 


particles which is due to nuclear reactions is negligibly ~ 


small (1078-1074), and the chief reduction in their 
effectiveness in traversing the thick target is due to 
deceleration resulting in the particles being brought 
completely to rest within the target material, while 
they pass completely through a thin target, losing only 

a portion of their energy on the way. The magnetic 
field can direct these particles onto the target until 
their initial energy is exhausted. The effectiveness of 
these particles diminishes for the same reasons, although 
the specific activity increases sharply over that in a 
thick target. If, after each passage, an energy AE, 

is imparted to the particle, they will pass through the 
given target v times, executing nuclear transformations 
with a constant energy Ep = 4.8 10° (HR)? Z4/a. In 
that case, we are dealing with a kind of resonance bet- 
ween energy losses in the material and energy incre- 
ment, Similar conditions may be set up in cyclic acce- 
lerators by choosing a target thickness AJ, such that 
energy losses in the target will be equal to the energy 
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increment per single revolution, i.e., AE) = 2 ZgVo 
cos 9 (¢ being the phase of the accelerating field at 
radius R, Vp the accelerating potential difference ex- 
pressed in million electron volts). It is better to have 
two targets each of Al,/2 thickness, placed symmet- 
rically about the center of the accelerator. The ex- 
pression for specific losses gives us 
AE, 2Vo cos 


aZ 


The number of passages v will not be very large 
when these conditions are observed and when extremely 
thin targets are used with a narrow accelerating slit 
width compared to the length of a stable orbit. How- 
ever, under real conditions, AZ, values are not very 
small, and it becomes necessary therefore to take into 
account multiple scattering of particles in the targets, 
resulting in a broadened beam and, if the focusing of 
the magnetic field is poor, impinging of particles on 
the walls of the chamber. To evaluate this effect and 
the focusing properties of the magnetic field, we may 
use the results of Casseles et al. [1] and formula (1). 
We then have, for the lower bound of the effective 
number of passages and for the total yield Begg = [No/ 
/ A] off’ Al,, respectively, 


g/cm’. (1) 


In 


d\2 aZ 
— . 3 
Lh g2V 9 COS In 183Z—"/3 (2) 


iz 


0-4 (Hd)? Eqn (Eo) 6 (Eo) 


where o (Ep) is the reaction cross section in barns; H 
is the magnetic field intensity in kilo-oersteds; n(E9) 
is the fall-off index of the magnetic field; d is the 
semi-height of the dees, in centimeters. This rela- 
tionship shows that the total yield is a function of the 
energy of the incident particle in terms of the factor 
Eon(E9) (Eo), for which reason it is profitable to use 
that orbit of the particular accelerator where this pro- 
duct has its maximum value. It is perfectly obvious 
that any further increase in the inhomogeneity of the 
magnetic field and in the height of the dees in the 
region of the stable orbit will result in increased Vege, 
and consequently,in increased total yield. 

We learn from (1) that the resonant target thick- 
ness at a specified orbit depends on the phase, which 
in turn is dependent on the accelerator conditions for 
resonance, so that restraint (1) may be fully satisfied 
by adjustment of the magnetic field (e.g., with 
respect to neutron emission), choosing Al,< (Als) 9-0 


Note that although less energy is lost in a thin 
target than in a thick one, the heat-removal mechanism 
is less effective, so that observation of the anticipated 
effect at high currents becomes more difficult without 
further improvements in the target. A rough calcula- 
tion of the allowable currents yields 


where the first term corresponds to emission, the 
second to thermal conductivity; T is the absolute tem- 
perature of the target; and p are the thermal con- 
ductivities of the target in 104 Mev/cm - sec and the 


, (4) 
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THE RANGE OF 14.7-MEV PROTONS IN ALUMINUM 


density in g/cm’, respectively; S is the beam area; i 
is the initial current, 

For a cyclotron with H = 15.5 kilo-oersteds; Vy = 
= 0.1 Mev; n = 0.25; gy = 60° [2] and for a copper target 
we have: Alg = 2.2 mg/cm’; Vege = 24; Vers J = 200-300 
Ha and Boge /B = 1-2, i.e., the yield from the thin 
target is slightly larger than that from a thick target, 
while the specific activities differ by as much as a 
factor of hundreds, 

Increased yield of nuclear reactions where thin 
targets are used may be achieved even in the absence 
of a magnetic field, by positioning the targets one 
behind the other and offsetting energy loss in each tar- 
get by some resonance method. Here, the scattering 
effect will be less important, and thermal loads will 
be insignificant. 

The method proposed would be particularly suitable 
when using low-energy particles,e.g., in obtaining 
neutrons by D-D, D-T, etc. reactions, when the energy 
of the incident particles is not large (S 0.5 Mev), 
since 4B/B = 1 at low energies, and the increase in 
number of passages (or the number of targets, where no 
magnetic field is used) would provide a net gain, 
assuming the foregoing conditions are observed, 

In conclusion, the author would like to take this 
opportunity to express his acknowledgement to A. A, 
Plyutto, N. I. Leont’ev, and P, N. Dmitriev for valuable 
discussions of the results. 
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Aluminum is the most frequently used material 
for measuring proton energies by the absorption method, 
since it may be prepared either as thin foil or massive 
blocks to a high degree of homogeneity and purity. 
The range-energy relation [1] for protons in 
aluminum is based on the Livingston and Bethe curves 


for low energies, and Smith's calculations for energies 
upward of 15 Mev. In the 2-4 Mev energy region, this 
relation agrees well with Wilson's measurements, 
Experimental results from recent researches [2-4] 
and our own data confirm the validity of the theore- 
tical curve in the 14,6-18 Mev energy range, devia- 
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tions from predicted behavior remaining well within 
1-2%, 

Protons of 14,7 Mev energy were obtained by the 
He® (d, p) He“ reaction, at a 90° angle to the path of 
the beam of He* ions, The accelerating voltage was 
150 kv. The current at the target was maintained at 
100 pa, The target employed was heavy-water ice, 
Protons were recorded by two proportional counters 
hooked up in a coincidence circuit. The two counters 
were mounted side by side within a common brass 
housing, and were themselves brass cylinders (80 mm 
high, 27 mm in diameter) with side windows (10 mm 
diameter) for admission of the proton beam. The coun- 
ters were filled with an 8% methane-argon mixture, at 
38 mm Hg pressure. 

The protons traversed the two counters in success- 
ion, at right angles to their fibers. The entrance win- 
dow was hermetically sealed by an aluminum disk 
1090 p thick. A disk with a stack of aluminum foil 
sheets of equal thickness was placed in a vacuum tube 
between the target and counters, The proton flow was 
monitored by two independent proportional counters of 
similar design, inserted into the coincidence circuit, 

The results of the measurements, shown in the 
diagram, showed that the mean range of 14.7-Mev 
protons in aluminum is 325.74 1 mg/cm’, the theore- 
tical value is ~ 330 mg/cm*, The way to find the 
maximum value of the possible error becomes clear 
from the geometrical structure of the diagram, The 
energy of the protons is computed from tabulated values 
for the mass of the neutron, proton, He®, and He*, The 
spread in energies is accounted for by the fact that the 
window of the far counter was visible from the center 


60 
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310 320 330 340 350 384,65 
Absorber thickness, mg/cm 


Absorption curve for protons traversing aluminum. 


Range of Protons in the 14.6-18 Mev Energy Interval, 
in Aluminum, 


Proton energy, |Range,mg/cm? |References 


320 {2 
329,71 
343, 5+-1% [3] 
477 ,0---0,5 [4] 


Our data 
10,2 3 
18-40,02 


of the target at an angle ~ 2°, and does not exceed + 
+ 10 kev. Some uncertainty in the energy of the pro- 
tons is a consequence of bombardment of the target 
by an unresolved beam of (He*)*tand (He*)* ions. 

However, assuming that all observed protons owe 
their existence only to doubly charged ions, the maxi- 
mum proton energy should be 14.73 Mev, instead of 
the accepted value of 14.7 Mev. 

Energy losses by protons traversing the target may 
be neglected, since He® ions of 150 kev energy are 
slowed down by an extremely thin surface layer. 

Assuming that 1 cm of air is equivalent (in stop- 
ping power) to 1,5 mg/ cm’ of aluminum, the theoreti- 
cally expected straggling value should, in our case, 
be ~ 1.68%[1, 5]; the empirical value is ~ 1.54% 
(cf. figure). 

The Table provides experimental data from the 
literature for the range of protons of 14.6-18 Mev 
energy, in aluminum, 

The authors gladly avail themselves of this oppor- 
tunity to thank O, I. Leipunskii for posing the problem 
involved and discussing the work, P. A. Yampol'skii 
for much valuable advice, V. I. Smirnov for his parti- 
cipation in the development of the associated electro- 
nic equipment, and L, A, Pets for his kind assistance 
in operating the accelerating tube. 
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Specimens of fissionable substances enclosed in 
cadmium containers were irradiated for 5 min in a 
stream of fast neutrons originating in the D (d, n) He® 
and T (d, n) He* reactions on a cascade generator, and 
were transferred to a counter. The geometry of the 
experiment and the mode of operation of the counting 
equipment and neutron source were the same as in 
previous experiments on relative- yield measurements 
[1], except for the manner of recording the counting 
rate, The total count of delayed neutrons Ny over a 
specified time interval per unit count of the monitor 
was recorded in this experiment. 

A fission chamber with a layer of the fissionable 
substance was then placed in the position corresponding 
to the position of the irradiated specimen, and the 
number of fission events N¢ in the same neutron flux 
was measured. As a result, Ng/ Ng ratios were obtained, 
to be used later in calculating the ratios of absolute 
yields of delayed neutrons of different elements fission- 
ed by neutrons of different energies. The specimens 
used in the experiments were in the form of disks 35 mm 
in diameter, 1 mm thick, and planar fission chambers 
with a layer of the fissionable material 35 mm in dia- 
meter, and thickness larger than the maximum range of 
any fission fragment, were used, 

It may be readily shown that the ratio of absolute 
yields of delayed neutrons from the same element 


fissioned by neutrons of different energies is, for a thin 
sample: 


Wal 


SEy 


a 


Here, indices Ey and E, refer to the corresponding 

neutron energies; |! is the absolute yield of delayed 
neutrons; ap is a correction taking into account anisot- 
ropic distribution of fission fragments averaged over the 
corresponding spectrum of neutrons responsible for fission, 
where a = R~ 1 (and R is the ratio of the yield of frag- 
ments at angles 0° and 90° with respect to the path of 


the neutrons); SE, and SE, are corrections accounting 
for attenuation of flux and double encounters within 
the specimen; is an abbreviation for the summation 
ty te 

( , in which Cj represents 
the relative, empirically measured, yield of the ith 
group of delayed neutrons (2 Cj = 1); Tj is the corres- 
ponding mean lifetime; t, and tz are the beginning and 
end of the measurement interval in time (it being 
assumed that activation of the specimen is brought to 
the saturation point). 

The — pe absolute yields of delayed neutrons 
from U™ and » Viz., s/Hg, when those substances 


are fissioned by monoenergetic neutrons, is found from 
the equation 


x [1- 9 Bs 
Of Bs (2) 


subscripts 8 and (5 + 8) referring to U™® and to 90% 
enriched U™, X3/ X5+) is the ratio of counting effi- 
ciencies for the first group of delayed neutrons when 
the respective specimens are transferred to the center 
of the counter unit (the change in counting efficiency 


‘is related to self-absorption); M denotes the mass of 


the specimen; and of®) 
tions of U™ and u™, 

The ratio xg/X(54g) is determined empirically by 
placing photoneutron sources Na + D,O and Na + Be, 
in the form of disks of the same diameter as the speci- 
mens, at the center of the counter, together with the 
specimens of and 

The ratio of absolute yields of U™* and Th™ upon 
fissioning of those nuclides by monoenergetic neutrons 
is determined from the relation 


are the fission cross sec- 


. 


where subscripts 8 and 2 refer to u™ and Th™, respec- 
tively. 

Ratios of absolute yields of delayed neutrons, for 
fission of U™ by thermal neutrons and by neutrons 
of 15,0 Mev energy, were also measured. In order to 
reduce the attenuation of the thermal neutron flux, a 
u* sample weighing 0.734 g was used in these experi- 
ments; the sample measured 35 mm in diameter and 
was placed in an aluminum container. A planar 
fission chamber made of aluminum witha layer of U™ 
(90% enriched), 1 mg/ cm? thickness, was used. 

With the ratio of absolute yields obtained from 
U™ fissioned by neutrons of different energies known, 
and the ratio of yields from U™® and U™ known for 
monoenergetic neutrons, we may find the absolute yields 
of delayed neutrons with the yield of delayed neutrons 
from fission of U™ by thermal neutrons taken as unity. 
The same holds for thorium. 

The values of the R coefficients are taken from the 
literature [ 2, 3], likewise the fission cross sections for 
energies 2,25, 3.8, and 15 Mev, which are taken from 
[4, 5, 6] respectively. Note that the median energy is 
indicated in [1] and in the present paper. The most 
probable energy value was,of course,used in the calcu- 
lations; instead of 2.4 Mev, 2.25 Mev; instead of 3.3 
Mev, 3.8 Mev (cf.,the spectra for the reactions involved 

The experimental results are entered in the Table. 
The tabulated data show that the absolute yields do 
not change as the energy of the neutrons responsible 
for fission increases to 4 Mev inclusive. Yields are 
almost doubled at neutron energies of 15.0 Mev. 


Absolute Yields of Delayed Neutrons, 


Neutron 
energy,Mev U238 


4,60+0, 28 
2,45--0,17 
2,58+-0, 14 


5, 1440, 38 
3, 18-£0, 26 
3,40-L0, 28 


The author would like to thank O. D, Kazachkov- 
skii and Yu. Ya. Stavisskii for their valuable counsel 
and kind attention to the work undertaken. 
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High-current accelerators presently in the design 
and construction stages are designed for the widest 
range of accelerated-particle energies. Interactions 
between the electron beam and various parts inside and 
outside the accelerator give rise to a powerfuly emis- 
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sion accompanied by intense neutron flux. The prin- 
cipal source of photoneutrons in the energy region 
extending to the region of meson production are the 
photonuclear reactions of types (y, n),(y, 2m), (7, Mp). 
The paramount contribution is accounted for by the 


15,0 1 ,86-+0,06 
3,3 0,99-+-0,04 
3 2,4 1,03+0,04 
i. 
2. 
3. 
4, 
5. 
6. 
** * 


(y, 1) reaction, which has a large cross section at 
energies of the order of 15-20 Mev. The beam is prac- 
tically completely absorbed as a result of the interaction 
between the electron beam on the one hand, and both 
accelerator parts and associated measuring equipment 
attached to the accelerator, on the other (e. g., inter- 
action with the Faraday cylinder used to measure 
electron-beam intensity). Photoneutron fluxes arising 
when the accelerators are operated at accelerated- 
electron energies in the 200-600 Mev range may 

attain values of the order of 10°-10" neutrons/ sec, 
according to data furnished in the literature [1, 2]. 
These figures back up the need to give special atten- 
tion to problems of shielding against neutrons during 
the operation of electron accelerators (this applies 
particularly to electron lineacs). 

In the high-energy region, approximation A from the 
avalanche theory [2, 3] may be found useful for esti- 
mating possible fluxes of photoneutrons originating in 
heavy elements during the operation of electron acce- 
lerators. The use of this approximation should not 
lead to large inaccuracies, since the basic restrictive 
inequalities used in the approximation are fulfilled, 
these being: Ey>> E, >> 8, where Ep is the energy of 
the accelerated energies; Ey is the resonant energy at 


which (y, n) reactions proceed effectively; 8 is the 
critical energy for the particular element. 

For elements in the medium and light range, the 
inequality E, >> 6 is not satisfied. Furthermore, 
electron accelerators are currently in wide use in 
various branches of science and engineering at particle 
energies of the order of 15-20 Mev, In this case, deter- 
mination of photoneutron yield requires the use of the 
more accurate results of avalanche theory, since B, 
e.g., for copper, is 16 Mev, Ep is 18 Mev, i.e., Eo is 
the same order as E,. In our work, calculations were 
extended to Pb, Cu, Al, and C, using the Belen'kii- 
Tamm equilibrium spectrum. The elements mentioned 
were used as the basic materials in working with 
accelerators (including Faraday cylinders, internals 
of the machine, shielding appurtenances). We also 
computed photoneutron yields from U and Bi, to corre- 
late the results obtained with results for elements of 
neighboring Z, 

Avalanche theory may prove useful in estimating 
the number of photoneutrons formed upon complete 
absorption of the photon and electron beams, In the 
case of a monokinetic electron beam, this quantity 
may be represented by the formula [1, 4): 


£0 
== de, (1) 


En 


where Q (Ep) is the number of photoneutrons arriving 
per single initial neutron, for photoneutrons formed 


Critical and Threshold Energies for Several Substances. 


Ref, for | 
Element | 3, Mev En, Mev relation 
Cc 68,7 18,7 [5] 
Al 36, 6* 14,0 [7] [6] 
Cu 16,4 8,3 (8) 
Pb 6,63* 7,3 [9] [12] 
Bi 6,54 7,4[5] {8} 
U 5,85 5,6 [10] {8} 
*Values of 8 for Pb and Al, cf. [4] also. 


upon complete absorption of a monokinetic electron 
beam of energy Eo; € = 2,29 E/B; Pp (Eo €)= 


=e= \ aa dx; J;— is the energy of photons in the 

shower; E,, is the threshold energy for the (y, n) reac- 
tion; Oy ,n (E) is the cross section of the (y, n) reac- 
tion; Op (E) = © pair ®) + 2Compton (E) is the total 
absorption cross section of the photons. Q(Ep) is cal- 
culated by integrating (1) using Simpson's formula; 
(E) is calculated from known formulas. Tabulated 
data are used, B values were found by the authors 
through the method outlined by Belen’kii [4]. The 
results of calculations. for Q(E9) are shown in Figs. 1 
and 2, 

In order to verify the numerical results, we ob- 
tained 0, , (E) for Cu, Al, and C by solving the 
integral equation (1) with respect to Oy sn {11}. Satis- 
factory agreement was found between the computed 
cross sections for (y, n) reactions and initial data, 
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no 10 15 20 25 


30 E,Mev 
Fig. 1. Dependence of photoneutron yields on energy 
of incident electrons for Cu, Pb, Bi, and U. Value 
per single fission: for Cu, 0.25 - 107% neutron/ electron; 
for Pb and Bi, 0.5 - 107° neutron/electron; for U, 107° 
neutron/ electron. 


| 
945 


Q(E,},neutron/ electron 


22 


22 Al EMev 
Fig. 2. Dependence of photoneutron yields on energy of 
incident electrons for Al and C. Value per single fission: 


for Al, 0.2 neutron/ electrons for C, 0.8 1075 
neutron/ electron, 


which confirmed the correctness of the program of 
calculations followed. 

The authors take this opportunity to express their 
appreciation to N. M. Balabanov, I. V. Gopak, G, A. 
Prokopenko, G, A, Redchenko for their kind assistance 
in performing the computational work. 
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Despite the large number of papers devoted to 
studies of interactions between neutrons and nuclei, 
the number of investigations in which the spectra of 
scattered neutrons are measured continue to be far too 
few. By studying the spectra of neutrons, it is possible 
to obtain unique determinations of energy levels and 
their partial excitation cross sections; a difficult task 
in measurements of gamma-ray spectra on account of 
cascade transitions, 
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In this paper, the authors present the results of 
measurements of spectra of neutrons scattered elasti- 
cally and inelastically on antimony and chromium nuc- 
lei, and scattered elastically on oxygen nuclei. 


The reaction D (d, n) He’ served as the neutron 
source, Deuterons were accelerated to 150 kev energy 
and directed onto a target of heavy ice. Scatterers 
were placed at a 13° angle to the direction of the deu- 
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teron beam, so that the energy of the neutrons was 2.94 
+ 0.1 Mev. 

The measurements were carried out in an annular 
geometry [1]. Scatterers of natural isotopic composition 
were chosen, each in the form of a torus of a round cross 
section, The chromium used was in oxide form, and 
placed inside a toroidal thin-walled container. Measure- 
ments were performed in four stages; the intensity of the 
direct beam was measured first, then the background, 
after which the intensity of the scattered neutrons was 
measured, and, finally, the background was determined 
again. The stability of the equipment performance was 
tested in each run of measurements by bursts of & par- 
ticles. The scattering angle was 90° in all cases. 

The neutron spectrometer used was a spherical 
ionization chamber combined with a linear accelerator 
and a 50-channel pulse-height analyzer [2]. In con- 
trast to [2], the chamber used to increase efficiency 
in [3] was filled with methane, which however had the 
effect of impairing energy resolution. Elastic and in- 
elastic scattering cross sections were computed in the 
manner indicated in [1]. 

The accompanying graphs show the neutron 
spectra obtained for neutrons scattered by antimony 
and chromium, Graph a depicts the spectrum of the 
direct beam; the dashed line shows the same spectra 
after introduction of corrections for wall effect [4]. 
Graph b gives the spectrum of neutrons scattered on 

antimony; graph c, the spectrum of neutrons scattered 
by chromium oxide; d, the same as c, but after account- 
ing for wall effect. 


OF LENINGRAD 


ATMOSPHERIC FALLOUT OF Sr® IN THE ENVIRONS 


It is clear from b that three excited levels of 
energies 0.85, 1.33, and 1.68 Mev are observed in 
antimony. The excitation cross sections of these levels 
are 22 2, 80 + 9, and 11 + 1 millibarnsAteradian, 
respectively. The elastic scattering cross section is 
147 + 16 millibarnsAteradian, 

In chromium (graphs c and d), a level at 1.82 
Mev is detected, in addition to the 1.45 level found in 
[5]. The corresponding differential inelastic scatter- 
ing cross sections are 59 + 6 and 24 + 3 millibarns/ 
/steradian, and the elastic scattering cross section 
is 208 + 28 millibarns/ steradian. 

Finally, the elastic scattering cross sections of 
neutrons at 2.9 Mev, for scattering on oxygen, is 47 + 
+ 5 millibarns/steradian, according to our data. 
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Determination of the amount of Sr” precipita- 
tion in atmospheric fallout in the environs of Leningrad 
have been carried out starting in 1954, The 1954-1955 
estimate was based on observations of the fallout from 


the total of nuclear explosion products [1]. Calculations 
were based on data from [2]. 
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Sampling of a large number of specimens of 
radioactive products descended from the atmosphere 
as fallout was carried out in 1956, It was assumed that 
the mean percentage of Sr” in those samples was the 
same as that in all active products contained in a year's 
fallout (the number of the latter is known from the data 


| 


Year's fall- ' Fallout resi- 


ut, mC/ d 
mC/km 


Total by Jan. 1, 1958 


in[1]). The amount of Sr contained in the 1956 total 
fallout can be estimated from it. 

In 1957, a radiochemical analysis of all of the 
fallout products trapped in a high-walled vessel with a 


OF NUCLEAR WEAPONS TESTING 
V. Santholcer 


Translated from Atomnaya Energiya,Vol. 7, No. 5, 


November, 1959 
Original article submitted May 25, 1959 


During 1956 the physics department of the medi- 
cal school at Hradce Kralove (Czechoslovakia), initially 
in collaboration with the Geophysics Institute and the 
Meteorology Division at Hradce Kralove, undertook 
systematic measurements of radioactive fallout. After 
the cessation of nuclear weapons testing, we also 
carried on sampling of radioactive fallout on a regular 
basis: at the present time, the schedule calls for 
sampling every two or three days, by a method descri- 
bed earlier [1]. A high-rimmed cuvette of plastic, with 
a collecting area of about 1 sq meter, is used in the 
sampling. A thin layer of water is maintained on the 
bottom of the vessel. The material collected is pro- 
cessed in the same manner as described in [1]. Ash is 
transferred while still moist to a low-rimmed alumi- 
num cup 3 cm? in area. After drying under a lamp, 
the cup and contents are placed in a conventional 
device for measuring 8 activity. Gamma activity is also 
measured for some samples, with a scintillation counter. 
The time behavior of the drop in activity of active 
specimens is monitored over several months, Excellent 
results are given by an automatic measuring device 


* 


RADIOACTIVE FALLOUT FOLLOWING THE CESSATION 


‘small amount of slightly acidic water on the bottom 
was carried out. Strontium and cerium were isolated, 
and the amounts of Sr” and Cc: ™ present were deter- 
mined. The amount of Ce contained in the 1957 fall- 
out was 55.7 mC/km?, 

Results of the Sr* analysis are presented in the 
Table. 
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recording the number of counts; the device was designed 
by J. Macek, a teaching fellow in our department. 


Samples collected from the 1 sq meter surface 
possess comparatively high weight (usually in tenths of 
a mg/ cm’), Self-absorption of 8 emissions must 
therefore be measured as accurately as possible, and 
over as large a population of samples as possible. 


On the basis of measurements of absorption in -alu- 
minum foil, it has been possible to approximately 
determine the energy of 8 radiation for our samples. 
This falls mainly into two energy groups: 0.6 and 2.1 
Mev. In samples of greater mass (which are predomi- 
nant), soft rays are practically absorbed entirely. 


The energy of y rays in radioactive fallout was 
determined approximately by scintillation counters, on 
the basis of absorption in copper and lead plates, and 
found to be 0.2 and 0.8 Mev. 


A Sr” standard obtained from the dosimetry section 
of the Institute for Nuclear Research in Prague (director: 
F, Behounek), served for absolute estimation of the 
activi ty of 8 -active samples. 
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Radioactive Fallout from November 1, 1958 to May 31, 
1959. 


Mean intensity | Mean intensi- 
for month, ty per day, 
mC/km? mC/km? 


November, 1958..... 8.1 0.27 
December, 1958..... 104,7 3.4 
January, 1959,...... 62.1 2.0 
February, 195)...... 75.1 2.7 
March, 1959. ...... 31.5 1.0 
April, 1959. ....... 45.8 1.5 
May, 1959.......6. 41.6 1,3 


The mean activity of atmospheric fallout for the 
current year (up to May 15, 1959) is approximately 
four times that for 1958, and approximately triple the 
mean activities recorded during the second half of 
1958. The highest activity (10.1 - 10~® curie /liter) 
was found in rain precipitated on April 6, 1959. 

Some data on radioactive fallout intensity are 
entered in the Table. 

Fallout precipitations were measured over several 
days, and the monthly activities of the fallout samples 
were derived from those measurements by simple 
summation processes, The fallout activity arrived at 
in this fashion for the entire half year was 327 mC/km?, 
The mean fallout activity for one day was 1.8 mC/km?, 
Low activities recorded for November,1958 are partially 
explainable by the low precipitation,meteorologically, 
during that month, 

For three months, from November 1, 1958 to 
January 30, 1959, the total amount of fallout activity 
was 175 mC/ km’; during the second three-month 
period, lasting up to June 30, 1959, it was 152 mC/ km?, 
However, instead of simple radioactivity sums, it 
would be more correct to determine the cumulative 
density of fallout activity [1]. By graphical methods, 
we found that the mean activity exponent of degree 
n = 1,13, or -1.64 for older samples, remained in force 
right up to the 3rd or 4th month of aging of the preci- 
pitations, For samples older than 100 days, we there- 
fore used the formula 


= A, (T—Ty)~ 114, 


for calculating the cumulative activity, A, being the 
activity per unit time; t = T-T the age of the sample. 
The data of nuclear testing Ty is determined 

approximately by graphical means (Figs. 1 and 2). In 
active samples, the activity is monitored constantly by 
measurements, so that the summands are true values, 
The same approach as described in [1] was applied to 
weak activity values. 
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From November 1, 1958 to January 31, 1959, the 
cumulative density of active precipitations came to 
90 mC/km’, the total from February 1 to June 30, 1959 
was 88 mC/km’, Since the radioactivity dropped from 
90 to 33 mC/km? in three months, the curnulative 
density one-half year after cessation of weapons testing 
was 121 mC/km’. 

Graphical determination of Ty and n conveniently 
pointed up the fact that the end points of the ordinate 
whose length was given by the difference in calendar 
dates, belonging to arbitrarily selected activities A 
and A/k,. define a straight line whose intersection with 
the time axis gives Ty while the angle formed by this 
straight line with the time axis is a function of k and 
the exponent n, 

It has been shown earlier [2] that the activity of 
radioactive fragments aged 100 days may be computed 
from the formula 


4/93 (Nov. 3, 1958) n=1,18 


‘ 


79731 10/31 11/30 12/31 1/31 2/28 
/ 1958 1959 


Fig. 1. Drop in activity of sample for radioactive 
fallout recorded, starting with November 3, 1958. The 
activity A (counts/ min) is plotted as ordinate, and the 
dates as abscissa. The date of origin of the activity 
(approximately October 11, 1958) is determined and 
confirmed graphically (curves 2 and 3), Straight line 
1 would correspond to the exponent n = 1.18, The 
exponent is found graphically from the triangle at 
left. Straight line 3 is the locus of difference in calen- 
dar dates (T ) associated with activities A and A/ 2. 
The departure of the empirical curve from the initial 
straight-line course of curves 2 and 3 is due to a com- 
paratively slight increment in activity contributed by 
previous nuclear weapons tests, 
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sample, whose activity at first dropped very rapidly) 
with n= 3-4, In March, 1959, the fall-off in radioacti- 
vity stabilized at n = 1.6, Gedeonov [4] indicated as 
he a likely explanation for the more rapid fall-off in 
activity the possible activity due to neutrons. Activity 
of such radioisotopes would diminish at a much faster 
rate than the activity of radioactive fragments. 
Samples of a certain minimum but sill apprecia- 
ble age remain rather active for a period of time 
a) following the cessation of nuclear testing. Even after 
several months had elapsed, most of the radioactive 
fallout samples were in a form such that they could be 
considered as isolated clouds of radioactive fragments. 
The relatively insignificant degree of mixing of atomic 
o clouds is pointed out in[5]. Time curves of decrease 
30 10/31 11/30 12/31 1/1 2/28 9/91 4/30 in activity for most of our samples could be taken as 
— — time curves for the fall-off in activity of radioactive 
fragments, whenever the increment in activity (a 
radioactive fallout, exponent n = 1,60, found by slight percentage attributable to initial activity) was 
gtaphical means. Sample 1 was taken on Jan, 23, at first negligible and appeared only after several 
1959, sample 2 on Feb, 16, 1959. The time of months had elapsed (cf. Fig. 1). 
origin of the activity was ascertained graphically 
to be in the neighborhood of September 25, 1959 LITERATURE CITED 
(see curves 5 and 6). Graphical controls are repre- 
sented by straight lines 3 and 4, 


~ 


Fig. 2. Fall-off of activity of two samples of 


1. V. P. Shvedov, V. A. Blinov, L. I, Gedeonov, and 
E. P, Ankudinov, Atomnaya Energiya 5, 5, 577 
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included a graphic representation of the curve of dimi- 3. J. Maly, J. Kutzendorfer, V. Machééek, V. 

nution in activity of B-active fragments, leading to n = Kourim, and V, Jetébek, Ceskoslovensky fysikalni 
= 1,5-1.7 for radioactive fragments three to four ¥asopis, 6, 661 (1956). 

months old, Our curves showing the time rate of 4, L. 1 Gedeonov, Atomnaya fnergiya 2, 3, 260 
decrease of active samples (18 samples),lend support (1957).* Bs 


for n = 1.64, for fragments more than three months old, 5. A. Sittkus, Beit. Phys. Atmosphare 30, 200 (1958). 
An anomalous.rapid increase in activity was ob- 


served only in a few samples (e.g., the Jan 23, 1959 
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ANALYSIS OF SOFT 8 RADIATION BY THE METHOD 
OF ABSORPTION IN AIR 


N. E. Tsvetaeva and L. A. Rozenfel'd 


Translated from Atomnaya Energiya, Vol. 1, No. 5, pp. 482-484 
November, 1959 
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Existing methods for analysis of 8 radiation by 
absorption in metal foils (of aluminum and gold) [1-4] 
are unsuitable for measuring the energy of soft 8 ra- 
diation (< 0.3 Mev) in a complex spectrum. In this 
case, good results may be had by using air as the absor- 
ber. 


Measurements of the absorption of 8 radiation in 
air were carried out in a vacuum desiccator into which 
the preparation and a T- 25 model BFL end-window 
counter with mica windows 1.3-2.5 mg/cm? were 
placed. The spacing between window and preparation 
could be varied from 1 to 11 cm. The pressure in the 
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desiccator was varied from atmospheric pressure to 20 
mm Hg with the object of varying the effective thick- 
ness of the layer of air. 

The absorption curves in air obtained for grag 
c™ (Ey = 0,155 Mev), s* = 0.167 Mev), Ca® (Ey = 
= 0.254 Mev), Co™ (Ep = 0.306 Mev), T1™ (Ey = 0.765 
Mev) display a strictly exponential shape if the spac- 
ing between preparation and counter is held to within 
8 cm. The ratio of the counting rate at zero to the 
counting rate at atmospheric pressure K = J, = 0/Jatm 


is a sensitive function of the energy of B radiation; the 
degree of sensitivity increases sharply as the length of 
the interval in air (Figs. 1, 2). 


The absorption curve for 6 radiation in air, for a 
two-component system and an air interval length =8 
cm, may be expressed in analytical form: 


_ Jsatm *1(Patm Featm ous 
(1) 


where Jys J atm are the counting rates at pressures p 


and atmospheric pressure, respectively,(subscripts 1 and 
2 referring to components 1 and 2); My, M2 are the 
absorption coefficients in air. 


=0,166 Mev, our 


106 
Rh 
3 E, Mev 


Fig. 1. Dependence of air coefficient K on 
energy of 8 radiation Eg: 1) distance between 
preparation and counter window h = 8 cm; 2) 
h=3 cm, 


By introducing the notation “atm = 2; 


Pe 
et (Patpyr P) — atm ; (PatnrP) 


we may transform (1) to the form 


Fatm—P PatnrP 


(2) 


Utilizing (2) as our initial equation, we derive a formu- 
la for determining the air coefficient for one of the 
components of B radiation in a two-component system: 


P atm 


(3) 


where pi, Pz, Ps are the values selected for the pressure, 


subject to the condition p2-py = Pz Joy Ino Ing are 
the counting rates at pressures pj, 


Taking into account the fact that x + y = 1, we 
derive from (2) an expression for determining the con- 
tent of a component; 


(4) 


ss 


K=J,_)/J atmos 


a. 


0712365678 MN 
h, cm 

Fig. 2. Dependence of air coefficients on 
the spacing between preparation and counter 
window h for different radioisotopes; 
(0.155 Mev); 2) s® (0. 167 Mev); 3) Ca® 
(0.254 Mev); 4) Co® (0.306 Mev); 5) T1™ 
(0.765 Mev). 
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Results of an Analysis of Mixtures Composed of Ca“ (Ey = 0.254 Mev; K = 3) and C“ (Ey = 0,155 Mev; K = 
= 10.7). 


Results of analysis 


ci4 Error in 


measure- 
Ca4s cu (Eo+ SE) Mae * vtsy (Eo+AE) Moe utdy ments, 


,% 


Content of components 


0,284 0, 60--0,03 


0,15179 | 0,32+0,05 4 


| 0,92+0,03 


| 1 


0,265+5 | 0,65-0,03 


0,147*) 0,270,05 1 


0,40-+-0 ,05 
0,52 0,48 


+908 0,50-40, 10 2 


eat positive and negative errors (AE) differ in absolute magnitude on account of the nonlinearity 
| of K= (Eo) (cf. Fig. 1). 


where K = Jp=0/ Jatm_ is the effective air coefficient Errors incurred in determining the air coefficient 


of the B radiation under analysis. (Ak) and the content of a component (Ay) become 
When studying a two-component system by means exact differentials of (3) and (4), where the variables 

of (3) and (4), it is well to assign a maximum length in the first case are J, Jz, and J3, and in the second case, 

to the air interval one at which the exponental charac- K and Kz. 


ter of the absorption curves would still prevail, i.e., 8 cm. 


1 
K,= +—atm x 
1\ Patm 
yK, "atm 


AK AK; (K—K)) 
(K;—K,)? 


Ay=+ 


To shed light on the possibilities inherent in the but at E, = 0.20 Mev and 10% content of the compo- 
approach outlined, (E2 + A E)/E,was plotted graphically nent to be determined. It is clear from Fig. 3 that the 


against the known energy value Ey and the energy value error in energy determination stays within 15% for the 
E, to be determined, at different parameters. 


Figure 3showsthe graphof (E,+ AE)/E, = f(E,) for energy ratio E,/E, = 1.3 and content = 50%, The error 
a different percentage content of the component to be diminishes rapidly as F,/E, increases, and at Ey/E, = 
determined (10, 20, 50, 100%) at E, = 1.4 Mev (solid = 1,7 the error is only 5%, The energy of a component 
curve). The dashed line indicates the same function, may be determined to within 10% accuracy at an energy 
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£,, Mev 
Fig. 4. Graph of(Ey+AE,)/E,=f (E,) at Ey = 0.14 
Mev, 


Fig, 3. Graph of (E, AE,)/E,= S(E,) at E, = const 
< Ey. 


ratio Ey/E, = 2, even for a small content (10-20%) of This approach has been verified and has demonstra- 
that component. ted excellent results in the analysis of mixtures compo- 
The conditions of the analysis are less favorable sed of Ca and C™ in different proportions (cf. Table), 
when the energy to be determined is higher than the and was successfully employed in the determination of 
known component, Figure 4shows that a low content of the purity of Zr” separated from Nb®, 
the component under analysis ( = 50%) combined with 
hardness is by a LITERATURE CITED 
ing error, which cannot be determined to sufficient 
accuracy by this approach at an energy > 0.4 Mev; E. Bleuler and W. Zunti, Helv. Phys, Akta 19, 375 
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lyzing 8 radiations by their absorption in air may be (1955). 
applied with greatest success to an analysis of the soft M, Forro, Z, Phys, 138, 3/4, 441 (1954). 
component at energies = 0.3 Mev, against a background 
of harder radiation, 


MEASUREMENT OF CONCENTRATION OF £ -ACTIVE 
GASES, USING END-WINDOW COUNTERS 
A. D. Turkin and L. M. Mikhailov 


Translated from Atomnaya Energiya, Vol. 7, No. 5, pp. 485-486 
November, 1959 
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End-window counters may also be used, in addi- centration of B -active gases in air. To find the 6 -ra- 
tion to already existing techniques, to measure the con- diation flux as a function of the concentration of the gas, 
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we considered the case where the end-window counter 
is placed at the center of the base of a cylindrical 
vessel filled with radioactive gas (Fig. 1). This arrange- 
ment, with certain assumptions, presents the possibility 
of determining by analytical means the relationship 
between 6 -radiation flux and the concentration of the 
gas. Account is taken of the absorption of 8 emissions 
of the radioactive gas in air and in the material of the 
entrance window of the counter. The following assum- 
ptions were made to find the relation between the B - 
particle flux through the entrance window of the coun- 
ter and the concentration of gas; 1) the radioactive gas 
is uniformly distributed over the entire cylindrical 
volume; 2) the absorption of 8 radiation proceeds 
exponentially; 3) scattering of 8 radiations back from 
the walls of the cylindrical volume and scattering in 
air do not occur; 4) the radii of the entrance window of 
the counter and of the cylindrical volume are in the 
ratio R/r = 5 [1] (this is a condition at which the 
counter may be regarded as a point-source), 

On the basis of the foregoing assumptions, the for- 
mula for the 6 -particle flux through the entrance 
window of the end-window counter is written as follows; 


(1) 


or 


(2) 


where q is the specific B activity of the gas in the 
measured volume (in curies per liter); S is the area of 
the counter entrance window (in cm); t_is the thickness 
of the counter entrance window (mg/ cm); lt is the mass 


Fig. 1. Geometry for deter- 
mination of 6 -radiation 
flux as a function of con- 
centration of gas. 


attenuation coefficient for B radiation (cm’/ mg); k is 
a conversion factor. 

Figure 2 shows values ofconversion factors k for five 
cylindrical volumes calculated for -radiation energies 
ranging 0.3-3 Mev, at a counter window thickness of 
5 mg/cm?, When any one of the indicated cylindrical 
volumes is used, the concentration of the 6 -active gas 
in air with a known energy peak in the 6 -emission 
spectrum is determined by measuring counting rate and 
solving the equation 


_ kN 
q= = 


(3) 


where 7) is a correction for the decay scheme of the 
measured gas[2]. The value of the factor k for a 
corresponding value of Eg is taken from Fig. 2. The 
value of N is determined in practice by taking two 
measurements of counting rate; the first measurement 
N, is carried out in a counter whose entrance window is 
screened by a light material 1-1.5 g/cm’ thick; the 
second measurement of N, is carried out using the same 
geometry minus the screen, In this second case, the 

B - particle flux N through the counter window is deter- 
mined as the difference N = No-Ny. Calculated values 
of k were verified empirically for 1.46 Mev energy of 

B radiation. In the experiments, a Sr” point-source 
was placed in the direction at right angles to the central 
axis of the end-window counter (a T-20 type counter , 
of window thickness 5 mg/cm’), The counting rate was 
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Fig. 2, Dependence of conversion factor 
k on peak energy in 6 -emission spectrum 
of radioactive gas: 1) H = 20 cmR = 10 
cm; 2) H = 30cm, R= 15 cm; 3) DZ-20 
chamber; 4) DZ-70 chamber; 5) H = 50 
cm, R= 25 cm, 
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measured as a function of distance between source and 
counter window, and as a function of the angle of 
incidence of the 6 -particle flux, This experiment was 
used to simulate cylindrical volumes uniformly filled 
with a 6 -emitting gas of known concentration, The 


accuracy of the measurements was 3%, Values of k obtain- 


ed by this method diverged from theoretically predic- 
ted values by as much as 10%, 
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Volume-distributed sources are more widely used 
in practice than point-sources, In the general case, 
the dose rate from a volume source will be a function 
of specific activity, specific composition of the y 
radiation, the source geometry, and the emitter mate- 
rial. For fairly large source dimensions, attenuation of 
y rays in the source leads to production of secondary 
y rays, Data on the yield of primary y -radiation from 
sources of the simplest geometries may be found in [1]. 
Our work involved an experimental investigation of 
the effect of multiple scattering in the source on dose- 
rate value, The volume sources used were a cylinder 
1.0 meter high and 1,5 meters in diameter, and a trun- 
cated cone described elsewhere [2]. The cylinder was 
filled with an aqueous solution of colloidal gold Au 
(Ey = 0.411 Mev) and a salt of Zn® (Ey = 1.12 Mev). 
The dose rate was measured with an air-equivalent 
ionization chamber placed along the axis at right 
angles to the gneratrix of the cylinder. The truncated 
cone was filled with an aqueous solution of salts of 
Na” (E, = 1.38 Mev; E, = 2.76 Mev), Zn® (E = 1,12 
Mev), and colloidal gold Au” (E = 0.411 Mev); the 
thickness of the active layer was varied from 2 to 100 
cm, and the dose rate was measured with the same 
ionization chamber. 


The dose rate for the truncated cone was computed 
from the formula 


2nK 
| 1 cos cos spy (il see ty) 


where Ky is the y constant; q is the specific activity; 
k is the attenuation coefficient in the source, Multiple 
scattering was accounted for by introducing into the 
formula a dose-rate build-up factor B (kl ; E): 


2nK.q 
xX [1—cos — (kl) tPo® (ki sec 


From experiments with the cylinder, we found that no 
relationship is detected between the B value and the 

detector-source separation for fairly large distances. 

The drawing shows the experimentally obtained curve 
for B (ki ; E) plotted against thickness of active layer. 
We see from insepction of the graph that; 1) B (kl :E) 
is a function of emitter energy and thickness of active 
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Build-up factor for scattered radiation in a volume 
source, as function of thickness of source and of energy 
of y quanta, 1) Na”, E, = 1.38 Mev, E, = 2.76 Mev; 
2) Zn, E = 1.12 Mev; 3) Au, E = 0,411 Mev. 
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layer, not exceeding 6 kv; 2) at k? > 6, B(kZ; E) is not where {Bj is the solid angle into which the source radi- 
dependent on the thickness of the active layer, but ates, 


solely on the energy of the y emitter; 3) B(kZ; E) varies 
inversely as the y radiation energy. LITERATURE CITED 

It is necessary, therefore, to take multiple scatter- y 
ing within the source into account in the calculations 1. G. V. Gorshkov, Gamma Radiation of Radioactive 
at kl < 6, while at ki > f the source may be viewed as Bodies (LGU, 1956). 
infinitely thick, with the dose rate found by the follow- 2. G,. V. Gorshkov and V. M, Kodyukov, Atomnaya 
ing formula: Energiya 5, 7, 71 (1958),* 
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ELECTROSTATIC METHOD FOR EXTRACTING ELECTRONS 
FROM A BETATRON 


A. A. Vorob'ev and B. A. Kononov 


Translated from Atomnaya Energiya, Vol. 7, No.5, pp. 487-488 
November, 1959 
Original article submitted February 19, 1958 


The electrostatic method for extracting electrons 
from a betatron [1-3] features several advantages over 
magnetic methods, The electric field may be mani- 
pulated with ease to effect appreciable deflection of 
electrons over a small interval of the electron path, 
The deflecting electric field may be localized to an 
extent such that the path of the accelerated electrons 
suffers no distortion before the electrons arrive at the 
deflecting condenser. No supplementary inductor wind- 
ings or high-power pulsed circuitry is needed then to 
broaden the orbit at the termination of an acceleration 
cycle, The height of the condenser may be made 
larger than the vertical dimension of the electron beam 
at the instant the orbit broadens at the end of the acce- 
leration cycle. 


The size and design of the deflecting condenser 
depend on the parameters of the betatron electron 
beam, specifically on the beam height and the pitch 

of its spiral trajectory as the orbit expands. Necessary 
information on the beam may be secured by means of 


Fig. 1. Result of photometric scanning of a 


probes of KCl or NaCl crystals placed at different crystal irradiated within betatron chamber; S) 
radial depths perpendicular to the beam trajectory. optical density of crystal darkening; x) distance 
The density of darkening of the crystal is proportional from edge of crystal (arrow pointing toward cen- 
to the number of bombarding electrons, making it rer of accelerator); y) vertical coordinate of 
possible to study the electron distribution in terms of beam ( y = 0 corresponds to the median plane 


the beam cross section (the optical density of the colo- of the accelerator). 
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ration is measured with a microphotometer). The spa- 
tial electron distribution of the 15 Mev Tomsk Poly- 
technic betatron is shown in Fig. 1. Beam height does 
not exceed 9 mm, and the number of electrons drops 
off rapidly as distance from the edge of the crystal 
increases, The loss of electrons on the plates of the 
deflecting condenser may be estimated by performing 
a graphical integration of the spatial figure formed, 
and the distance between plates may be assigned on 
that basis. The pitch of the spiral trajectory at radius 
R = 160 mm, corresponding to magnetic field index n = 
= 1, is 2mm, and 2.5 mm at R= 165 mm. When the 
range of values R= 160-165 mm is used, a cutoff plate 
0.05 mm thick corresponds to 16% electron losses, and 
the extraction efficiency is ~ 80%, 

The deflecting condenser for the 15 Mev betatron 
at Tomsk was designed as shown in Fig. 3, on the basis 
of the electron distribution over the beam cross section. 
The cutoff plate (50 y-thick tantalum foil), stretched 
between stainless steel guides, is grounded. Another 
plate of stainless steel is placed under high voltage. The 
edges of the plates are chamfered and precision ground 
to impart greater electrical strength to the interval. Op- 
erating voltage is kept under 40 kv, and is sufficient to 
deflect electrons of 15 Mev energy. The deflecting con- 
denser is mounted on a glass base. 

In choosing the voltage to be fed to the deflecting 
condenser, it should be kept in mind that the electric 
fringing field of the condenser is capable of perturbing 
the electron orbits at the moment of injection. This 
perturbation is so strong that when the injector is placed 
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Fig. 2. Curves of distribution of beam current I 
along coordinate x, necessary for determining 
the pitch of the spiral and electron losses on the 
condenser cut-off plates, for condenser set at 

R = 158 and R= 165 mm. Condenser electrode 
cross sections shown hatched. 


on the outside of the equilibrium. orbit, the beam-capture 
process is badly impaired at voltages across the condenser 
as low as~ 2 kv. By shielding the condenser properly, 
normal operation of the betatron is successfully engineered 
at deflector-plate voltages up to 6 kv. At higher voltages, 
say 9 kv, the number of electrons captured in the accele- 
ration process is negligible (1000 times less than nominal 
value). More involved shielding could of course result in 
satisfactory radiation intensity at the output at condenser 
voltages ranging in tens of kilovolts. Losses of accelera- 
ted electrons within the condenser would thereupon be 
increased, and extraction efficiency would deteriorate in 
the process. It is advisable to drop any attempt at using 
constant or variable voltage supplies to the condenser 

and to use pulsed supplies. 

There are no special requirements to be met by a 
pulsed voltage generator; it is sufficient that the pulse 
have the necessary amplitude and that the voltage re- 
main constant when the accelerated electrons arrive at 
the condenser. We investigated two pulse generators: in 
the first, the pulse-shaping circuit discharged into a pulse 
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Fig. 3, Design of deflecting condenser (dimen- 
sions given in millimeters), 


I, relative units 


Fig. 4. Electron current I as a func- 
tion of displacement of condenser 
through angle ¢ of optimized posi- 
tion, 
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transformer through a thyratron, and in the second, the 
circuit voltage was delivered to a step-up transformer 
with rapid-saturaring core (a high-voltage peak trans- 
former). If it is required to vary the energy of the ex- 
tracted electrons by releasing them at different points in 
the acceleration cycle, then the peak transformer should 
be supplied through a phase converter. This supply cir- 
cuitry is extremely simple and reliable, which is impor- 
tant when the betatron is to be used by personnel lacking 
any special technical qualification in that area. 

The equipment described is capable of extracting 
an electron beam from the 15 Mev Tomsk Polytechnic 
betatron. The beam parameters are as follows: mean 
current 4 + 107° amp (current per pulse 0,27 1a), cross 
section at 2 cm distance from discharge window (of 
aluminum foil 0.03 mm thick) is 6 X 10 mm, beam 
divergence over the first 17 cm of the path in air, in 
the vertical and horizontal planes, 5 and 8.5°, respective- 
ly, extraction efficiency ~ 60%, The 20% drop in 
efficiency below that indicated earlier is due to the fact 
that the cutoff plate of the condenser became slightly 
warped during rolling, Use of a higher quality foil 


would undoubtedly result in increased extraction effi- 
ciency. 

In the course of our research, we had the opportunity 
to observe the output of electrons in the absence of any 
voltage across the condenser, These electrons emerged 
in the form of a well-focused beam 6 mm across. The 
beam intensity depended on the position of the deflect- 
ing condenser (Fig. 4); the intensity dropped to zero when 
the condenser was removed from the operating area of 
the accelerator chamber. The observed electron output 
is explained as the interaction between the electron 
beam and the current induced in the plates. The extrac- 
tion efficiency was small ( ~ 6%), but still higher than, 
say, the efficiency attained when the beam is deflected 
by means of a magnetic shunt. 


LITERATURE CITED 


1. K. and H, Reich, Z, Phys, 126, 5, 383 (1949), 

2. HH, Pollock, Brit. J, Radiol. 26, 307, 368 (1953), 

3. A. A. Vorob‘ev, A. Kononov, V, A. Moskalev, 
and L. S. Sokolov, Izvestiya Tomsk. Politekhn. 
Inst, 82, 149 (1956). 


News of Science and Technology 


DIRECT CONVERSION OF NUCLEAR ENERGY TO ELECTRIC ENERGY 


An experiment on the conversion of thermal 
energy from a nuclear reactor to electric energy while 
bypassing the steam cycle was successfully completed 
for the first time at Los Alamos Scientific Laborato- 
ries, with the development of a thermionic converter 
("plasma thermocouple”) [1]. 

The thermionic converter is a vacuum diode in 
which one of the electrodes is heated by some suitable 
means, Thermoelectric current is due to the difference 
in the work functions of thermoelectrons coming from 
electrodes at different temperatures, and this current 
develops electric energy in its passage through the ex- 
ternal circuit between the two electrodes. 

The magnitude of the thermoelectronic current 
is limited by the space charge of the electrons in the 
interelectrode space. Several techniques have been 
suggested [2-4] for minimizing space-charge effects, 
from reducing the interelectrode space to 0.01 mm, 
which is impracticable, to much more promising appro- 
aches of setting up a plasma in the interelectrode 
space, 

The most efficient method for creating the plasma 
is that of contact (or resonance) ionization of a gas 
with a low ionization potential (e.g. ,Rb, Cs) introduced 
into the interelectrode space, Another method involves 
ionizing a gas occupying the interelectrode space by 
means of radioactive radiation. Both methods suffer 
from certain flaws. The first envisages the presence 
of surfaces heated to temperatures above 2000° C in 
the thermionic converter, while the second calls for 
a very large flux of radioactive particles. The effi- 
ciency of a thermionic converter, i.e., the conversion 
efficiency or ratio of electric energy produced to ther- 
mal energy expended, is calculated below on the assump- 
tion that heat losses include only losses by radiation, 
For the case of optimum choice of load resistance in 


the external circuit, and assuming no space charge, the — 


conversion efficiency is expressed as follows [5]: 
n=( ) 
[44 


where ®, is the work function [of the collector]; © is 
the Stefan-Boltzmann constant; € is the thermal emis- 
sivity; A is the thermal-emission constant in Richard- 
son's equation, and subscripts 1 and 2 refer, respectively, 
to the cathode, or emitter, and the anode, or collector. 


The efficiency of a thermionic converter may reach 


very high values, For example, for existing available 
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materials (having their corresponding values for , €, 
A), calculations show that 25-30% conversion efficiency 
may be had at an emitter temperature of 1100° C and 
collector temperature of 500° C. The total conversion 
efficiency may possibly be increased by utilizing heat 
taken from the collector. 

The conclusions of the thermionic converter 
theory were confirmed empirically by use of a special 
arrangement in which the cathode was heated by elec- 
tron bombardment [1, 6]. Figure 1 shows the emitter 
temperature dependence of the voltage across the elec- 
trodes of the thermionic converter with the external 
circuit open. Inspection of the graph demonstrates 
clearly that the open-circuit voltage developed is 
largely dependent on the cesium vapor pressure within 
the arrangement. 

A thermionic converter used for a demonstration 
of direct conversion of heat energy to electric energy 
is depicted in Fig. 2, The thermoelectron emitter in 
this case is a rod made of a solid solution ZrC ;: UC. 
Uranium in the lower portion of the rod is 94% enriched 


Voltage, open-circuit, v 


1 l i 

1600 2000 2400 
Emitter temperature, °K 

Fig. 1. Voltage as function of temperature with 
external circuit open, at different cesium vapor 
pressures, in a facility providing for electric 
heating of the emitter: 1) without cesium; 2) 
Pc, = 10" mm Hg; 3) Pcs = 10°? mm Hg; 4) Pcs= 
=4-+10°° mm Hg. 
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Fig. 4. Thermionic converter proposed by Martin 
for spaceship rockets [5]. This is comprised of a 
solid cylinder of uranium carbide in a graphite 
matrix with thoriated-tungsten cathodes and oxy- 
genated-tungsten anodes, concentric about the 
reflector. The anode functions as a radiator of 
excess heat, secondarily. This thermionic converter 
will weigh 560 kg and develop 27 kw electric po- 
m wer at 9.7% efficiency. The area of the cylindri- 

scale cal reactor core is 1800 cm”; cathode-anode spac- 

ing is 2.54 mm, 
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Fig. 2. Thermionic converter pow- 
ered by nuclear energy, operating in 
the Los Alamos reactor; 1) insula- 
tor material; 2) cooled collector; 

3) incandescent uranium fuel ele- 
ment; 4) plasma; 5) circulating oil 
coolant; 6) vacuum insulation; 7) 
liquid cesium pool. 


Fig. 5. Westingh 
conductor 8-9. stinghouse ther 


insulator moelectric fuel element 
hot junction [5]. This is a “Yankee” 

, type fuel element (uranium 
dioxide clad with stainless 
steel) enclosed in a layer 
of mixed-valence thermo- 
electric material, Alter- 
nating layers of p- and n- 
type material in the vertical 
direction form an array of 


elements whose emfs add in 
14,0 series. These elements con- 
f vert 5% of the heat from 
3,5 cladding of fission into electrical ener- 
fuel element gy; the remaining heat passes 
I / a into a conventional steam 
+ 
cycle, 
5 j doo a with U™ isotope. The collector is made of stainless 
3 : 5 steel, copper-edged. It is cooled by a circulating oil 
Boast Pe 41,5 flow. The space charge is attenuated by the cesium 
v rd 3 ions; cesium ionization occurs at the emitter electrode 
3 a x 41,0 a at temperatures 2000° C and up. The thermionic con- 

o - bs verter was placed in the core of the 5 Mw Omega- West 
pile. The open-circuit voltage and short-circuit cur- 
le rent values for the converter, as functions of reactor 

0 20 40 60 80 700 power, were determined during a reactor exercise 
% of peak reactor power (Fig. 3), Failure of a remote-control switch forestalled 


study of the functioning of the converter under condi- 
tions close to optimum power. However, by analogy 
with results of the experiment using electron bombard- 
ment to heat up the cathode, it was possible to verify 


Fig. 3. Voltage and current produced in a thermio- 
nic converter placed in a pile as functions of reactor 
power; 1) voltage; 2) current. 
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the fact that about 30 w power per 1 cm’ emitter sur- 
face area was attained, Calculation of the number of 
fissions occurring within the fuel rod made it possible 
to estimate the conversion efficiency. The estimate 
was 5%, The thermionic converter was operated 
through two cycles of 5 hr each, the voltage remain- 
ing unchanged with the external circuit open, and the 
current dropping 20% in the first hour with the circuit 
shorted. 

The use of a thermionic converter is not the only 
way to achieve direct conversion of heat energy into 
electrical energy. Semiconductor thermoelements, 
the theory of which has undergone extensive develop- 
ment [7], may be successfully employed. Semicon- 


ductor thermoelements cannot yield such high efficien- 


cies as thermionic converters, but they have their own 


advantages (such as greater reliability in operation and 
bypassing any need for complicated plasma- generating 


devices), At the present time, several designs have 


A 3 Mev LINEAR ELECTRON ACCELERATOR 


been proposed for thermionic and thermoelectric con- 
verters, Two of these are shown in Figs. 4 and 5. 


L. M. 
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O. Val'dner, P. Dmitrovskii, D. Zorin, Yu. Mizin 


A new 3 Mev lineac with a high current of acce- 
lerated electrons was commissioned in January,1958 


at the Moscow Engineering and Physics Institute. Some 


research experiments carried out on this machine dur- 
ing 1958 demonstrated that it is an effective instru- 
ment for solving various scientific and engineering 
problems associated with the use of a powerful beam 
of fast electrons, 


The accelerator may be used successfully in 
various fields in science and engineering: in experi- 
mental physics as a source of high-intensity ionizing 
radiations, in chemical engineering for irradiation of 
plastics to improve their properties, in metallurgy for 
nondestructive inspection of parts and activation ana- 
lysis of alloys and pure metals, in the food industry, 
and also in radiotherapy. 


The accelerator is powered by high-frequency 
power from a magnetron oscillator. The length of 
the accelerating waveguide is 1.18 meters, The elec- 
trons are injected into the guide at 50 kev energy. 
Electrons with total energy of 2.8 Mev at the peak of 
their spectrum have been obtained at the accelerator 
output, The halfwidth of the electron energy spectrum 
is ~7%, The energy is measured by a magnetic 
spectrum analyzer with a resolution of not less than 
2%, Future plans call for stepping up the energy of 
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the accelerated electrons to 3 Mev and higher,by 
increasing the high-frequency power fed into the 
accelerating waveguide, The current of accelerated 
electrons at the output is 230 ma per pulse, with an 
average current of 280 wa and pulse duration of 3 
Hsec, The pulse-repetition rate of the current may 
vary over a 200-400 cps range. The diameter of the 
output electron beam is less than 4 mm, 

Focusing of the electrons is achieved by means 
of a constant longitudinal magnetic field which is 
set up by a system of eight coils, and subject to 
control from 0 to 900 oersteds. The magnetic field 
is maintained constant along the length of the wave- 
guide, to within 5% accuracy. Almost 100% focusing 
of the electrons is effected at a field strength of 450 
oersteds, The power required by the focusing system 
at this point is 1.5 kw. The working vacuum is 10~° 
mm Hg. 

The accelerator functions satisfactorily over a 
frequency bandwidth of + 2.5 Mc. When the frequency 
varies greatly, the rate of energy decrease is 0.06 
Mev/ Mc on the average. All of the input parameters 
(frequency, high-frequency power, energy, injection 
current, etc.) are easily maintained constant during 
accelerator operation to within a range such that the 
output parameters of the electron beam may be kept 
constant, 
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Accelerator performance is remote-controlled 
with relative ease from a special control panel, 


BRIEF COMMUNICATIONS 


USSR, The main items of construction work are 
going up at a fast tempo at the site of the Novo-Voro- 
nezh nuclear electric power station, All of the com- 
munication lines have been laid and traffic has star- 
ted on a 20-km railway line and on the automobile 
highway. An electric power transmission line has 
been put into operation. A prefab-concrete plant, 
woodworking combine, concrete plant, machine shops, 
and several other units are functioning. 

The machine hall of the electric power station 
and the combined accessory works have been built, 
Concrete is being poured into the foundations of the 
main reactor building. 

USSR, The design of a steam turbine rated at 
70,000 kw, for the Novo-Voronezh nuclear power sta- 
tion, has been completed, and the first turbine unit 
is being manufactured at the Kirov factory in Khar*kov. 


The accelerator is designed for continuous opera- 
tion over anextended period of time. 


The turbines will operate on saturated steam at 29 
atmos pressure and a temperature of 231°C, Steam 
moisture will not exceed 12%, 

USSR, An exposition displaying the achievements 
of the USSR in the field of the peaceful uses of atomic 
energy was open to the public during September and 
October,1959 in Tashkent, The opening of the exposi- 
tion was scheduled to coincide with the inauguration 
of a conference on the peaceful uses of atomic ener- 
gy held at Tashkent in the early part of September 
1959, In the volume and scope of the various areas of 
application of atomic energy evidenced, the Tashkent 
exposition was the largest of any similar Soviet travel- 
ing exposition seen either in the USSR or abroad. 

Poland. Design of a reactor intended for research 
on radiation chemistry has been completed. 
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N. N. Shumilovskii and L. V. Mel'ttser, Fundamen- 
tals of the Theory of Automatic Control Devices Using 
Radioactive Isotopes. Moscow, USSR Academy of Sciences 
Press, 1959, 144 pp., 6 rubles, 40 kopeks. 

This book represents the first attempt on record to 
generalize the results obtained by different scientific 
research organizations in devising instruments whose 
operation is based on the utilization of radioactive iso- 
topes and nuclear radiations, for process control and 
automation, 

The book consists of an introduction and nine 
chapters. 

The introduction offers a concise exposition of the 
fundamental properties of radioactive emissions, a 
brief description of radiation detectors, and a discussion 
of basic measurement circuits, as well as several special 
problems associated with errors in measurements encoun- 
tered in the use of automatic control devices incorpo- 
rating radioactive isotopes. 

Eight chapters of the book deal with problems in 
the theory of measurement of thickness (density), relay 
devices, level measurement, liquid flowrate and gas 
flowrate measurement, gas pressure Measurement, and 
monitoring of the composition of matter. 

Each chapter provides a detailed consideration of 
the basic factors affecting errors in measurement, and 
derives formulas relating the characteristics of the ra- 
diation source to the instrumental constants and the pa- 
rameters of the object sensed, making it possible to 
evaluate the activity of the source, speed of response, 
and accuracy of measurement of the instrument, etc. 

The concluding chapter discusses some problems 
in determining minimum source activity in the instru- 
ment, as a function of the dynamic properties of 
the instrument. 

The material discussed is of considerable interest 
to engineering and technical workers engaged in design- 
ing and using automatic control devices based on the 
use of radioactive isotopes. 

However, it must be noted that the book would 
have still greater value if the reference material contain- 
ed within it were expanded to include thickness ranges 
of materials, densities and concentrations of various 
substances,as well as other parameters susceptible to 
measurement by sources of radioactive radiations, ab- 
sorption coefficients of various forms of radiation in 
matter, etc, 


Some of the subdivisions, e.g., "Automatic Gamma- 
Defectoscopy," don't do full justice to the work done 
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in this field. The authors make no mention of the find- 
ings reported by I, Fakidov at the Ural Branch of the 
USSR Academy of Sciences, and others in the field of 
the automation of gamma-ray flaw detection using 
scintillation counters, although these techniques are 
highly interesting. Nevertheless, despite the insufficien- 
cies noted here, the book is on the whole a highly valu- 
able and interesting contribution, which for the first 
time generalizes ten years of experience of work in the 
field of devising instruments of this type, and also in- 
cludes a detailed bibliography on related questions. 


bP. 


P, A. Andreev, A. A. Kanaev, and E. D, Fedorovich, 
Liquid-Metal Coolants for Nuclear Reactors, Leningrad, 
Sudpromgiz, 1959, 384 pp., 14 rubles, 35 kopeks. 

This book provides a survey of foreign and some 
Soviet work on studies of the properties of liquid metals, 
Chapters 1 and 2 of the first section are devoted to the 
physical and chemical properties of liquid metals and 
their alloys. Chapter 3 outlines the basics of fluid 
mechanics and heat transfer in liquid metals; chapter 4 
describes the corrosion properties of various structural 
materials of nuclear reactors used in liquid- metal en- 
vironments, 

The second section, "Design and Performance of 
Liquid-Metal Systems,” consists of three chapters, des- 
cribing the basic principles underlying the design of 
liquid-metal systems and their subassemblies, equipment 
used in the systems (pumps, heat exchangers, accesso- 
ries, measurement devices), and operating rules govern- 
ing liquid-metal systems. 

The third section of the book surveys Soviet research 
on liquid-metal reactor coolants. 

The book assembles and analyzes a large amount of 
theoretical and experimental material on liquid metals; 
it is of unquestionable interest to scientists, engineers 
and technicians engaged in the field of reactor design 
and related fields, 
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